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. Sumnury 

vThe  theoretical  method  presented  for  the  calculation  of  the  linearised  values  of  generalised  airforce 
coefficients  on  a harmonically  oscillating  flat  plate  wing  in  subsonic  flow  is  a further  development  of  the  lifting 
surface  theory  of  Multhopp  type,  and  has  been  justified  mathematically.  The  method  of  spanwise  integration  is 
refined,  and  the  number  of  upwash  points  taken  may  be  greater  than  the  number  of  basic  functions  used  in  the 
approximation  for  the  loading.  The  refined  method  of  spanwise  integration  leads  to  improved  accuracy  in  the 
evaluated  values  of  upwash,  especially  near  the  perimeter  of  the  planform.  By  taking  more  upwash  points  than 
the  number  of  basic  loading  functions  more  accurate  values  should  be  obtained  for  the  generalised  airforce 
coefficients  for  a given  number  of  basic  functions. 

The  procedure  has  been  programmed  in  ICL  1900  FORTRAN.  The  program  has  been  used  to  obtain 
results  for  tapered  swept  and  rectangular  wings.  For  the  particular  examples  of  tapered  swept  and  rectangular 
wings  the  convergence  of  the  results  is  studied.  /|^ 
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1.  Introduction 


Methods  for  calculating  generalised  airforce  coefficients  on  a harmonically  oscillating  flat  plate  wing  in 
subsonic  flow  based  on  linearised  theory,  have  been  in  the  course  of  development  for  many  years.  The 
linearised  theory  is  used  to  set  up  an  integral  equation  relating  the  unknown  loading  distribution  to  the  known 
upwash  distribution  on  the  wing.  Basically  two  methods  have  been  in  the  course  of  development  for  solving  the 
integral  equation  numerically.  In  the  one  the  loading  distribution  is  replaced  by  a distribution  of  concentrated 
loads  on  certain  lines  and  is  known  as  the  doublet  lattice  method,  whereas  in  the  other  the  loading  distribution 
is  replaced  by  an  approximation  which  is  continuous  over  the  wing  except  in  the  neighbourhood  of  its  leading 
edge  and  is  known  as  the  lifting  surface  method.  There  are  also  methods  which  are  not  based  on  the 
above-mentioned  integral  equation,  for  example  the  vortex  lattice  method.  In  this  paper  we  concern  ourselves 
exclusively  with  the  lifting  surface  method. 

The  lifting  surface  theory  of  Multhopp’  was  for  steady  flow  and  required  the  loading  to  be  approximated  to 
by  a polynomial  in  the  wing  coordinates  4 and  ri  multiplied  by  a function  of  f and  tj  which  took  into  account  the 
known  singular  behaviour  of  the  loading  at  the  edges  of  the  wing.  This  approximation  to  the  loading  was 
substituted  into  the  integral  equation  to  get  an  approximation  to  the  upwash,  and  this  approximation  to  the 
upwash  was  equated  to  the  known  upwash  at  a set  of  points  on  the  wing  equal  in  number  to  the  number  of 
unknown  coefficients  in  the  expression  for  the  approximation  to  the  loading.  A set  of  linear  equations  for  the 
unknown  coefficients  was  thus  obtained  and  this  set  could  be  solved.  The  approximation  to  the  loading  was 
then  known  and  an  approximation  to  any  required  generalised  airforce  coefficients  could  be  obtained  from  it. 

Multhopp’s  method  was  extended  to  low-frequency  harmonic  oscillations  in  Ref.  2 and  to  general- 
frequency  harmonic  oscillations  by,  among  others,  Acum*,  Richardson"*  and  Davies’.  It  was  found,  however, 
that  in  all  these  methods,  in  which  the  chordwise  integration  was  carried  out  first,  the  spanwise  integral  was 
evaluated  numerically  by  too  coarse  a method  and  this  resulted  in  inaccurate  estimation  of  the  approximation 
to  the  upwash  at  points  near  to  the  leading  and  trailing  edges  of  the  wing  with  consequent  loss  of  accuracy  in  the 
results  for  generalised  airforce  coefficients.  Garner  and  Fox*  refined  the  method  of  numerical  integration  of 
the  spanwise  integral  and  applied  their  refinement  to  the  case  of  low-frequency  harmonic  oscillations.  Long’ 
applied  this  same  refinement  to  the  case  of  general-frequency  harmonic  oscillations.  Zandbergen,  Labrujere 
and  Wouters*  also  refined  the  method  of  numerical  integration  of  the  spanwise  integral  in  a manner  somewhat 
different  from  that  of  Ref.  6 and  for  the  particular  case  of  steady  flow.  Lehrian  and  Garner’  then  extended  the 
refinement  of  Ref.  8 to  the  case  of  general-frequency  harmonic  oscillations.  A method  in  which  the  spanwise 
integral  is  evaluated  first  is  that  of  Hewitt  and  Kellaway*’.  The  numerical  integrations  in  Ref.  10  are  all  carried 
out  accurately  so  that  no  refinements  are  necessary. 

All  the  above  methods  require  that  the  approximation  to  the  upwash  be  equated  to  the  known  upwash  at  a 
set  of  points  on  the  wing,  their  number  being  the  same  as  that  of  the  unknown  coefficients  in  the  expression  for 
the  loading,  and  then  these  unknown  coefficients  may  be  determined.  There  are  other  methods  of  determining 
these  unknown  coefficients  and  one  of  these  methods  is  discussed  by  Davies".  In  this  method  the  coefficients 
are  determined  by  equating  integrals  involving  the  approximation  to  the  upwash  to  corresponding  integrals 
involving  the  known  upwash.  This  process,  theoretically,  leads  to  the  generalised  airforce  coefficients  being 
obtained  with  the  highest  possible  precision  for  a loading  approximation  of  a particular  form.  The  said 
integrals  are  evaluated  numerically  for  the  present  paper,  but  the  number  of  integration  points  may  exceed  the 
number  of  unknown  coefficients  in  the  expression  for  the  loading.  If  the  number  of  integration  points  equals 
the  number  of  unknown  coefficients  in  the  expression  for  the  loading  then  the  method  of  equating  the 
approximation  to  the  upwash  to  the  given  upwash  at  a set  of  points  on  the  wing  is  retrieved.  Furthermore  the 
refinement  of  Garner  and  Fox*  is  used  in  evaluating  the  spanwise  integral  in  the  integral  equation,  although 
this  is  modified  to  some  extent  in  that  the  parameter  q which  determines  the  number  of  spanwise  integration 
points  may  depend  on  the  location  of  the  upwash  point  concerned. 

The  process  is  described  in  detail  in  the  present  paper.  A program  has  been  written  in  ICL 1 900  FORTRAN 
to  calculate,  by  using  this  process,  the  generalised  airforce  coefficients  and  the  loading  distribution  for  a wing 
oscillating  harmonically  at  general  frequencies  in  subsonic  flow.  The  procedure  for  using  the  program  is 
described  in  Ref.  1 7 . Calculations,  using  the  program,  have  been  carried  out  for  swept  t-apered  and  rectangular 
wings.  The  results  obtained  are  given  here  and  their  convergence  is  studied. 

2.  Theoretical  Considerations 

2.1.  Preliminary  Formulae 

We  refer  points  of  space  to  an  inertial  right-handed  Cartesian  coordinate  system  Oxyz  where  O is  the  origin 
of  coordinates  and  Ox,  Oy,  Oz  are  the  axes  of  x,  y,  z coordinates,  positive  z being  upwards.  We  introduce  a 
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thin  wing  W into  the  space  and  consider  it  to  be  vibrating  in  such  a way  that  the  position  of  any  material  point 
on  the  surfaces  of  W is  near  to  a fixed  point  which  is  the  mean  position  of  that  surface  point.  We  consider  a flow 
of  fluid  in  the  space  about  the  wing  W,  which  at  large  distances  ahead  of  the  wing  is  uniform  and  horizontal  with 
speed  V in  the  direction  of  the  positive  x-axis.  We  then  have  the  problem  of  determining  the  flow  about  the 
vibrating  impervious  surfaces  of  the  wing  W and  we  take  this  to  be  potential  flow  except  across  the  wake 
surface  which  extends  as  a sheet  of  vorticity  downstream  from  the  trailing  edge  of  W.  Once  the  flow  is  known 
the  pressure  forces  acting  on  the  wing  can  be  determined. 

When  the  wing  W is  very  thin  and  is  vibrating  in  such  a manner  that  all  the  points  on  its  surfaces  are  always 
very  near  to  the  plane  z = 0,  the  governing  equations  of  motion  of  the  fluid  may  be  linearised.  Let  the 
orthogonal  projection  of  the  mean  positions  of  the  material  points  on  the  surface  of  W on  to  the  plane  z = 0 
define  the  area  S,  which  we  shall  call  the  wing  planform.  Then,  as  far  as  the  aerodynamic  problem  is  concerned, 
the  wing  is  replaced  by  the  area  S in  the  plane  z = 0 and  the  boundary  condition  of  the  fluid  not  penetrating  the 
surfaces  of  W is  replaced  by  given  fluid  speed  distributions  normal  to  the  top  and  the  bottom  surfaces  of  the 
planform  area  5.  The  wake  becomes  a flat  surface  in  the  plane  z = 0 extending  from  the  trailing  edge  of  S to 
infinity. 

There  are  two  material  points  on  the  surfaces  of  W,  one  on  the  top  and  one  on  the  bottom  surface,  whose 
mean  positions  have  an  orthogonal  projection  onto  the  point  (x,  y,  0)  of  the  planform  area  S.  Let  the 
components,  in  the  direction  of  the  positive  z-axis,  of  the  displacement  of  these  points  at  time  t from  the  point 
(x,  y,  0)  on  the  planform  area  S be  denoted  by  Z+(x,  y,  t)  and  Z-(x,  y,  t)  respectively  for  the  point  on  the  top 
surface  and  for  the  point  on  the  bottom  surface  of  W.  According  to  linearised  theory,  the  fluid  speed 
lV+(x,  y,  t)  normal  to  the  top  of  the  planform  area  S and  measured  positive  upwards  is  given  by 

W^x,  y,  0=  (1) 

Similarly  the  fluid  speed  W_(x,  y,  t)  normal  to  the  bottom  of  the  planform  area  S and  measured  positive 
upwards  is  given  by 


W-(x.y,  0 = 


^,aZ-(x,  y,  f)  ^ dZ-(x,  y,  t) 
dx  dt 


(2) 


The  perturbation  velocity  potential  in  the  fluid  can  be  split  up  into  the  sum  of  two  constituents,  one  of  which 
is  synunetric  about  the  plane  z = 0 and  the  other  of  which  is  antisymmetric  about  the  plane  z = 0.  The  speed 
distribution  normal  to  the  planform  area  5,  corresponding  to  the  symmetric  velocity  potential  constituent  is 
the  same  in  magnitude  but  opposite  in  sign  on  the  top  and  bottom  of  the  planform  area  S,  whereas  the  speed 
distribution  normal  to  the  planform  area  5,  corresponding  to  the  antisymmetric  velocity  potential  constituent 
is  the  same  in  magnitude  and  sign  on  the  top  and  bottom  of  the  planform  area  S.  Corresponding  to  the 
symmetric  velocity  potential  constituent  there  is  no  net  pressure  loading  across  the  planform  area  S but 
corresponding  to  the  antisymmetric  velocity  potential  constituent  there  is  a net  pressure  loading  across  the 
planform  area  S,  and,  by  the  principle  of  superposition,  this  is  the  total  pressure  loading  across  the  planform 
area  S.  This  pressure  loading  will  give  rise  to  generalised  airforces  on  S and  these  can  be  taken  to  be  the 
linearised  values  of  the  corresponding  generalised  airforces  acting  on  the  wing  W.  Accordingly,  to  determine 
these  generalised  airforces  we  need  deal  only  with  the  antisymmetric  velocity  potential  constituent. 

Let  us  now  write 


Z+(x,  y,  r)  = Z(x,  y,  t)+Zi(x,  y,  t)  (3) 

and 

Z-(x,  y,  /)  = Z(x,  y,  /)-Zi(x,  y,  f).  (4) 

The  function  Z|(x,  y,  t)  describes  the  wing  thickness  distribution,  which  will  not  normally  be  changing  with 
time,  whereas  the  function  Z(x,  y,  t)  describes  the  position  of  the  camber  surface,  which  will  be  changing  with 
time.  The  antisymmetric  velocity  potential  constituent  depends  exclusively  on  Z(x,  y,  t)  whereas  the  symmet- 
ric velocity  potential  constituent  depends  exclusively  on  Zi(x,  y,  r).  Since,  as  stated  above,  we  need  deal  only 
with  the  antisymmetric  velocity  potential  constituent  we  may  disregard  the  thickness  function  Z|(x.  y,  r) 


4 


i 


henceforth,  and  consider  only  the  displacement  constituent  function  Z(x,  y,  ()■  The  corresponding  upward 
component  of  the  fluid  velocity  normal  to  the  planform  area  S on  its  top  and  bottom  is  W(x,  y,  t)  where 

• ,5) 

dx  dt 

The  quantity  W{x,  y,  t)  is  called  the  upwash  on  S. 

For  a vibrating  wing  the  displacement  function  Z(x,  y,  t)  can  be  given  as  a linear  combination  of 
independent  modes  of  oscillation.  Thus  we  can  write 

Z(x,y,t)=l  I Mx,y)bdt)  (6) 

k-l 

where  / is  some  typical  length  of  the  planform  5,  (kix,  y)  is  the  modal  function  and  bk(t)  is  the  generalised 
coordinate  for  the  mode  number  k,  both  of  which  are  non-dimensional.  The  modal  functions  (k(x,  y), 
need  to  be  a complete  set  of  functions  for  (6)  to  be  valid,  in  general,  and  to  be  of  practical  value 
the  summation  in  (6)  must  be  truncated  to  a finite  number  of  terms.  This  truncation  may  entail  an  error  but  if 
the  number  of  terms  retained  is  sufficiently  large  the  resulting  error  is  negligibly  small. 

The  generalised  coordinate  bk{t)  is  a real  function  of  time  t,  which  we  shall  assume  to  consist  of  a linear 
superposition  of  harmonic  constituents 

bk(<o)e^  + bt(a>)c-^'  (7) 


over  a range  of  values  of  circular  frequency  a>,  where  btio))  is  the  complex  conjugate  of  the  complex  number 
bk((o).  Since  our  aerodynamic  problem  has  been  linearised  we  may  take 

bkit)  = bk(<o)e‘^  (8) 

to  carry  out  the  determination  of  the  generalised  airforces  at  the  circular  frequency  <o.  The  generalised 
airforces  for  the  problem  when 

6*(r)=Af(a,)e-“'  (9) 

are  the  complex  conjugates  of  those  for  bk{t)  given  by  (8).  The  generalised  airforces  corresponding  to  the  real 
function  bk(t)  given  by  the  expression  (7)  are  then  the  real  quantities  obtained  by  adding  the  two  complex 
conjugate  generalised  airforces  corresponding  to  (8)  and  (9). 

If  we  substitute  for  Z(x,  y,  t)  from  (6)  into  (5)  and  use  the  expression  (8)  for  bk(t)  we  get 


W(x,  y,t)=V  I at(x,  y;  v)bk(<o)e‘‘ 

k-l 

where  P is  the  number  of  terms  to  be  retained  in  (6), 

ak{x,  y;  v)=  I + w{k(x,  y) 


is  the  frequency  parameter  corresponding  to  the  circular  frequency  <u. 
Corresponding  to  the  upwash 

Va*(x,y;.«)**(«)e‘-' 
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in  the  mode  number  k,  there  is,  across  the  planform  area  S,  a normal  pressure  force  per  unit  area  in  the 
direction  of  the  positive  z-axis,  called  the  aerodynamic  loading  at  the  point  (x,  y,  0)of  the  planform  S at  time  t. 
We  can  write  this  loading  in  the  form 

pV'='/*(x,y;^,Afcc)5*(w)e“'  (14) 

where  p is  the  density  of  the  fluid  in  the  uniform  flow  far  upstream  of  the  wing  and  Moo  is  the  Mach  number 

Moc  = - (15) 

a 

where  a is  the  speed  of  sound  in  the  uniform  flow  far  upstream  of  the  wing. 

On  using  the  governing  linearised  partial  differential  equation  for  the  perturbation  velocity  potential,  the 
boundary  condition  of  prescribed  upwash  on  S,  and  the  condition  of  no  loading  on  the  wake,  we  can  set  up  an 
integral  equation  relating  the  upwash  on  S to  the  loading  on  S (see  Ref.  5).  For  the  mode  number  k this  takes 
the  form 


«t(x,  )';*')  = 1 1 4(xo,  yo;  v,  ^—j^;  v,  exp  dxo  dyo  (16) 


where,  for  subsonic  flow, 

and 


-w/  du  , ,2Moo(AfooX  + R)  _(  iv(-x  + M^R\] 

‘ ■ /?(?ry^V^*n'TV  1-Ml  )\ 


R(x^+y^) 


R=y/{x^  + (l-Ml)y^}. 


(17) 


(18) 


For  dynamical  analyses  of  oscillating  wings  in  an  airstream  we  need  to  know  the  generalised  airforce 
coefficients  Qjkii',  M)  which  are  given  by  the  expressions 


0/k(v,  Mco)  = p 1 1 £Mo,  yo)lk(xo,  yo; M®)  dxo  dyo 


(19) 


in  the  linearised  approximation.  It  is  the  main  purpose  of  this  paper  to  discuss  the  numerical  evaluation  of 
OjkivjMco)- 


2.2.  Approximate  Solution  of  the  Integral  Equation 

The  integral  equation  (16)does  not  have  a unique  solution,  but  if  we  impose  the  condition  that  the  loading  at 
the  trailing  edge  of  S vanishes  then,  generally,  the  solution  becomes  unique.  As  a consequence  of  imposing  this 
condition  the  loading  acquires  a certain  behaviour  near  the  edges  of  the  planform  5 and  this  behaviour  is 
known.  We  shall  seek  an  approximate  solution  of  the  integral  equation  (16)  which  has  this  known  behaviour 
near  the  edges  of  the  planform  5. 

We  introduce  parametric  coordinates  on  the  planform  S by  means  of  the  formulae 


^o  = -7^{xo-XI.(yo)} 
c(yo) 


(20) 
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where  s is  the  semi-span  of  the  planform  S,c(y)  is  the  chord  length  and  XL{y)  is  the  x coordinate  of  the  leading 
edge  at  the  spanwise  position  y of  the  planform  S,  as  shown  in  Fig.  1. 

The  integral  equation  (16)  may  now  be  written  as 


We  now  take  an  approximation  fk(xo,  yo)  to  /*(xo,  yo;  v,  M^)  which  is  given  by  the  formula 

/'t(A:o.yo)  = ;;^exp(^^)  f X (22) 

c{yo)  ' I ' 1-1  p-i  ’ fo 

The  approximation  U(xq,  yo)  to  Ikixo,  yo;  v,  Af*,)  vanishes  at  the  trailing  edge  of  5 and  has  the  required 
behaviour  near  the  edges  of  S.  The  values  of  the  coefficients  Ak-i,p  are  such  that  the  simultaneous  linear 
equations  (33)  are  satisfied. 

The  approximation  tk{xo,  yo)  in  formula  (22)  has  been  expressed  in  terms  of  n chordline  base  functions  fo  '. 
/ = 1, 2, . . . , and  of  m spanwise  base  functions  T/g"*,  p = 1, 2, . . . , w.  The  formula  (22)  will  be  put  into  a 
different  form,  which  is  more  convenient  for  numerical  evaluation  than  is  the  formula  (22)  which  involves 
merely  simple  monomials  fo  ^o  ’■ 

The  points 


^i">  = |[l-cos(|^,r)], 


are  a set  of  n distinct  points  in  (0, 1).  We  form  interpolation  polynomials  /iJ"’(fo),  r=  1, 2, . . . , n,  based  on 
these  points.  These  interpolation  piolynomials  are  given  by  the  expressions 


M"'(fo)=_n(^7^),  r=l,2,...,n 


and  have  the  property 


where  is  Kronecker’s  delta 


a<"»(^!"')=5„ 


I f r = /■ 

'-■lo  <“) 

The  function  fJT',  i = 1, 2, . . . , can  be  expressed  as  a linear  combination  of  the  /ti"’(^o),  r = 1, 2, . . . , n, 
because  these  are  a set  of  n linearly  independent  functions  each  of  degree  (n  - 1)  in  ^o- 
The  points 


Tj^'"‘=cos(^j^y  p = l,2 m 


are  a set  of  m distinct  points  in  (-1, 1).  We  form  interpolation  polynomials  gi'"*(tjo),  s = 1, 2, . . . , m,  based  on 
these  points  and  given  by  the  expressions 


Es  (^o)  n i — (w«)_ —(«))’ 

p-i  Vp  ' 


5 = 1, 2, ....  m 


which  have  the  property 


= (29) 

where  S,p  is  Kronecker’s  delta.  The  function  rjg"*,  p = 1 , 2, . . . , m,  can  be  expressed  as  a linear  combination  of 
the  gi”\vo),  s = l,2,...,m,  because  these  are  a set  of  m linearly  independent  functions  each  of  degree 
{m  - 1)  in  Tjo. 

Therefore  /*(xo,  yo)  of  formula  (23)  can  be  expressed  in  the  different,  but  equivalent,  form 

4(xo.yo)=-r^exp(^)  i Z (30) 

C(yo)  ' < ’ Co 

where  the  unknown  coefficients  of  formula  (22)  have  now  been  replaced  by  the  new  unknown 

coefficients  Bit;,.,,  which  are  linear  combinations  of  the  Akxp-  We  need  determine  only  the  coefficients  B*;,., 
in  order  to  ..now  the  function  ft(xo,  yo)- 

The  n distinct  points  ^5"’,  / = 1, 2, . . . , n,  in  (0, 1)  and  the  m distinct  points  p = 1, 2, . . . , m,  in  (-1, 1) 
have  the  advantage  over  other  choices  in  that  simple  expressions,  concerning  integration,  which  have  been 
developed  in  Appendix  A,  may  be  used  in  the  ensuing  part  of  the  Report.  Other  choices  of  points  may  be  just 
as  good  as  far  as  the  numerical  accuracy  of  the  final  results  is  concerned,  but  it  is  possible  to  have  an 
unfortunate  choice  of  the  points  ^o  in  (0, 1 ) and  t/o  in  (- 1 , 1 ) resulting  in  ill-conditioning  of  sets  of  simultaneous 
equations  and  the  numerical  accuracy  of  the  final  results  is  poor  when  only  a moderately  small  number  of 
significant  figures  is  used  in  the  calculations. 

If  we  use  the  approximation  4 (xo,  yo)  from  (30)  for  k (xo,  yo;  v,  M^)  in  the  integral  equation  (2 1)  we  shall  get 
a corresponding  approximation  a/k(x,  y)  to  the  upwash  function  a^ix,  y ; v)  which  is  given  by 

C ” m / — ii/x\ 

ak(x,  y)  = yS^  t.^Bk-,r.sUrAx,y,f,M^)expy—j-j  (31) 

where 


UrAx,  y ; Moo)  = ^ g$'"'(7}o)Vr^  dr,o 

r=  1, 2, . . . , s = 1, 2, . . . , m.  (32) 

Following  Ref.  1 1 we  determine  the  unknown  coefficients  Bt;,.s,  r = 1, 2, . . . , « ; s = 1, 2, . . . , w,  from  the 
set  of  mn  linear  simultaneous  equations 


/»i'’\l-f)\/^{ak(x,y;»')-at(x,y)}exp(^^«f^  = 0, 

/■  = 1, 2, . . . , n;p  = 1, 2, . . . , m.  (33) 

The  corresponding  approximation  Ojk  to  the  generalised  airforce  coefficient  0,*  is  then  obtained  from 
formula  (19)  on  replacing  /*(xo,  yo;  v,  MA)  by  /It(xo,  yo)  and  is 


Oik  = p I j CAxo,  yo)/lk(xo,  yo)  dxo  dyo. 


(34) 


If  we  substitute  for  dk  (x,  y ) from  (31)  into  (33)  then  we  can  write  the  set  of  mn  linear  simultaneous  equations 
(33)  in  the  alternative  form 


7 i I >l'i.p.r..Bk.r.,  = dkxp,  / = 1, 2, . . . , n; p = 1, 2, . . . , m, 

/ /—I 


(35) 
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where 


Ok.i.p  - j ^ gp'"‘(r?)Vl-i7^  dtf  a*(;c,  y;  if)  exp  df 

/ = 1,2,  ...,rt;p  = l,2,  ...,m,  (36) 

and 

= I ^ gn^)Vl-7,^  dT,  jj  y ; »-.  Moo)  df, 

/=  1, 2, . . . , rt;p  = 1, 2, . . . , /«,  /■  = 1, 2, . . . , «;  j = 1, 2, . . . , /n.  (37) 

We  assume  that  the  wing  planform  S is  symmetric  about  the  x-z  coordinate  plane  and  that  the  modes  of 
oscillation  are  either  all  symmetric  modes  or  all  antisymmetric.  Then  we  must  have 


c(-y)=c(y) 

(38) 

XLi-y)=XL{y) 

(39) 

ik(x,-y)  = K^k{x,y) 

(40) 

where  k — +1  for  the  case  in  which  all  the  modes  are  symmetric  and  k = — 1 for  the  case  in  which  all  the  modes 
are  antisymmetric. 

We  note,  from  formula  (27),  that 


_,(m)  (m)  _ 1 ^ 

Vm-p-hi  ~~  y}p  t p — ■ 1, 2, . . , , nt. 

(41) 

Then,  it  follows  from  (28)  that 

g^m-.+i(v)  = gi^X-v),  s = 1, 2, . . . , m. 

(42) 

from  (32)  that 

Ur.m-,*i(x,  y ; If,  Moo)  = C/,.,(x,  -y ; v,  Af®),  r = l,2,  ...,«;s  = l,2, .. 

• , m. 

(43) 

because  the  kernel  function  K(x,  y;  if,  M)  is  symmetric  in  y,  and  from  (37)  that 

»+l  “ ^i,p\r,at  / — 1,  2,  . . . , /T  J P = 1,  2,  , , . , ftl,  T = 1,  2,  . . . , W J 5 = 1,  2,  . . , , ffl.  (44) 
From  (40)  and  (1 1)  we  get  that  the  reduced  upwash  function  a^ix,  y ; if)  satisfies 

aic(x,-y  ;p)  = /ca* (x,  y;  if),  (45 ) 

and  then,  from  (36),  we  get  that 

^k;i.m-p+l  ~ l(^k-,i,pt  / ~ 1,  2,  . . . , n J p = 1,  2,  . . . , fit,  (46) 

Because  of  (44)  and  (46),  the  coefficients  r = 1,  2, . . . , n ; s = 1, 2, . . . , m,  which  satisfy  the  set  of  mn 

linear  simultaneous  equations  (35),  must  satisfy  the  relations 

^k;r.m— »+l  KSk;r,tt  r = 1,  2,  , . , , rt  J S = 1,  2,  , , , , m.  (47) 
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We  now  define  the  integer  ntn  by  means  of  the  formula 

mH  = 4{2m+K~K{-l)'”}.  (48) 

We  note,  from  (46)  and  (47),  that  when  m is  an  odd  integer  and  »(  = -1  we  must  have 

0k...m„*l=0,  (49) 

and 

(50) 

For  all  values  of  m and  k,  only  wih  values  of  the  .r.s  are  therefore  unknown  and  consequently  the  set  of  mn 

linear  simultaneous  equations  (35)  may  be  replaced  by  the  set  of  linear  simultaneous  equations 


where 


and 


Iri 


mu 

1 '^i.p 

,r,s^k',r. 

^kii.py 

1 = 1,2,.. 

.,n\p- 

1,2,. 

■ . , niH, 

(51) 

s®*  1 

Ck-.r,s  — 

s,  r=l. 

2, . . . , n;s 

= 1,2,.. 

■ , 

-1, 

(52) 

1 

\^2Bkxm„ 

m even 

Ck\r, 

s2Bk-,r,m„ 

m odd,  K 

= -l. 

(53) 

[Bk-,r,mH 

m odd,  K ■ 

= +l, 

r = l. 

2, . . . , «, 

tl-J  + l). 

i = 

1, 2, ...,« 

;p  = l,2,.. 

= 1,2, 

. . . , n ; 5 = 1, 2, . . 

. .,Wh.  (54) 

Furthermore,  the  set  of  linear  simultaneous  equations  (51)  may  be  written  as  the  matrix  equation 


ym[a]=[»fc]  (55) 

where  ['F]  is  a square  matrix  of  order  rtiHH  x mnn  consisting  of  elements  which  are  the  quantities  'F, 
j = 1,2, . . . , n;  p = 1, 2, . . . , mn,  r = 1, 2, . . . , n;  s = 1, 2, . . . , wh,  [G]  is  a column  matrix  consisting  of  the 
rtinn  elements  Q;,.,,  r=  1, 2, . . . , r = 1, 2, . . . , mn,  and  [^t]  is  a column  matrix  consisting  of  the  m/pi 
elements  / = 1, 2, . . . , n ; p = 1, 2, . . . , mu-  The  arrangement  of  the  elements  Ck-,.,  in  the  column  matrix 
[Cic  ] is  immaterial  and  so  is  the  arrangement  of  the  elements  the  column  matrix  [0*],  but  once  these  two 

arrangements  have  been  specified  the  arrangement  of  the  elements  in  the  square  matrix  ['I^]  is 

determined. 

To  be  definite  we  may  specify  the  following  arrangement  for  the  elements  in  the  matrices  [C*],  [0*]  and  ['?]. 
The  column  matrix  [C*]  has  the  element 

Ck.,.«  T=\,2,...,n-,s  = \,2,...,mH,  (56) 

in  the  n(mH  — s)+/th  row.  The  column  matrix  [^*]  has  the  element 

Bk\i,pj  t = 1,  2, . . . , n , p 1,2,...,  (57) 

in  the  n(mH -p)+n-i  + 1st  row.  Consequently  the  square  matrix  ['I']  has  the  element 

I = 1,  2, . . . , n;  p = 1, 2, . . . , wh,  r = 1, 2, . . . , n ; s = 1, 2, . . . , wh,  (58) 

in  the  n(»jH-p)  + « ~/  + 1st  row,  and  «(»»«  — 5)  + rth  column. 
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2.3.  Evaluation  of  the  Generalised  Forces 


or’- 1’' 


If  we  use  the  numerical  integration  formulae  (A-16)  and  (A-35)  with  n integration  points  chordwise  and  m 
integration  points  spanwise  to  evaluate  Xr.r.’  and  from  formulae  (60)  and  (61)  respectively  we  get 

Xi.r..  - I Z V;(xl'J'"\  y'-^)  exp 

1-1  P-1  ' / ' 

= V/(^‘rr y exp  r=  1,2 n ; r = 1 , 2, . . . . m,  (73) 


= f/^”^GlrW(xjrr\ y^r'-.  exp  (yxlr ').  « = l,  2, . . . , « ; P = l,  2 m. 

On  substituting  from  (73)  and  (74)  respectively  into  (66)  and  (67)  respectively  we  then  get 


Ci;r.s  ^ Ci(xir\  yi""),  r=l,2,...,n;s=l,2,...,m. 


ak;i.p-a*(x/.p'"\yp'"’; »').  / = l,2,...,«;p  = l,2,...,w.  (76) 

The  formulae  (73)  and  (74)  are  only  approximate  formulae,  but  they  may  be  very  good  approximations  to 
the  correct  values  at  high  values  of  n and  m if  c(y ),  a*  (x,  y ; p)  and  i,(x,  y ) are  continuous  functions  with  only  a 
few  undulations  over  the  planform  5.  There  are  instances,  e.g.  control  surface  rotation,  in  which  their  accuracy 
is  not  as  good  as  may  be  desired  unless  n and  m are  unpractically  high.  If  it  is  considered  that  the 
approximations  (73)  and  (74)  are  sufficiently  accurate  then  formulae  (75)  and  (76)  may  be  used  for  {/■,,,  and 
«k;i.p  respectively  instead  of  the  expressions  (66)  and  (67)  with  ^ic;<.p  being  obtained  accurately  by 

'imerical  integration  of  the  integral  relations  (60)  and  (36)  using  formulae  which  are  much  more  accurate  than 
tormulae  (73)  and  (74),  or,  indeed,  by  carrying  out  the  integrations  analytically  if  this  is  possible. 

Let  Hi]  be  the  row  matrix  of  mnn  elements  with  the  element 

(i-.r.,,  r=l,2,  ...,n;s  = l,2,  ...,mH,  (77) 

in  the  n(mH~s)+rth  column.  Let  [o*]  be  the  column  matrix  of  elements  with  the  element 

ak  xp,  1 = 1,2 n ; p = 1,  2, . . . , mw,  (78) 

in  the  n(wH-p)+«  -/  + 1st  row.  Let  [E)  be  the  diagonal  matrix  of  order  wh«  x mutt  with  the  element 

^"’Gi'"’  exp  (^xirj , r = 1 , 2 n;s  = 1,2 m„. 
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(79) 


in  the  n(mH  - s)+ rth  row  and  column.  Let  [D|  be  the  diagonal  matrix  of  order  mz/n  x mnn  with  the  element 


¥ 


exp  (jxirj.  / = 1. 2, . . . . n ; P = 1,  2, . . . , m„, 


(80) 


in  the  n(mH-p)  + n-i  + 1st  row  and  column.  Then  we  can  write  formula  (66)  and  (67)  as  the  equivalent 
matrix  formulae 


and 


[fl*]  = [£)][a*] 

respectively.  On  substituting  from  (81)  and  (82)  into  (65)  we  get 


(81) 

(82) 

(83) 


which  is  the  final  formulation  for  the  approximation  (5,*  to  the  generalised  airforce  coefficient  O/i. 

Suppose  that  there  are  P modes  of  oscillation  corresponding  to  the  index  k in  the  summation  on  the  right 
hand  side  of  (6)  taking  the  values  k = \,2, . . . ,P,  only,  just  as  in  formula  (10).  Let  [f]  be  the  matrix  of  order 
P X mnn  obtained  by  arranging  the  row  matrices  [^;],  j=\,2,...,P,  sequentially  one  below  another.  Let  [a] 
be  the  matrix  of  order  rrinn  x P obtained  by  arranging  the  column  matrices  [a*],  Jt  = 1, 2, . . . , P,  sequentially 
one  alongside  another.  Let  [(5)  be  the  square  matrix  of  order  P x P which  has  the  element  [<^/*]  in  the  /th  row 
and  kth  column,  / = 1, 2, . . . , P;  fc  = 1, 2, . . . , P.  Then  we  get  immediately  from  formula  (83) 

(84) 


2.4.  Evaluation  of  the  Loading 

We  can  obtain  the  approximations  tk{xo,  yo),  k = 1, 2, . . . , P,  to  the  loading  distribution  hixo,  yo;  M)  in 
the  mode  k of  oscillation  directly  from  formula  (30).  By  using  the  definitions  (52)  and  (53)  in  formula  (30)  in 
order  to  replace  the  coefficients  by  Q;,.,  we  get 

/'*(^o,yo)  = :^exp(=^)  I ? (85) 

2C(yo)  ^ ’ Co 


where 

ki'””'(T7o)  = g^J”\r)o)+ Kg'r\—no),  s = 1, 2 niH.  (86) 

The  relations  (42)  and  (47)  have  also  been  used  in  obtaining  the  formula  (85). 

Let  [F(xo,  yo)]  be  the  row  matrix  of  mnn  elements  with  the  element 

sciyo)  ^ ~ r=l,2, ...,«;  s = 1,2,  ...,mH,  (87) 

in  the  n(mH-s)+nh  column  and  let  [/)k(xo,  yo)]  be  the  matrix  consisting  of  the  single  element  fk(.Xo,  yo)-  Then 
we  can  replace  the  formula  (85)  by  the  matrix  formula 

[/t(JCo,  yo)]  = 2 yI^(-*o,  yo)][C*].  (88) 

If  we  obtain  [C*]  from  formula  (55)  and  express  [tf*]  in  terms  of  [a*  ] by  means  of  formula  (82),  then  we  can 
write,  instead  of  (88) 

l/lk(A:o,  yo)]  = i[F(xo,  yo  )]['?]■'[£>][«*]  (89) 

which  is  a formula  for  determining  fkixo,  yo)- 
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In  particular,  if  we  put  in  formula  (89) 


•*0  ^ ‘*i,p 

yo  = yp"’ 


where  y},"’  are  defined  in  formulae  (68)  and  (69),  with  i and  p replaced  by  r and  s,  we  get 

[/;(jr!.r‘.  y|r^] = j[F(xir\  yS,"’')im-‘[o][a*i.  (91) 

We  note  that  the  row  matrix  [F(x!.p'"’,  yp’"^)]  has  only  one  non-zero  element,  namely  the  one  in  the 
n (niH  -p)+ith  column.  This  fact  comes  from  applying  the  properties  (25)  and  (29)  of  the  functions  h and 

g,iVo)  to  the  functions  fci'""’(T/o)  of  (86)  and  the  elements  (87)  of  [F(xo,  yo)]- 

Now  let  [/*]  be  the  column  matrix  of  rtinn  elements  with  the  element 

/'t(4r’.yp’"‘),  / = l,2,...,«;p  = l,2 (92) 

in  the  n (ntn  -p)+ith  row  and  let  [F)  be  the  matrix  of  order  ttiun  x niipi  whose  n (mn  - p ) + fth  row  is  the  row 
matrix  F(xj^,  yjr^).  Actually,  according  to  the  above,  the  square  matrix  [F]  is  a diagonal  matrix.  It  can  be 
written  as  the  product  of  diagonal  matrices 

[F]  = [//][E]  (93) 

where  the  diagonal  matrix  [E]  has  been  defined  immediately  before  expression  (79)  and  [H]  is  the  diagonal 
matrix  of  order  rtinn  x ninn  with  the  element 

1^)  / = 1, 2, . . . , n ; p = 1, 2, . . . , mH,  (94) 

in  the  «(/«//— p)  + fth  row  and  column. 

From  equation  (91)  for  / = 1, 2, . . . , n ; p = 1, 2, . . . , Wh,  we  then  deduce  the  matrix  equation 

[/;]  = |[//][E]m-*[D][a,].  (95) 

If  we  solve  equation  (95)  for  [^]~‘[I>l[aik]  and  substitute  the  result  into  formula  (89)  we  get 

[Lixo,  yo)]  = [F(xo,  yo)l[F]-‘[//]-'[r»]  (96) 

which  is  an  alternative  formula  to  formula  (89)  for  determining  [/“(.(xo,  yo)]-  Formula  (96)  is  useful  if  one 
computer  program  is  constructed  to  evaluate  the  rrinn  values  /*(Xj^p'"\  yp”'),  / = 1, 2, . . . , n ; p = 1, 2, . . . , ntn, 
for  each  value  of  k from  equation  (95)  and  a second  computer  program  is  constructed  to  evaluate  /'*(xo,  yo)  at 
given  values  of  Xo  and  yo  using  these  mun  values  /lt(x,'.p'”’,  yp'”’), » = 1,  2, . . . , n;  p = 1,  2, . . . , mn,  for  each 
value  of  k in  formula  (96). 

We  notice  that  the  matrix  product  [E]['I']“'[D]  occurs  both  in  formula  (84)  for  the  generalised  airforce 
coefficients  and  in  formula  (95)  for  the  loading  at  the  loading  points.  It  was  to  achieve  this  that  the  product 
formula  (93)  for  [F]  was  introduced. 

If  [/^  is  the  matrix  of  order  ntHfixP  obtained  by  arranging  the  column  matrices  [4],  Ic  = 1,  2, . . . , F, 
sequentially  one  alongside  another,  and  if  [/“(xo,  yo)]  is  the  row  matrix  of  P elements  obtained  by  arranging  the 

elements  /fc(xo,  yo),  fc  = l,2 P,  sequentially  one  alongside  another,  then  we  get  immediately  from 

formulae  (95)  and  (96)  the  formulae 

[n  = i[//][F]m-'[D][«]  (971 


[^(^0,  yo)]  = [F(xo,yo)][F]-'[//]-'[f]. 


J 


3.  Namerical  faitesnitioii 

The  quantities  4>i.p.r.s  defined  in  formulae  (37)  as  double  integrals,  are  to  be  evaluated  numerically,  and  we 
do  this  by  evaluating  the  double  integrals  using  the  Gaussian  numerical  integration  formula  (A- 16)  over  N 
chordwise  points  f in  (0, 1)  and  the  Gaussian  numerical  integration  formula  (A-35)  over  M spanwise  points  rf 
in  (-1, 1).  This  process  leads  to  the  formula 


i 1, 2, . . . , /t , 1,2,...,  /n,  /'“1»2,  — 1,2, .,,,  w,  (^9) 


where  (see  Appendix  A) 


6 -2[i- 


cos 


(100) 


= 2^0-^n,  7=1, 2,...,  a;  (102) 


Gy*^=|  ^ gjr*^(r;)v'l-Tj^  dr) 


"m+I^^  ('jV^VI,  7-1.2 M. 

(103) 

= c(yr^l  -^n  + XL(yS*'^  7=1,2,.. 

and 

..N;/=l,2,.. 

.,Af, 

(104) 

yr=sr,r. 

(105) 

We  define  the  integer  Mh  by  means  of  the  formula 

MH  = k{2M+K-{-l)'^K}. 

(106) 

Then,  if  we  use  the  relations  (41),  (42)  and  (43)  in  (99),  and  use  the  definitions  (54)  we  get 


where 


/»1 /»1 


' = 1,2 n ; p = 1,  2, . . . , /Hh,  r = 1,  2, . . . , n ; s = 1,  2 ntn,  (107) 


Lpu 


L 


pMh  - i 


= g'p”Wn+Kgl”\-r,f^ 
g^W^-g^X-v^i^) 


«r’(0) 


J 

A/ even 
M odd,  x = — 1 
M odd,  »f  = +1. 


(108) 

(109) 


1 


j 


I 

i 

i 

i 

i 

j 

) 

1 


j 

j 


I 


i 
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We  may  now  write  for  the  matrix  appearing  in  formula  (55), 


m=[L][t/] 


(110) 


where  [L]  is  the  matrix  of  order  mnny-MnN^  with  the  element 


. / = 1, 2, . . . , p = 1,  2, . . . , Wiv,  /=  1. 2 N-,J=\,2,...,Mh 


(111) 


in  the  «(mH—p)+n  — i + 1st  row  and  N{MH—J)+N~I+lst  column,  and  [U]  is  the  matrix  of  order 
MnNy.  ntHn  with  the  element 


UrAxTj'^\  yf‘\  yT\  v,  m*). 


r=l,2,  s = 1,2, 


(112) 


in  the  N(Mh-J)+N-I+  1st  row  and  rt(mH-s)+rth  column. 

We  note  that  the  matrix  ['P]  given  by  (1 10)  is  non-singular  only  if  MnN^mun.  We  shall  take  Mh  ^ mn  and 
N^n  so  that  this  condition  is  satisfied  and  then  ['I']  can  be  inverted  to  give  [^]”‘  required  in  formula  (65). 
Now,  from  (32),  we  get 

,r  , ^ f*' gi'"'(»7o)'^l->7o  Knv,  N J 

Ur,,(x,  y;  V,  Moo)  = -2  ^2 "r  (^,  TJ,  7]o\  P,  Ma,)  dr}o  (113) 

s J-i  \'n~Vo) 


where 


V,  Vo;  K f M„)  dio-  (i  14) 

To  evaluate  Ur.s(x,  y : v,  Mao)  from  (113)  we  write,  following  Ref.  12, 


,, , ...  r^g^r\vo>n^o.. 

l/...(x,y;.,M„)  = pJ__ 

x|/t"’(f  7j,  T/o;  p,  Mao)-ii"K^,  r?,  t?; M<=o)-(vo-v)^-^^&\4,  V,  Vo;  v,  A/oo)j  I d-no  + 

M.)  f dr,,- 

s i-x  (v-Vo) 

/'r  a f*' gi”\vo)^T^o  . 

— ^ — "r  (f  V,  Vo;  P,  Mao)  ; r drjo.  (1 1 5) 

s Lapo  (V-Vo) 


Now 


/ ~^—^li''^(irVrVo;P,Mao)-/i'’\^,V,V;i',ll^'=o)-(VO-V)f~^r"^(^,V,Vo;P,Mao)]  j (116) 

(V-Vo)  LdT/o 


becomes  logarithmically  infinite  when  rjo  - r?,  so  that  a numerical  integration  process  for  evaluating  directly 
the  first  integral  on  the  right  hand  side  of  formula  (1 1 5)  must  not  have  an  integration  point  at  rjo  = rj-  Also  it  is 
not  easy  to  estimate  the  accuracy  of  such  a numerical  integration.  We  can  replace  the  function 
Vr  Vo,  p.  Moo)  in  the  first  integral  on  the  right  hand  side  of  formula  (115)  by  a function 
/•,"’(f,  rj,  Tjo;  p,  Mac)  so  that  numerical  integration  of  the  resulting  integral  is  straightforward  using  a Gaussian 
numerical  integration  process  and  analytical  integration  of  the  difference  between  the  two  integrals  can  be 
carried  out.  This  is  achieved  by  using  the  known  analytical  behaviour  of  v,  Vo;  p,  Mao)  for  r/o  near  to  rj 
(see  Ref.  5,  Appendix  4). 
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We  introduce  the  functions  //"’(f,  tj,  tjo;  i'.  A/®)  by  means  of  the  formulae 

17,  »}o;  J/,  Mco)=t]r''\^,  tj,  Tfo; Moo)+F‘,"’(^,  rj; login -tjoI,  r=  1,  2, . . . , n,  (U7) 


where 


r=1.2,....n.  (118) 

The  lowest  order  logarithmic  contribution  to  r\,  170:  v,  M®)  is  then  completely  separated  as  the  second 
term  on  the  right  hand  side  of  (117).  It  follows  that  the  function 

’I’  Afoo)-  n.  n ; a/®)-  (no-  ’i-  ^o;  j',  A/®)j  |, 

/•=1.2 « (119) 

is  finite  when  no  = n if  J^rfy)  and  c“{y)  are  finite.  The  function  (1 19)  is  then  finite  for  all  170  in  (-1, 1)  and 
therefore  is  straightforward  to  deal  with  as  part  of  an  integrand  of  an  integral  which  is  to  be  evaluated  using  a 
Gaussian  numerical  inte^ation  process. 

If  we  substitute  for  t/,  no; A/®)  from  (117)  into  the  first  integral  on  the  right  hand  side  of  (1 1 5)  we 
obtain 

y;  Af®)=p J ^ - |/‘r"’(f  n. no:  Af®)-/^,"*(f  t?. n;  y. A/®)- 

~(no-n)[^/i"’(^.  n.  no; A/®)]  | di7o+ 

L3no  Jtio-i* 

/2 

+ -if^r"'(f  n ; A/»)  g‘,'"'(no)  log  In  - nol'/i-no  </no  + 

j j_i 

n.  n; Af®)  [ drjo- 

s J-i  (n-no) 

l^r  d If**  e^'”Xv  ) 

— j — n.  no;  i',  A/®)  [ l - r^l  drio.  (120) 

s Ldno  (n-no) 

It  is  observed  that  the  third  and  fourth  integrals  on  the  right  hand  side  of  (120)  are  principal  value  singular 
integrals. 

The  first  integral  on  the  right  hand  side  of  (120)  may  be  evaluated  numerically  by  means  of  a Gaussian 
numerical  integration.  The  obvious  formula  is  the  one  with  Vl  -tjo  as  weight  function,  but  we  shall  find  it 
necessary  to  use  also  other  weight  functions  in  order  to  justify  the  final  form  (156) for  Ur.,(x,  »<,  A/®)  in  all 

circumstances.  These  weight  functions  are  V(1  - t7o)/(1  + T70),  '/(I  + ’7o)/(l  ->70)  and  l/Vl  -tjq.  If  we  use  A 
integration  points,  the  application  of  the  Gaussian  numerical  integration  corresponding  to  each  of  these 
weight  functions  gives  the  result 

“K,?.  (W?? 

Lot^O  JT)0«y|J 
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The  values  of  Pa  and  p = 1, 2, . . . , A,  for  the  different  weight  functions  may  be  obtained  from  equations 
(A-86),  (A-88),  (A-91)  and  (A-94)  of  Appendix  A. 

For  the  weight  function  Vl-rjo  these  values  are 


Pa  = (A+1) 


(122) 


and 


/(A)_  ^(A) 
sp  Vp 


For  the  weight  function  >/(l  -77o)/(l  + 170)  these  values  are 

Pa  = (A  + §) 
and 


(123) 


(124) 


,(A)_  ^(2A) 

4p  Vp 


-““(sfe)-  


For  the  weight  function  V(1  + t/o)/(1  — ^to)  these  values  are 

Pa  = (A  + J) 

and 


(125) 


(126) 


/(A)_  „(2A) 

ip  V2p-l 


= P = 1-2,...,A. 


For  the  weight  function  l/Vl  -770  these  values  are 


Pa  = A 


(127) 


(128) 


and 


)-(A)_„(2A-1) 
ip  ~V2p-l 


= cos 


p = 1.2,....A. 


(129) 


Formula  (121)  is,  in  general,  only  an  approximate  formula  but  for  the  weight  function  •Jl-vo  it  would  be 
exact  if  the  function  (119)  were  a polynomial  of  degree  «2A-m  in  tjo-  For  the  weight  functions 
'/(l“Tjo)/(l  + i7o)  and  V(1  + t7o)/(1-i7o)  it  would  be  exact  if  the  function  (119)  were  a polynomial  of  degree 
«2A-m-l  in  770  and  for  the  weight  function  l/'Jl-rtl  it  would  be  exact  if  the  function  (119)  were  a 
polynomial  of  degree  <2A  — m — 2 in  770.  However,  since  the  integration  weights  are  positive  numbers  and  the 
function  (119)  is  finite  for  all  170  in  (-1, 1),  the  right  hand  side  of  (121)  will  converge  to  the  exact  value  of  the 
integral  on  the  left  hand  side  of  (12 1 ) as  A tends  to  infinity.  We  can  take  A to  be  large  enough  for  formula  (121) 
to  be  sufficiently  accurate  for  all  practical  purposes. 

We  note  that  the  form  of  the  right  hand  side  of  (121)  can  be  used  only  if 


p = l,2,  ...,A, 


(130) 
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but  if  this  is  not  the  case  the  right  hand  side  of  (121)  is  easily  modified  and  the  resulting  formula  involves  the 
second  derivative  rj,  t/o;  M^)  for  tjo  = for  the  appropriate  value  of  p. 

We  can  further  write 

J ^ g^rXvo)  log  h - Tjol'/ 1 - Tjo  dr)o 

= 1^  log  |r7-’?ol'/l  -rjo  dVo  + g'.^Xv)  ^ log  |tj -tjo|>/1  - rjo  dr/o 


= -f  i log h -^<^1(1  -(^5.^y}  + ^gl'"'(r,){Tj^-|-log2} 

“ap-1  2 


(131) 


where  one  of  the  above  Gaussian  integration  formulae  with  A integration  points  has  been  used  to  give 


I ^ {gi”’(»Jo)-  gi'" ’(»»)}  log  |i7  - Tjo|Vl-Tjo  drjo 


= f I {gi’"*(^i.^‘)-gi""(r,)}logh-^‘''1{l-(f‘"y}. 

r^Ap-i 


(132) 


The  formula  (132)  is  only  approximate  but  since  the  integration  weights  are  positive  numbers  and  the 
function  {gi’"’(T/o)-gi'"’(r/)}  log  |tj  - 7jo|  is  finite  for  all  tjq  in  (-1,  l)the  right  hand  side  of  (132)  will  converge  to 
the  exact  value  of  the  integral  on  the  left  hand  side  of  (132)  as  A tends  to  infinity.  We  can  take  A to  be  large 
enough  for  the  formula  (1 32)  to  be  sufficiently  accurate  for  all  practical  purposes.  Actually,  the  integral  on  the 
left  hand  side  of  (131)  can  be  evaluated  analytically,  but  it  is  simpler  to  get  its  value  from  the  expression  on  the 
right  hand  side  of  (131),  and  then  the  accuracy  is  consistent  with  the  accuracy  of  formula  (121). 

On  substituting  from  (121)  and  (131)  into  (120)  we  get 

U (x  V V M \ = ——  y ii£— LLx 

Ur.s{.X,  y , P,  nJx,)  p 2 2-  , (aK2  ^ 

r\  s p-i  (7j-4p  ) 

X [/<">(f  n,  Af«)+/=1'’>(f,  r?; M^)(v  log  k + 

‘‘"-‘-rj, — on?? — J' 

S L3i7o  Jno-»t 

. r f g\'"\Vo>f^O  . TT  a g<'">(^<A)){l  - (^<Ay}l 


V, p,Moo)gr(v)[^(v^-k-\og2)~  z iogi»,-^n{i-(^<^y}] 


= -n-5  L 7 7(ak2 (f  rt.^p  : r'.M®)- 

r\  S p-i  (V~(p  ) 


(f  r,,r,,.,M«)[J  ^ (,-2,0)^  J" 

^ M M - , U ll  rf"'«*'"’(^o)Vr^  . rr  A g<r>(f','>Kl-(0^}l 

Ji,  — on?) — J 

+ i logll-rt'lll-K?)’)!.  (133) 

S L Z r\ pm\  J 


19 


Again, 


i: 


I - vl 


(V-Vo) 

»+i 


dvo 


f*'{gr(vo)-grM-(vo-v)gr^(v)}  r. — ? . , 

= } rj Vl-i7oai7o  + 

J-i  (v~Vo) 

‘”’1.  ‘”'1. 

■p.pt.  ^ 

- ^g^s^Xv)- 


(134) 


where  the  principal  value  singular  integrals  jtj  '/(I -i7o)/(rj  - nof  drjo  and  jlj  N/(l-7jo)/(r7  - Vo)  dvo  have 
been  replaced  by  their  respective  values  -tt  and  irrj.  Similarly 


i: 


gi'"’(’7o)V'l-T?o 


(v 

-Vo) 

=r 

j-i 

TT 

A {g! 

"n, 

L — 
► = 1 

dvo 


iv-Cn 


Ta5^ {1  - (^'p'^Vl  + irvgi^Xv  )• 


(135) 


In  equations  (134)  and  (135)  one  of  the  above  Gaussian  integration  formulae  with  A integration  points  has 
been  used  to  give 

r'{g'r\vo)-gi"'\v)-(vo-v)g'r\v)}  r. — 5 . 

J-, 

p ^ / A\y\2  U Up  ) I 

i A p-1 


iv-CXi^ 


(136) 


and 


f"MgrW)-gr(r7)} /, 


^ ^ {grW-g^r\v)}^^ 


P -1 


^ bP 


The  function  g^r^'iv)  is  the  derivative  with  respect  to  77  of  gj^^rj).  The  formula  (136)  is  exact  if 

A ^2^  — 1 


for  weight  function  Vl  -770,  if 

A^lm  -2 

for  weight  functions  '/(I -77o)/(l +770)  and  V(1 +77o)/(l -770),  and  if 

A>|m 

for  weight  function  1/Vl  -770. 


(137) 


(138) 


(139) 


(140) 
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The  formula  (137)  is  exact  if 


for  weight  function  Vl  -tjo.  if 


A^irn  -5 


(141) 


(142) 


for  weight  functions  '/(I  -7jo)/(l  +t?o)  and  '/(1  + t/o)/(1  — tjo).  and  if 


(143) 


for  weight  function  1/Vl— t)o. 

Hence,  if  one  of  the  conditions  (141),  (142),  (143)  holds,  as  appropriate  to  the  weight  function  used,  then 
both  formulae  (134)  and  (135)  are  exact  and  we  can  substitute  them  into  the  right  hand  side  of  (133)  to  get  an 
alternative  formula  which  is  precisely  equivalent  to  (133), 


IT  I ^ 2 

UrAx,y,p.M^)=^—-^  I ^ 
r \ S p-i 


(v-cy 


17, 17 ; r,  Mco)gi’"\v)\  1 Z 

5 L r\ p«i 


(v-cyi 


-rr—  —{Ir  (1, 17.  ■'Jo;  v,  Moo)gs  (170)}  ^ MJ  I TwT  + 

S LaT7o  J,Io-.l  L ^^Ap=l  (17 -Sp  1 j 


+ 7rW’(f  r7;.',A/oo)g‘r‘(ri)x 


x[|(T,^-i-log2)-^  i log (0^)1  (144) 

L * A P“1  -* 

and  this  is  the  numerical  integration  formula  that  we  use  for  evaluating  (x,  y ;v,M).  However,  we  need  the 
values  Urj(x,y  ; v,  M)  only  for 


y = yT\  /=1,2 M, 

in  order  to  evaluate  the  elements  (1 12)  of  the  matrix  [t/],  and  these  values  of  y correspond  to 

n=nT\  /=1,2,...,M. 


(145) 


(146) 


With  the  values  of  1)  given  by  (146)  a simplification  of  the  expression  (144)  occurs  with  certain  special  values  of 
A,  thesespecial  values  of  A being  specific  to  the  weight  function  Vl  -T70,  V(1-t7o)/(1+t7o).  "^(l  ■^•i7o)/(1  ~i?o) 
or  l/'/l  -T70  being  used.  In  fact  we  have  (see  Appendix  B) 


{1-(01 

{■nr-O 


= 0 


(147) 


_ j_  A {i-im 

Pap-.  {-nT'-Cy 


= Pa 


(148) 


in  the  following  five  cases; 

Case  (i)  Any  positive  integer  M,  any  positive  integer  J,  1 and  weight  function  1/Vl  -T70. 


Values  of  A are  given  by 

A = a(M+l)  (149) 

where  a is  any  positive  integer. 

Case(ii)  Any  odd  positive  integer  Af,  any  even  positive  integer  /,  2=sy«M-l,  and  weight  function 

i/VT^. 

Values  of  A are  given  by 

A = i(2a-1)(M+1)  (150) 

where  a is  any  positive  integer. 

Case(iii)  Any  even  positive  integer  M,  any  even  positive  integer  J,  2^/ssAf,  and  weight  function 
Values  of  A are  given  by 


A = 5(2a-l)(A/+l)-5 


(151) 


where  a is  any  positive  integer. 

Case(iv)  Any  even  positive  integer  M,  any  odd  positive  integer  J,  l=s/«Af-l,  and  weight  function 
y/l-Vo/l+Vo- 
Values  of  A are  given  by 


(152) 


where  a is  any  positive  integer. 

Case  (v)  Any  odd  positive  integer  M,  any  odd  positive  integer  7,  1 M,  and  weight  function  -Jl-rtl. 
Values  of  A are  given  by 


A = i(2a-1)(M+1)-1 


(153) 


where  a is  any  positive  integer. 

If  we  substitute  from  (147)  and  (148)  into  (144)  for  the  values  of  y given  by  (145)  we  get 


“ f>  7^  ^ / (mV  ►(AK2  (v*  » ip  , Moo) 

M\  S p— I (17  J »P  ) 

T,r’.  Tjr*;  MAgr\r,'r)  + ir‘-FT%  r,r-,  P,  Mo.)g'-r\riD^ 

s s 


xfi((r,ry-i-iog2)-^  i iogii7r’-^nx(i-(^rv}l-  (154) 

L r\p^\  -* 

It  is  to  be  noted  that  the  quantity  [(d/dTjo){/J"’(f,  t\T\  Vo',  v,  Moc)gi'"*(7}o)}],o-,<*‘>  does  not  occur  in  formula 
(154). 

Formula  (154)  is  valid  for  the  following  cases: 

Case  (i)  Any  positive  integer  M and  any  positive  integer  7, 1 « 7 =s  M,  when  A is  given  by  ( 1 49),  is  given  by 

(128)  and  p = 1, 2. ....  A,  are  given  by  (129). 

Case  (ii)  Any  odd  positive  integer  M and  any  even  positive  integer  7, 2«7«M-1,  when  A is  given  by  ( 1 50), 

P\  is  given  by  ( 1 28)  and  p = 1 , 2 A,  are  given  by  (129). 

Case  (iii)  Any  even  positive  integer  M and  any  even  positive  integer  7, 2 « 7 « M,  when  A is  given  by  ( 1 5 1 ),  Pa 
is  given  by  (126)  and  fp^\  p = 1, 2, . . . , A,  are  given  by  (127), 

Case  (iv)  Any  even  positive  integer  A/ and  any  odd  positive  integer  7, 1 «7«Af- 1,  when  A is  given  by  (152), 
Pa  is  given  by  (124)  and  p = 1, 2, ....  A,  are  given  by  (125). 

Case  (v)  Any  odd  positive  integer  A/ and  any  odd  positive  integer  7, 1 «7«A/,  when  A is  given  by  (153),  Pa  is 
given  by  (122)  and  i'p^\  p = 1, 2, . . . , A,  are  given  by  (123). 
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By  using  formula  (A-55)  we  can  put  the  result  (154)  into  yet  a different  form,  in  which  we  do  not  have  to 
distinguish  between  the  cases  (i),  (ii),  (iii),  (iv)  and  (v).  If  we  put 

M = q(M+l)-\  (155) 

then  the  result  (154)  can  be  expressed  in  the  form 

Ur  Ax,  yT\  f,  MA) 


-iL  y /<’•)(£  „(*).  „ jvf  f 

^ K^iVJ  iVp  y^i^<x>)Ss  \Vp  ) I / (MK2 

s p-i  J-1  (Vo-Vj  ) 

+ ‘-pF^r"\l  i',  M«)g<'">(7,n[f  ((r?S*^)"-i-log  2)- 

- ! (r,r-r,f¥loghS*^-r,n 

p-i  J-i  (Vo-Vj  ) J 


(156) 


If  4 is  a positive  even  number  then  we  put 

fl  = 2a,  (157) 

where  a is  a positive  number,  and  formula  (156)  follows  directly  from  case  (i)  above. 

If  4 is  a positive  odd  number  then  we  put 

q = 2a-\,  (158) 

where  a is  a positive  number.  Formula  (156)  then  follows  immediately  from  case  (ii)  above  when  M is  any  odd 
positive  integer  and  /,  2=s/«Af— lis  any  even  positive  integer,  from  case  (iii)  above  when  M is  any  even 
positive  integer  and  /,  2 « / « A/  is  any  even  positive  integer,  from  case  (iv)  above  when  M is  any  even  positive 
integer  and /,  l«/«Af-l  is  any  odd  positive  integer  and  from  case  (v) above  when  M is  any  odd  positive 
integer  and  /,  1 ^ ^ M,  is  any  odd  positive  integer. 

The  formula  (156)  is  precisely  the  formula  derived  by  Garner  and  Fox*  and  now  it  has  been  demonstrated 
that  the  numerical  procedure  converges  and  provides  a value  for  f4.,(jr,  yj^;  v.  Mo)  which  is  correct  when  q 
becomes  indefinitely  large.  Garner  and  Fox*  considered  quasi  steady  flow  only  but  their  formula  is  valid  for 
general  frequencies  of  harmonic  oscillation  and  was  applied  to  this  case  by  Long’.  For  a given  value  of  q the 
corresponding  value  of  a is  obtained  from  (157)  when  q is  an  even  positive  integer  and  from  (158)  when  q is  an 
odd  positive  integer.  The  number  of  integration  points  A is  determined  then  from  (149)  when  q is  an  even 
positive  integer  and  from  (150),  (151),  (152)  or  (153)  when  q is  an  odd  positive  integer. 

The  value  of  q used  in  evaluating  UrAxu'^,  y/*^; Moo)  from  formula  (156)  can  be  taken  individually  for 
each  combination  (r,  s,  /,  /),  that  occurs,  as  a value  which  gives  adequate  accuracy  for  U,Ax^u*^,  yj^^',  v,  Mx). 
It  is  to  be  observed  that  a higher  value  of  q is  needed  for  a given  accuracy  in  UrAxu'^,  yj^^',  v,  Mao)  as  the 
point  approaches  the  leading  edge  or  trailing  edge  of  the  planform  than  is  otherwise  needed.  In  our  work 

that  follows  we  shall  take  q to  depend  only  on  I and  denote  it  by  qi.  This  will  enable  us  to  take  higher  values  of  q 
for  points  near  the  leading  or  trailing  edges  of  the  planform  than  for  others  further  away  from  the 
leading  or  trailing  edges,  if  we  so  desire,  without  unduly  complicating  the  arrangement  of  the  calculation,  but 
we  note  that  a more  general  variation  of  q could  be  taken. 

We  would  like  to  mention  here  that  Zandbergen,  Labrujere  and  Wouters  in  Ref.  8 take 


A = a(A/+l)-l  (159) 

where  a is  any  positive  integer,  and  use  the  weight  function  Vl  -tjo  for  the  numerical  integration  in  formula 
(121).  Provided  that  condition  (130)  holds,  their  final  formula  is  our  formula  (133),  but,  as  we  have  indicated, 
this  is  precisely  equivalent  to  (144)  if  condition  (141)  holds. 

In  this  case 


and  correspondingly  the  formula  (133)  must  likewise  be  modified.  The  second  derivative 
[(a^/di7o)^r"’(f  »?J^,  Vo; as  well  as  the  first  derivative  [(a/3i7o)4"’(f,  Vj^;  Vo;  v, 
appear  in  these  modified  formulae,  in  contrast  to  formula  (156),  where  neither  derivative  appears.  That  the 
second  derivative  appears  for  all  /=  1, 2, . . . , A/,  is  a direct  result  of  A having  been  given  the  formula  (159), 
and  it  is  only  when  A has  this  formula  that  the  second  derivative  appears  for  all  /. 

The  modified  formula  (133)  would  be  exact  with  the  A integration  points  given  by  formula  (159)  if  the 
function  (119)  were  a polynomial  of  degree  2fl(Af+l)— w— 2 in  vo-  The  formula  requires  the  values  of 
Vj^,  Vo;  V,  Mco)  at  the  a(A/+ 1)- 1 integration  points  t?o  = p = 1, 2, ....  A,  together  with  the 
values  [(aVa77o)/^"’(f  9o; Mao)]„o-nT  [(3/39o)4"\f,  »?/  , r?o;  J',  Af«.)]no-V}‘’.  a{M+\)+\ 

values  in  all.  

On  the  other  hand,  if  we  take  A to  be  given  by  (149)  and  use  the  weight  function  1/v1-tj5  for  the  numerical 
integration  in  formula  (121),  then  the  formula  (154)  is  valid  and  would  be  exact  if  the  function  (1 19)  were  a 
polynomial  of  degree  «2a(A/+  l)-m  - 2 in  vo-  The  formula  requires  the  values  of  4"'(f  17  v.  Moo), 

p = 1, 2, . . . , A,  where  is  given  in  formula  (129),  and  /i"*(f,  vf*^>  vf*^;  Af*),  i.e.  a{M  + 1)+ 1 values  in 

all.  We  conclude  that  formula  (154)  for  case  (i)  above  is  of  comparable  accuracy  with  the  modified  formula 
(133)  for  the  same  value  of  a,  although  less  work  is  required  in  the  evaluation  of  (154)  than  in  the  evaluation  of 
the  modified  formula  (133).  All  the  same,  the  modified  formula  (133)  is  somewhat  better  conditioned  than 
formula  (154)  because  the  absolute  values  of  the  quantities 

(164) 

attain  greater  values  in  (154)  for  a given  value  of  a than  they  do  in  the  modified  formula  (133). 

We  note  further  that  formula  (154)  for  case  (i)  above  with  a = 1 is  the  formula  used  by  Hsu'^  Also  the 
modified  formula  (133)  with  weight  function  Vl  - tjJ  and  with  a = 1 in  the  formula  (159)  for  determining  A is 
the  formula  suggested  by  Multhopp  in  Appendix  V of  Ref.  1. 


l-(0' 


24 


¥ 


T 


If  A/  = m and  ^ = 1 so  that  M = M=m,  then  formula  (156)  is  the  formula  for  Ur,,(x,  yf*^;  v.  A/®)  which  is 
normally  used  in  lifting  surface  theory.  The  corresponding  value  of  A and  the  weight  function  for  the 
integration  formula  (121)  are  those  appropriate  to  case  (ii),  (iii),  (iv)  or  (v)  above,  whichever  has  to  be  used,  and 
this  depends  on  whether  the  numbers  A/  and  / are  even  or  odd.  The  numerical  integration  formulae  (121)  and 
(131)  are  not  of  good  accuracy  in  this  case  because  the  value  of  A given  by  any  one  of  the  formulae  (150),  (151), 
(152)  or  (153)  is  not  large  enough.  The  numerical  integration  formulae  (134)  and  (135),  on  the  other  hand,  are 
precise.  Consequently  the  formula  (156)  may  not  give  good  accuracy  for  t4,,(x,  yi"*;  v,  Af*)  in  this  case.  A 
significant  improvement  in  accuracy  for  UrAx,  yj^\  v,  Af*)  should  be  obtained  if  q is  increased  to  2,  and  this  is 
equivalent  to  using  Hsu’s  formula'^  By  taking  the  value  of  q high  enough  we  can  obtain  t/„(x,  yi"';  v,  MA  to 
any  accuracy  that  we  like. 

If  now  we  take 


q=qi,  7=1,2,  ...,fV, 


and 


M,  = q,(M+\)-l,  /=1,2,...,N, 

then,  analogously  to  formula  (156)  we  have  the  formula 
Tf  \jf  \ 

y y j y |/y  A^oo) 


=4  f 

S ps=l 


(165) 


(166) 


X I V,"'(iT',  vr-,  v,  M^)g^r\vr)x 

j-l  (Vo-Vj  ) s 

X [f  ((^ry-^-iog  2)-f  wr-  log  i^r- x 


(Vo-v'Tf  ‘‘H 


where 

iT>=i-^r 

Quite  analogously  to  the  derivation  of  formula  (167)  we  derive  the  numerical  formula 

rr  /=WM)  (M). 

^ _(M)  _(A/,).  „ w N„(m)  v 

2 iVJ  iVp  1 Aiooj^m— 5+1  Wp 


(167) 


(168) 


(vo-vr^^ 

+ p/=^">(fT\  vr-,  A7„)g|:i...(r,nx 
x[f((T,S'^^-i-log2)-?^(,,r-r,r'yioghS^>-r,r'’|x 


f'"' gir^'WKi -’Jo  , 1 

L,  (v,-Wh‘'  ‘‘n 


(169) 
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Then,  on  making  use  of  (42),  we  get  from  (167)  and  (169) 


u"’,  yi"';  M«)+  Kt/,,„-.*,(xi;^*'\  ir-,  V,  A/oc) 

= 4 x'  (if‘.'"*(»?p"'V^«'."’*(-Tjr'V,"*(fT’,  np"'*;  Woo)x 

5 p»  I 

J-i  iVo-Vj  ) * 

x[^((T,ry-i-log2)-  ? (r,r’-r,r'f  log|r,r-T,r''|x 

12  p«i 

X (170) 

J-1  (Vo~Vj  ) J 

which  is  the  formula  we  use  in  getting  the  elements  (112)  of  the  matrix  [(/)  appearing  in  equation  (1 10). 
We  note  that 

?((„r)^-^-iog2)-  ^ (vr-v^'fiogivr-vr^x 

2 p-i 

(171) 

J-i  (T/o-Tjy  ) 

appears  to  converge  to  zero  as  qi^<x),  as  far  as  one  can  judge  from  numerical  results,  but  the  present  author 
has  not  succeeded  in  proving  this  analytically.  If  this  is  true  then  we  could  miss  out  the  coefficient  of 
A^oo)  in  formula  (170)  and  still  get  convergence  to  the  correct  value  as  q,  -» c»,  but  we  must 
expect  the  convergence  to  be  slower  than  that  of  formula  (170). 

To  apply  formula  (170)  we  must  still  evaluate  Vj^\  v,  M^)  numerically  from  (1 14)  and  the 

process  for  doing  this  is  described  in  Appendix  C,  where  the  use  of  Chebyshev  polynomials  is  recommended 
and  illustrated. 


Examples 

We  give,  in  this  section,  a selection  of  results  of  calculations  carried  out  on  four  planforms.  The 
approximations  (5,*  to  the  generalised  airforce  coefficients  O/k  (j',  M®)  are  obtained  from  formula  (84)  and  the 
approximations  fk(.Xo,  yo)  to  the  loading  functions  /t(xo,  yo;  v,  M<x>)  in  the  mode  k of  oscillation  at  points 
(xo,  yo)  on  the  wing  planform  are  obtained  from  formula  (98). 

The  four  planforms  considered  are  a tapered  swept  wing  of  aspect  ratio  2,  a tapered  swept  wing  of  aspect 
ratio  6,  and  rectangular  wings  of  aspect  ratios  2 and  8.  Diagrams  of  these  planforms  are  given  in  Figs.  2, 3, 4a 
and  4b. 

The  tapered  swept  wings  have  trailing-edge  motivators  as  shown  and  each  is  taken  to  oscillate  in  rotation 
about  the  hinge  line  at  its  leading  edge.  The  loading  in  this  mode  of  oscillation  has  a logarithmic  singularity  at 
the  hinge  line,  but,  despite  this,  an  approximation  to  the  loading  in  the  form  (22)  is  permissible,  except  in  the 
immediate  neighbourhood  of  the  hinge  line,  and  the  corresponding  approximations  to  the  generalised  airforce 
coefficients  obtained  from  formula  (34)  should  be  acceptable.  The-approximations  (75)  and  (76)  are  not  valid 
for  the  motivator  mode  so  that  we  must  use  tlie  exact  expressions  (66)  and  (67)  together  with  the  associated 
formulae  (60)  and  (36). 


4.1.  Tapered  Swept  Wing  of  Aspect  Ratio  2 

The  planform  of  this  wing  is  illustrated  in  Fig.  2.  The  x coordinate  XL(y)  of  the  leading  edge  at  spanwise 
position  y is  given  by  the  formula 

XL(y)  = '«^|y|  -s«y«s,  (172) 
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and  the  chord  length  c(y)  at  spanwise  position  y is  given  by  the  formula 


c(y)  = ^(2^3 + 3)-j(4V3- 2) 


(173) 


The  leading  edge  and  trailing  edge  of  the  wing  have  a discontinuity  of  slope  at  the  centre  line  of  the  wing  and, 
in  order  to  enhance  the  convergence  of  the  numerical  results,  it  is  necessary  to  change  the  shape  of  the  leading 
and  trailing  edges  of  the  wing  so  that  there  are  no  discontinuities  of  slope.  This  change  of  shape  is  made  in  the 
region  -yR«y  =£yR,  where  0<yR<s,  and  is  known  as  rounding  of  the  leading  and  trailing  edges  in  the 
neighbourhood  of  their  central  portions.  It  is  desirable  to  arrange  that  XtCy ) and  c”(y)  are  continuous,  as  well 
as  Jt/.(y),  x'tiy),  c(y),  c'(y),  because  a more  accurate  numerical  estimate  for  y ; if,  Moc)of  formula  (113) 
is  then  obtained.  The  change  of  shape  is  achieved  by  taking 


( -yR^y^ya 

^ ('^yR|A|  yR«|y|«s 


(174) 


and 


where 


-(2V3  + 3)-^(4>/3-2)/(A)  -yR«y«yR 

J(2V3  + 3)-j(4V3-2)lAl  yR^lyl^i 


(175) 


A=—  (176) 

yR 

and  /(A ) is  an  arbitrary  even  function  of  A with  continuous  second  derivative  f'(A ) in  - 1 « A « 1 which  is  such 
that 


/(1)=1, 

(177) 

Ai)=i. 

(178) 

r(i)=o. 

(179) 

The  origin  of  coordinates  at  the  apex  of  the  tapered  wing  is,  in  general,  not  on  the  rounded  leading  edge  of  the 
modified  wing. 

Here  we  shall  follow  Hewitt"*  and  take 


/(A)— i6  + i6A^~l^A^  + i6A*’.  (180) 

Further  we  shall  take 

^ = sin  = 0 1950903  (181) 

s V16/ 

and  we  shall  take  the  typical  length  / to  be  the  geometric  mean  chord  c of  the  planform,  i.e., 

l = c = s.  (182) 
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The  motivators  shown  in  Fig.  2 are  known  when  the  coordinates  (A^i,  Ti)  and  (X2,  Tj)  of  the  inboard  and 
outboard  extremities  respectively  on  the  leading  edge  of  the  motivator  on  the  starboard  side  of  the  wing  are 
known.  These  coordinates  are  given  by 


A,=i(3-t-2>/3)c, 

(183) 

II 

(184) 

Ar2  = i(9-l-4v^)c, 

(185) 

Y2  = C. 

(186) 

All  the  edges  of  the  motivators  are  straight  and  the  port  motivator  is  the  mirror  image  of  the  starboard 
motivator  in  the  plane  of  symmetry  of  the  wing  planform. 

The  leading  edge  of  the  starboard  motivator  has  the  equation 

x=x„(y)  y.^y^Tj  (187) 

and  the  leading  edge  of  the  port  motivator  has  the  equation 

(188) 

The  explicit  form  of  the  function  Xniy)  may  be  obtained  from  the  information  given  above. 

Three  modes  of  oscillation  are  considered.  These  are  specified  by  giving  the  functions  ik(x,  y),  introduced  in 
formula  (6),  for  A:  = 1, 2 and  3.  These  functions  are  taken  in  this  example  to  be 


6(-*,y)=l. 

(189) 

II 

1 H 

(190) 

and 


{3(x,y)={'‘~-^)ff(x-x„(y)XH(y-Y0  + H(-y-y>)-ff(y-y2)-H(-y-y2)]  (191) 

where  H(x)  is  Heaviside’s  unit  function 


x<0 

x>0. 


(192) 


Approximations  (5i/  to  the  generalised  airforce  coefficients  have  been  evaluated  for  /=  1,2,3; 

/■  = 1, 2, 3,  when 


»'  = 0-32506  and  Afct.  = 0-78060. 


We  write  Oi/  in  the  form 


(193) 


where  O'l/  and  are  real  quantities. 

The  numerical  values  of  O!,,  <5j,  / = 1,  2, 3;  / = 1,  2,  3,  obtained  with  m=M  and 

q>  = q,  /=1,2 N,  (194) 

for  a selection  of  values  of  the  parameters  m,  n,  N and  q are  given  in  Table  1 . Examination  of  Table  1 reveals 
that  the  Oti  and  OJi  need  a fairly  high  value  of  q in  their  evaluation  for  them  to  be  reliable  for  a given  set  of 
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parameters  m,  n and  N.  As  n and  N are  increased  with  m = M = 1 5 the  values  O'a  and  <5,y  seem  to  converge. 
This  behaviour  is  illustrated  graphically  in  Figs.  5a  to  5r.  The  convergence  with  increasing  q and  increasing  n 
and  N appears  to  be  best  for  the  motivator  mode  of  oscillation,  mode  3. 

Approximations  /*  (x,  y ) to  the  loading  also  have  been  evaluated  for  ^ = 3 for  the  above  values  of  t'  and  A/*. 
We  write  fk(x,  y)  in  the  form 

tk(x,y)=^f'k(x,y)+ifk{x,y)  (195) 

where  y)  and  fk(x,v)  are  real  o^uantities. 

Numerical  values  of  /^(x,  y)  and  fjix,  y)  are  given  in  Table  2 for  m = 15,  n = 10,  A/=  15,  JV=  10,  q/  = 8, 
/=  1, 2, . . . , 10.  The  coordinates  (x,  y)  have  been  transformed  to  the  coordinates  ((,  tj)  by  means  of  the 
transformation  (20).  The  loadings  f'jix,  y)  and  f^x,  y)  are  given  at  the  set  of  chordwise  lines 

T}  = Tjy  = cos|^^-j^|  /=1.2 8.  (196) 

The  locations  of  the  hinge  4 = ^h(v)  3t  the  chordwise  lines  within  the  motivator  span  are  given  by  the  following 
set  of  numbers; 


i 


/=  4 5 6 7 8 

7jy  = 0-55557  0-70711  0-83147  0-92388  0-98079 

fH(^/)  = 0-74654  0-73387  0-71750  0-69853  0-68158 


Curves  of  t3(x,  y)and  f^x,  y)  for  a range  of  values  of  ^ along  each  chordwise  line  rj  = tj„/=  1, 2, . . . , 8,  are 
given  in  Figs.  6a  and  6b.  The  curve  for  each  value  of  17  has  its  own  origin,  which  is  marked. 

It  is  to  be  noticed  that  f'jix,  y)  along  a chordwise  line  within  the  motivator  span  has  a deep  minimum  near  to 
the  hinge  line.  The  actual  theoretical  real  part  of  the  loading  becomes  infinite  like  log  |^-fH(»7)|  at  the  hinge 
line.  Curves  for  both  f'jix,  y)  and  /'^(x,  y)  have  a number  of  undulations  along  them  and  this  is  due  to  the 
truncation  of  an  infinite  series,  implied  in  formula  (22),  for  a function  which  has  a logarithmic  singularity  at  the 
hinge  line. 


4.2.  Tapered  Swept  Wing  of  Aspect  Ratio  6 

The  planform  of  this  wing  is  illustrated  in  Fig.  3.  The  x coordinate  XL(y)  of  the  leading  edge  at  spanwise 
position  y is  given  by  the  formula 


-ti.(y)  = 3('^  + 2)|y|  -J«y«s,  (197) 

and  the  chord  length  c(y)  at  spanwise  position  y is  given  by  the  formula 

c(y)  = -3|y|+^^  -^«y«s.  (198) 


Here  again,  as  with  the  tapered  swept  wing  of  aspect  ratio  2 we  round  the  leading  and  trailing  edges  in  the 
neighbourhood  of  their  central  portions.  The  functions  Xl  (y ) and  c (y ) above  are  thereby  modified  to  the  forms 


and 


U('/3+|)yR|A|  yR«|y|=ss, 


, . f2J-3yR/(A)  -yR^y^yn 

Ui-sy/ilAl  y„«|y|«j. 


(199) 


(200) 


where  A is  given  by  formula  (176). 
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In  this  example  we  shall  follow  Zandbergen,  Labrujere  and  Wouters^  and  take 

/(A)  = i + A'-i|Ap. 


(201) 


We  shall  again  take  yu  to  be  given  by  formula  (181)  and  the  typical  length  I to  be  the  geometric  chord  c of 
the  planform,  i.e.. 


I = C = 3S. 


(202) 


The  motivators  shown  in  Fig.  2 are  known  when  the  coordinates  (Xi,  Yi)  and  (X2,  Y2)  of  the  inboard  and 
outboard  extremities  respectively  on  the  leading  edge  of  the  motivator  on  the  starboard  side  of  the  wing  are 
known.  These  coordinates  are  given  by 


^ (2^5  97  \. 


^ /7n/3  91 

V 10  ioo)‘^’ 


Y-^ 

10- 


(203) 

(204) 

(205) 

(206) 


All  the  edges  of  the  motivators  are  straight  and  the  port  motivator  is  the  mirror  image  of  the  starboard 
motivator  in  the  plane  of  symmetry  of  the  wing  planform.  The  ratio  of  motivator  chord  to  wing  chord,  in  this 
example,  is  independent  of  the  spanwise  location  of  the  chord  and  has  the  value  0-3. 

Three  modes  of  oscillation  are  considered  and  these  are  again  defined  by  equations  (1 89),  (190)  and  (191 ). 

Approximations  to  the  generalised  airforce  coefficients  (?„(»',  M)  have  been  evaluated  for  /=  1,  2, 3; 
j=  1,2,2,  when 

j^  = 3-1569  and  Moo  = 0-4. 

We  write  (5,,  in  the  form  (193).  The  numerical  values  of  O'a,  Oh,  i = 1,2,2,  j = 1,2,  3,  obtained  for  a selection 
of  values  of  the  parameters  m,  n,  M,  N and  qi,  I = 1,2, ..  .,N,  are  given  in  Table  3.  For  all  evaluations  we  have 
again  taken  all  the  q,  equal  as  in  formula  (194). 

Results  for  On,  / = 1,  2;;  = 1, 2,for  this  same  tapered  swept  wing  of  aspect  ratio  6 with  the  same  rounding  of 
the  leading  and  trailing  edges  and  the  same  frequency  parameter  and  Mach  number  have  been  presented  by 
Lehrian  and  Garner’.  Lehrian  and  Garner’s  results  have  to  be  multiplied  by  -6  to  make  them  compatible  with 
the  present  results  because  of  their  different  non-dimensionalising  factor  in  the  definition  of  On.  Their  results, 
after  multiplication  by  -6,  are  presented  in  Table  4. 

Lehrian  and  Garner’  use  an  extension  to  oscillatory  flow  of  the  steady  flow  method  of  Zandbergen, 
Labrujere  and  Wouters*.  In  their  development  of  this  method,  the  number  of  upwash  points  is  equal  to  the 
number  of  basic  loading  functions,  so  that  M = m and  N=  n.  The  parameter  a determines  the  number  of 
spanwise  integration  stations  to  be  used  in  the  numerical  evaluation  of  the  spanwise  integral  in  the  integral 
equation  much  as  the  parameter  q determines  the  number  of  spanwise  integration  stations  to  be  used  in  the 
present  method  for  the  numerical  evaluation  of,  effectively,  the  same  spanwise  integral.  The  amount  of 
computation  necessary  to  get  results  with  q = 2a  in  the  present  method  is  comparable  with  that  required  in  the 
method  of  Lehrian  and  Garner  and  the  results  also  may  be  expected  to  be  of  comparable  accuracy,  and  this  is 
borne  out  by  examination  of  the  respective  results  obtained.  Lehrian  and  Garner  used  effectively  the  same 
values  (66)  and  (60)  to  obtain  (i-,,,  as  were  used  to  get  the  present  results  but  they  used  values  of  0*.^, 
corresponding  to  the  approximations  (76)  rather  than  those  obtained  from  (67)  and  (36).  The  differences  in 
results  arising  from  this  cause  are  expected  to  be  small. 

In  the  present  method  we  may  take  M^m  and  N'^n.  The  effect  of  taking  M>m,  N = n can  be  seen  by 
examining  Table  3.  The  comparatively  large  differences  in  the  values  of  0^  for  m = 14  and  m = 15  when 
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M=m  are  considerably  reduced  when  M>m.  The  differences  in  values  of  On  for  m=22  and  m=23  when 
M=  m are  not  so  large  as  those  for  m = l4  and  m = 15,  but  these  differences  also  are  reduced  when  M>  w. 
The  differences  in  values  of  On  for  w = 30  and  m = 31  when  M=m  are  again  smaller  than  are  those  for 
m=22  and  m = 23.  This  is  an  indication  that  the  results  are  converging  as  m increases,  and  that  the 
convergence  is  more  rapid  when  M>m. 


4.3.  Rectangular  Wings  of  Aspect  Ratios  2 and  8 

The  planforms  of  these  two  wings  are  illustrated  in  Figs.  4a  and  4b.  The  origin  of  coordinates  is  taken  to  be  at 
the  middle  point  of  the  leading  edge.  The  typical  length  / is  taken  to  be  the  wing  chord  c.  Two  modes  of 
oscillation,  heave  and  pitch  about  the  leading  edge,  are  considered.  These  are  specified  by  the  functions 
{k(x,  y)  introduced  in  formula  (6)  and  defined  for  k = 1 and  2 respectively  in  equations  (189)  and  (190),  with 
c = c. 

Approximations  0,,  to  the  generalised  airforce  coefficients  Oni^Mx)  have  been  evaluated  for  /=  1,2; 
/ = 1, 2 when 

j/=10,  M*  = 0-8  and  q/  = 32,  /=  1, 2, . . . , N, 


for  a selection  of  values  of  the  parameters  m,  n,  M,  N.  We  write  On  in  the  form  (193).  The  numerical  values 
obtained  for  O'n,  Oh,  i = 1. 2;  / = 1, 2,  are  given  in  Table  5 for  the  rectangular  wing  of  aspect  ratio  2 and  in 
Table  6 for  the  rectangular  wing  of  aspect  ratio  8.  For  both  sets  of  results  there  is  evidence  of  convergence  as  m, 
n,  M and  N are  increased,  although  it  would  appear  that  the  convergence  is  more  rapid  with  odd  values  of  m 
than  with  even  values  of  m.  We  shall  consider  the  values  obtained  for  m = 19,  n = 8,  A/=  19,  8 to  be  the 

best  estimates  of  the  values  of  O'n,  Oh, « = 1,  2;  / = 1, 2.  Values  of  O'n,  Oh/  are  in  close  agreement  with  these 
best  estimates  for  smaller  values  of  m,  n,  M and  N.  To  quantify  the  closeness  of  agreement  of  On  with  (?„  we 
introduce  a measure  en  of  percentage  difference  by  means  of  the  formula 


en 


= 100x 


{o'iif+Aohf  J 


(207) 


and  the  arithmetic  mean  e of  over  / = 1, 2;  y = 1, 2, 

e =4(eii  + ei2  + e2i  + «22).  (208) 


A selection  of  values  en  for  the  wings  of  aspect  ratios  2 and  8 are  given  in  Tables  7 and  8 respectively. 

The  values  of  On  are  all  in  poor  agreement  with  the  best  estimates  On  when  n = 2.  For  m = 4, « = 2,  Af  = 4, 
Af  = 2,  the  mean  p is  10-9  for  the  wing  of  aspect  ratio  2 and  101  for  the  wing  of  aspect  ratio  8.  The  percentage 
differences  do  not  change  greatly  when  m,  M and  N are  changed  with  n kept  at  2,  even  when  m,  M and  N are 
all  increased  substantially. 

There  is  a distinct  improvement  in  the  values  of  On  when  n is  increased  to  4.  For  m = 4,  n = 4,  Af  = 4,  N = 4, 
the  mean  e is  0 042  for  the  wing  of  aspect  ratio  2 and  018  for  the  wing  of  aspect  ratio  8.  For  m = 4,  n = 4, 
Af  = 19,  N = 8 the  mean  e is  0-032  for  the  wing  of  aspect  ratio  2 and  0-27  for  the  wing  of  aspect  ratio  8.  Thus 
changing  (A/,  N)  from  (4, 4)  to  (1 9, 8)  with  (m,  n ) kept  at  (4, 4)  has  caused  a small  improvement  in  the  values  of 
On  for  the  wing  aspect  ratio  2 and  a deterioration  in  these  values  for  the  wing  of  aspect  ratio  8.  This  latter 
suggests  that  m = 4 is  not  large  enough  for  the  wing  of  aspect  ratio  8. 

For  m=9,  n=4,  M = 9,N  = 4 the  mean  e is  0-022  for  the  wing  of  aspect  ratio  2 and  0- 10  for  the  wing  of 
aspect  ratio  8.  For  m = 9,  n = 4,  A/=  19,  N = 8 the  mean  e is  0-018  for  the  wing  of  aspect  ratio  2 and  0-08  for 
the  wing  of  aspect  ratio  8.  Thus  changing  (Af,  N)  from  (9, 4)  to  (19, 8)  with  (m,  n)  kept  at  (9, 4)  has  caused 
merely  a marginal  improvement  in  the  values  of  On-  There  are  no  further  distinct  improvements  in  increasing 
m to  14  or  to  19  while  keeping  n=4. 

Further  improvements  in  the  values  of  On  are  obtained  when  n is  increased  to  6.  For  m=4  the  percentage 
differences  for  n = 4 and  n = 6 are  of  the  same  order,  but  for  higher  values  of  m the  percentage  differences  are 
lowerforn  =6than  they  are  forn  =4whenA/=  19. Form  =9,/i  =6,  Af  = 9,  V=6themeane  is0-032forthe 
wing  of  aspect  ratio  2 and  0-12  for  the  wing  of  aspect  ratio  8.  For  w = 9,  n = 6,  A/  = 19,  ^ = 8 the  mean  e is 
0-001  for  the  wing  of  aspect  ratio  2 and  0-000  for  the  wing  of  aspect  ratio  8.  Thus  changing  (A/,  N)  from  (9, 6) 
to  ( 1 9, 8)  with  (w,  n ) kept  at  (9, 6)  has  produced  almost  complete  agreement  of  the  On  with  the  best  estimates. 

For  w = 14,«=6,  A/=14,N  = 6the  mean  c is  0-  0 1 5 for  the  wing  of  aspect  ratio  2 and  0-038  for  the  wing  of 
aspect  ratio  8.  For  m = 14, « = 6,  A/=  19,  Af  = 18  the  mean  e is  0-014  for  the  wing  of  aspect  ratio  2 and  0-037 
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for  the  wing  of  aspect  ratio  8.  Thus  changing  (M,  N)  from  (14, 6)  to  (19, 8)  with  (m,  n)  kept  at  (14, 6)  has 
practically  no  effect  on  the  values  of  On-  We  may  note  that  almost  complete  agreement  of  the  On  with  the  best  i 

estimates  has  not  been  achieved  with  (m,  «)  kept  at  (1 4, 6),  but,  all  the  same,  the  percentage  differences  c/,  are  i 

very  small.  i 

It  is  as  well  to  remark  here  that  our  results  have  been  obtained  with  o;  = 32,  / = 1 , 2, . . . , Af,  whereas  results  : 

* ^ 

from  earlier  theories,  such  as  those  of  Refs.  3, 4 and  5,  are  comparable  with  the  results  we  would  get  by  taking  j 

<?/=!,  /=  1, 2, . . . , N;  M=m  and  N=n.  For  m=M  = A,  N=n=A  and  = 1,  /=  1, 2, . . . , N,  we  get  i 

Oii  = 0-84678-i3-2052,  (5,2  = -3-2858-/3  1810,  (52i  = 0-90492 -(0-83073,  (^22  = -0-51381 -/2  0731  I 

for  the  wing  of  aspect  ratio  2 and  (5ii  = -1  1040-/13-627,  (5,2  = — 16-484-/7-6979,  (^21=  ] 

1-7608-/4-5769,  (^22  = -4-5283-/6-2760  for  the  wing  of  aspect  ratio  8.  The  corresponding  percentage 

differences  en  are  £,,  = 2-5,  £,2  = 3-1,  £2,  = 5-0,  £22  = 5-4  for  the  wing  of  aspect  ratio  2 and  £,,  = 16-6,  i 

£,2=  17-7,  £2,  = 21-7,  £22  = 26-3  for  the  wing  of  aspect  ratio  8.  The  magnitudes  of  the  £,,  are  of  the  order  one  i 

hundred  times  the  corresponding  values  in  Tables  7 and  8 where  the  qi  are  changed  to  ^/  = 32,  / = 1 , 2, . . . , A/.  - j 

To  get  Sii  of  the  same  order  as  those  in  Tables  7 and  8 with  <7/  = 1,  /=  1,  2, . . . , /V,  M=m  and  /V=  n an  j 

undesirably  large  value  of  m would  be  needed  when  n = 4.  The  values  of  the  On  with  <7,  = 1,  / = 1, 2, . . . , N, 

m = M = 4,  N = n =A  are  quite  unacceptable  for  the  aspect  ratio  8 wing.  j 

To  appraise  the  convergence  of  the  results  as  M and  N are  increased  we  examine  in  Tables  5 and  6 the  values  j 

of  O'n,  O'n  when  m = 4,  n = 4.  For  the  wing  of  aspect  ratio  2 the  convergence  of  all  the  O'n,  On  is  rapid  for  N i 

fixed  and  M increasing  through  4, 9, 14, 19.  Then  with  M held  at  19  the  convergence  is  rapid  for  N increasing  . 

through  4,  6,  8.  For  the  wing  of  aspect  ratio  8 convergence  of  these  quantities  is  on  the  whole  slower.  Thus,  I 

while  0\\,  Ou  continue  to  converge  quite  rapidly,  the  convergence  of  0'22,  Oh  for  A/ fixed  and  M increasing  i 

through  4,  9,  14,  19  is  less  rapid,  and  to  achieve  convergence  to  five  significant  figures  the  value  of  M would  i 

have  to  be  increased  beyond  19  when  N = 4,  6 and  8.  Convergence  to  four  figures,  however,  has  practically  j 

been  achieved  at  M = 19,  and  also  at  A/=  8 when  M is  held  at  19.  The  behaviour  of  the  other  O'n,  Olj  is  ] 

intermediate  between  the  behaviours  of  O'n,  Oh  and  0'22,  Oh-  j 

The  same  pattern  of  convergence  may  be  observed  with  other  values  of  m and  n kept  fixed  while  M and  N ] 

are  increased  although  there  are  more  numbers  available  in  any  convergence  sequence  for  the  lower  values  of 
m and  n.  With  M and  N being  increased  while  m and  n are  fixed  there  are  converged  values  of  O'n,  On  1 

appropriate  to  m spanwise  loading  functions  and  n chordwise  loading  functions.  Convergence  with  respect  to  j 

m and  n being  increased  then  leads  to  the  best  estimates  of  these  values.  J 

We  may  conclude,  for  both  rectangular  wings,  that  increasing  n from  n = 2 effects  a substantial  improve-  i 

ment  in  the  results  for  O'n,  On,  / = 1, 2;  j = 1, 2,  even  while  keeping  m fixed  at  m = 4,  whereas  increasing  m j 

from  w = 4 and  keeping  n fixed  at  n = 2 hardly  effects  any  improvement  at  all.  For  practical  purposes  values  of  j 

O'n,  On,  / = 1, 2;/  = 1, 2,  can  be  obtained  to  sufficient  accuracy  with  n = 4 for  both  wings  and  respectively  with 
m = 4 and  m = 9 for  the  wings  of  aspect  ratios  2 and  8. 

4.4.  Discussion 

Results  for  the  generalised  airforce  coefficients  on  a tapered  swept  wing  of  aspect  ratio  2 show  that  a more 
refined  numerical  integration  of  the  spanwise  integral  than  the  one  used  by  Multhopp'  must  be  used  in  order  to 
obtain  acceptable  results.  Values  of  M and  N equal  to  m and  n respectively  were  used  for  the  results  on  the 
tapered  swept  wing  of  aspect  ratio  2,  but  for  a tapered  swept  wing  of  aspect  ratio  6 values  of  M>m  were 
considered.  For  the  tapered  swept  wing  of  aspect  ratio  6 values  obtained  with  M^^lm  showed  better 
convergence  on  the  whole  than  did  those  with  M=m.  In  particular,  the  values  of  the  generalised  airforce 
coefficients  for  m = 15,  n =6,  A/=  30,  V = 6 are,  on  the  whole,  closer  to  the  values  for  m = 30, « =6,  Af  = 30,  j 

A/  = 6 than  are  those  for  m = 15,  n =6,  M=  15,  N = 6. 

A comprehensive  set  of  results  for  rectangular  wings,  of  aspect  ratios  2 and  8 respectively,  illustrate  the 
nature  of  convergence  of  values  of  the  generalised  airforce  coefficients  when  m,  n,  M and  N are  increased. 

Generally,  the  effect  of  increasing  the  values  of  M and  A7  is  to  give  convergence  of  results  for  given  m and  n, 
and  then  the  effect  of  increasing  the  values  of  m and  n is  to  give  the  final  converged  results. 


5.  Concluding  Remarks 

In  the  lifting  surface  theory  developed  in  this  paper  the  loading  is  represented  approximately  as  a linear 
superposition  of  a finite  number  of  known  linearly  independent  elementary  functions,  there  being  a 
combination  of  n chordwise  by  m spanwise  of  them,  as  in  a number  of  other  theories.  For  given  n and  m,  there 
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are  N ^ n chordwise  by  Af  ^ m spanwise  integration  points  at  which  the  upwash  is  evaluated.  The  mathemati- 
cal analysis  leads  to  the  same  kind  of  refined  process  of  numerical  evaluation  of  the  spanwise  integral  in  the 
integral  equation  as  was  introduced  by  Garner  and  Fox  in  Ref.  6. 

The  results  obtained  for  generalised  airforce  coefficients,  using  the  Gamer  and  Fox'*  refined  process  of 
numerical  evaluation  of  the  spanwise  integral  in  a lifting  surface  theory  are  a considerable  improvement  over 
those  obtained  by  a simple  extension  of  Multhopp’s  steady  flow  method  to  oscillatory  flow.  The  results 
obtained  with  N = n and  M = m aiexi  good,  particularly  for  rectangular  wings,  that  it  is  not  easy  to  detect  any 
further  improvement  obtained  by  taking  N > « and  M>m.  Nethertheless  results  for  the  example  of  a tapered 
swept  wing  aspect  ratio  6 do  show  evidence  of  a further  improvement  being  obtained  when  M>m.  Thus,  for  a 
given  M,  if  one  allows  m<M  it  may  be  possible  to  use  a lower  value  of  m than  if  one  insists  on  having  m=M, 
! and  get  results  of  comparable  accuracy,  with,  incidentally,  less  storage  space  and  execution  time  on  the 

computing  machine. 
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LIST  OF  SYMBOLS 

Speed  of  sound  in  undisturbed  main  flow.  Also  a positive  integer  which  deter- 
mines A from  one  of  (149),  (150),  (151),  (152)  or  (153)  | 

Defined  by  formula  (A- 19) 

The  integral  (A-75).  Also  a point  within  (0,  oo) 

Approximation  given  by  formula  (A-79)  to  the  integral  A 

Coefficients  appearing  in  formula  (C-46) 

Coefficients  appearing  in  formula  (22) 

Generalised  coordinate  for  mode  number  k (see  formula  (6)) 

Quantity  defining  amount  of  harmonic  constituent  in  bkit)  (see  formula  (7)) 

Complex  conjugate  of  bk{<o) 

The  integral  (A-76)  ■ 

1 

Approximation  given  by  formula  (A-83)  to  the  integral  B 
Coefficients  appearing  in  formula  (C-47) 

Coefficients  appearing  in  formula  (30) 

Chord  length  of  planform  S at  spanwise  position  y (see  Fig.  1)  ^ 

Coefficients  appearing  in  formulae  (C-44)  and  (C-45)  j 

Matrix  appearing  in  formula  (55)  and  defined  immediately  afterwards  j 

Coefficients  defined  in  formulae  (52)  and  (53)  j 

Defined  by  formula  (A-38) 

Diagonal  matrix  with  (80)  as  general  diagonal  element 
Coefficients  appearing  in  formula  (C-44) 

Diagonal  matrix  with  (79)  as  general  diagonal  element 
The  functions  defined  by  formula  (C-62) 

Coefficients  appearing  in  formula  (C-45) 

The  function  defined  by  equation  (C-23) 

Defined  by  formula  (A-73) 

Matrix  appearing  in  equation  (93)  and  defined  immediately  before 
The  function  defined  by  equation  (C-30) 

The  functions  defined  by  equations  (C-34)  and  (C-35) 

Coefficients  appearing  in  equation  (C-75) 

Row  matrix  with  (87)  as  general  element 
Quantities  defined  by  formula  (118) 

The  function  defined  by  equation  (C-24) 

Set  of  interpolation  polynomials  defined  by  equations  (28) 

Set  of  interpolation  polynomials  defined  by  equations  (A-27) 

Set  of  interpolation  polynomials  defined  by  equations  (A-62) 

The  function  defined  by  equation  (C-31) 
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Gn(a)  The  functions  defined  by  equations  (C-36) 

Gr(A ) Coefficients  appearing  in  equation  (C-75) 

Gp"*'  Quantities  defined  in  formulae  (72) 

G/**'  Quantities  defined  in  formulae  (103) 

G/'*  Quantities  defined  in  formulae  (A-32) 

op  Quantities  defined  in  formulae  (A-67) 

h(a)  The  function  defined  in  formula  (C-67) 

/i«(a)  The  function  defined  by  equation  (C-70) 

hi''\€o)  Set  of  interpolation  polynomials  defined  by  equations  (24) 

hP(()  Set  of  interpolation  polynomials  defined  by  equations  (A-8) 

[//]  Diagonal  matrix  with  (94)  as  general  diagonal  element 
H(a)  The  function  defined  by  equation  (C-65) 

Quantities  defined  in  formulae  (71) 

Quantities  defined  in  formulae  (102) 

//j'*  Quantities  defined  in  formulae  (A- 13) 

/ Integral  introduced  in  formula  (A- 14) 

7^'*  Approximation  to  / 

/i(a)  Modified  Bessel  function  of  the  first  kind  and  first  order 

7t"’(^,  V.  Vo',  V,  Moo)  Quantities  defined  in  formula  (1 14) 

V,  Vo',  t>.  Moo)  Quantities  defined  in  formula  (117) 

/(a)  The  function  defined  in  equation  (C-22) 

J Integral  introduced  in  formula  (A-33) 

7^'*  Approximation  to  / 

Approximation  to  J 
Approximation  to  J 
Approximation  to  J 

y;  p,  Mooj  Quantity  defined  in  formula  (C-9) 

k(a)  The  function  defined  in  formula  (C-67) 
k„(a)  The  function  defined  by  equation  (C-71) 

kP"\vo)  The  function  defined  by  equation  (86) 

K Integral  introduced  in  formula  (A-68) 

K(a)  The  function  defined  by  equation  (C-66) 

A^i(«)  Modified  Bessel  function  of  the  second  kind  and  first  order 
X y \ 

y;  V,  Ma> I Subsonic  kernel  function  defined  by  equation  (17) 

/ Typical  length  of  the  planform  S 
lk(x,  y;  p.  Moo)  Loading  function  introduced  in  expression  (14) 

fk(xo,  yo)  Approximation  (22)  to  the  loading  function  lk(xo,  yo',  p,  M) 
[/It]  Column  matrix  with  (92)  as  general  element 
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[^k(xo.  yo)]  Matrix  consisting  of  the  single  element  fk(xo,  yo) 

lr(4)  Orthogonal  polynomials  of  degree  r in  ^ satisfying  equation  (A-1) 

[Z.]  Matrix  with  (1 1 1)  as  general  element 

/.i(a)  Modified  Struve  function  which  is  related  to  the  Struve  function  |//|.i((a)  by 
means  of  formula  (C-18).  The  function  is  defined  by  equation  (C-21) 

Lpj  Quantities  defined  in  formulae  (108)  and  (109) 

m Number  of  spanwise  loading  functions  (see  formula  (22)) 

ntH  Integer  defined  by  formula  (48) 

Moo  Mach  number  defined  by  formula  (IS) 

M Integer  defined  by  formula  (155) 

M Number  of  integration  points  spanwise  for  evaluation  of  (A>.p:r.i  from  the  numerical 

formula  (99) 

Mh  Integer  defined  by  formula  (106) 

M/  Set  of  N positive  integers  defined  by  formulae  (166) 

M,{A)  Coefficients  appearing  in  equation  (C-74) 

n Number  of  chordwise  loading  functions  (see  formula  (22)) 

/V  Number  of  integration  points  chordwise  for  evaluation  of  il't.p.r.s  from  the  numeri- 

cal formula  (99) 

Nr(A ) Coefficients  appearing  in  equation  (C-74) 

P(a)  Function  defined  in  formula  (C-15) 

P\  Quantities  defined  in  equations  (122),  (124),  (126)  or  (128) 

q A positive  integer  for  determining  M from  formula  (155) 

qi  A set  of  N positive  integers  introduced  in  formula  (165) 

<?,*(»<’,  Moo)  Generalised  airforce  coefficients,  defined  by  formula  (19) 

Oik  Approximation  (34)  to  the  generalised  airforce  coefficient  0,t(i',  Moo) 

[(}]  Matrix  appearing  in  formula  (84)  and  defined  immediately  before 
0\(v)  Function  defined  by  formula  (B-46) 

R Quantity  defined  by  formula  (18) 

R\(rj)  Function  defined  by  formula  (B-39) 

s Semi-span  on  the  planform  S 
S Wing  planform 

S(a)  Function  defined  by  formula  (C-12) 

Sa(i7)  Function  defined  by  formula  (B-22) 

/ Time 

Tr(rj)  Qrthogonal  polynomials  of  degree  r in  t;  satisfying  equation  (A-56).  Also 
Chebyshev  polynomial  defined  in  formula  (C-43) 

[ U]  Matrix  with  (1 1 2)  as  general  element 

C4.5(jt,  y;  M)  Quantity  defined  by  formula  (32) 

V Speed  of  main  flow.  It  is  in  the  positive  x direction 

W(x,  y,  I)  Quantity  defined  by  equation  (5) 
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Quantity  defined  by  equation  (1) 

Quantity  defined  by  equation  (2) 

Rectangular  cartesian  coordinates  of  a point  relative  to  a frame  fixed  with  respect 
to  the  mean  position  of  the  wing 

X coordinate  of  the  leading  edge  of  the  planform  S at  spanwise  position  y (see  Fig. 

1) 

Quantities  defined  by  formulae  (69) 

Quantities  defined  by  formulae  (70) 

Quantities  defined  by  formulae  (104) 

Quantity  defined  by  formula  (C-10) 

Quantities  defined  by  formulae  (68) 

Quantities  defined  by  formulae  (105) 

Quantity  obtained  from  equations  (3)  and  (4) 

Quantity  obtained  from  equations  (3)  and  (4) 

Displacement  from  the  plane  of  S of  a point  on  the  wing  top  surface 
Displacement  from  the  plane  of  5 of  a point  on  the  wing  bottom  surface 
Reduced  upwash  function,  defined  by  formula  (11) 

Approximation  (3 1 ) to  the  upwash  function  o/t  (x,  y ) corresponding  to  the  approx- 
imation tk(xo,  yo)  to  the  loading  function  /jt(xo,  yo;  i',  A/) 

Matrix  appearing  in  formula  (84)  and  defined  immediately  before 

Quantities  defined  by  formulae  (A-80) 

Quantities  defined  by  formulae  (A-84) 

Euler’s  constant 

Orthogonal  polynomials  of  degree  r in  17  satisfying  equation  (A-21) 
Kronecker’s  delta,  defined  in  formula  (26) 

Function  defined  by  equation  (C-32) 

Quantity  defined  by  formula  (208) 

Function  defined  by  equation  (C-33) 

Quantities  defined  by  formula  (207) 

Modal  function  for  the  mode  number  k as  it  appears  in  equation  (6) 

Set  of  A points  defined  in  equations  (123),  (125),  (127)  or  (129) 

Parametric  coordinate  defined  by  formula  (20) 

Parametric  coordinate  defined  by  formula  (20) 

Set  of  m points  in  (-1, 1)  defined  by  formulae  (27) 

Set  of  M points  in  (-1, 1)  defined  by  formulae  (101) 

Set  of  r points  in  (-1, 1)  defined  by  formulae  (A-25) 

Matrix  appearing  in  formula  (55)  and  defined  immediately  afterwards 
Quantities  defined  by  formulae  (A-7) 

Quantities  defined  by  formulae  (36) 
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_ I +1  For  symmetric  modes  of  oscillation 
~ I - 1 For  anti-symmetric  modes  of  oscillation 

A Number  of  integration  points  used  in  the  numerical  evaluation  of  the  integral 

(121) 

finia)  Function  defined  by  equation  (C-68) 

p Frequency  parameter  defined  by  formula  (12) 

Pni<x)  Function  defined  by  equation  (C-69) 

4 Parametric  coordinate  defined  by  formula  (20) 

Parametric  coordinate  defined  by  formula  (20) 

Set  of  n points  in  (0, 1)  defined  by  formula  (23) 

4^^  Set  of  N points  in  (0, 1)  defined  by  formula  ((100) 

1/'^'  Set  of  N points  in  (0, 1)  defined  by  formula  (168) 

Set  of  r points  in  (0, 1)  defined  by  formula  (A-6) 
p Density  of  the  fluid  in  the  main  flow 
Variable  defined  in  equation  (B-1) 

Quantities  defined  by  formulae  (A-26) 

Quantities  defined  by  formulae  (A-61) 

Lt,]  A row  matrix  appearing  in  formula  (63)  and  defined  immediately  afterwards 
xV  Quantities  defined  by  formulae  (A-60) 

Xi\r.s  Quantities  defined  by  formulae  (60) 

iAi.p;r.j  Quantities  defined  by  formulae  (37) 

'Pi.p;r.j  Quantities  defined  by  formulae  (54) 

[’P]  Matrix  appearing  in  formula  (55)  and  defined  immediately  afterwards 
0)  Circular  frequency  of  harmonic  oscillation 

' Indicates  differentiation 

" Indicates  differentiation  twice 


38 


REFERENCES 


No.  Author(s)  Title,  etc. 

1 H.  Multhop  Methods  for  calculating  the  lift  distribution  of  wings.  (Sub- 

sonic lifting-surface  theory.) 

A.R.C.  R.  & M.  No.  2884  (1950) 

2 H.C.  Garner  Multhopp’s  subsonic  lifting-surface  theory  of  wings  in  slow 

pitching  oscillations. 

A.R.C.  R.  & M.  No.  2885  (1952) 


3 W.  E.  A.  Acum 


4 J.  R.  Richardson 


5 D.  E.  Davies 


Theory  of  lifting  surfaces  oscillating  at  general  frequencies  in  a 
subsonic  stream. 

A.R.C.  R.  & M.  No.  3557  (1959) 


A method  for  calculating  the  lifting  forces  on  wings  (unsteady 
subsonic  and  supersonic  lifting  surface  theory). 

A.R.C.  R.  & M.  No.  3157  (1955) 


Calculation  of  unsteady  generalized  airforces  on  a thin  wing 
oscillating  harmonically  in  subsonic  flow. 

A.R.C.  R.  & M.  No.  3049  (1963) 


6 H.  C.  Garner  and  D.  A.  Fox Algol  60  programme  for  Multhopp’s  low-frequency  subsonic 

lifting-surface  theory. 

A.R.C.  R.  & M.  No.  3517  (1966) 


7 G.  Long 


An  improved  method  for  calculating  generalized  airforces  on 
oscillating  wings  in  subsonic  flow. 

A.R.C.  R.  & M.  No.  3657  (1969) 


8 P.  J.  Zandbergen,  T.  E.  Labrujere 
and  J.  G.  Wouters 


A new  approach  to  the  numerical  solution  of  the  equation  of 
subsonic  lifting  surface  theory. 

NLR  Report  TR  G.49  (1967) 


9 Doris  E.  Lehrian  and  H.  C.  Garner  . . Theoretical  calculation  of  generalised  forces  and  load  distribu- 

tion on  wings  oscillating  at  general  frequency  in  a subsonic 
stream. 

A.R.C.  R.  & M.  No.  3710  (1971) 

10  B.  L.  Hewitt  and  W.Kellaway A new  treatment  of  the  subsonic  lifting  surface  problem  using 

curvilinear  coordinates.  Part  I — Details  of  the  method  as 
applied  to  regular  surfaces  with  finite  tip  chords  and  a prelimi- 
nary set  of  numerical  results. 

BAC  Preston  Division  Report  Ae.290  (S.  & T.  Memorandum 
20/68)  (1968) 

An  application  of  Flax’s  variational  principle  to  lifting-surface 
theory. 

A.R.C.  R.  & M.  No.  3564  (1967) 


11  D.  E.  Davies 


12  D.  E.  Williams Three-dimensional  subsonic  theory. 

AGARD  Manual  on  Aeroelasticity,  Vol.  II,  Chapter  3 (1961) 


12  D.E.  Williams Three-dimensional  subsonic  theory. 

AGARD  Manual  on  Aeroelasticity,  Vol.  II,  Chapter  3 (1961) 

13  P.T.  Hsu  Some  recent  developments  in  the  flutter  analysis  of  low  aspect 

ratio  wings. 

Proceedings  of  the  National  Specialists’  Meeting  on  Dynamics 
and  Aeroelasticity,  Fort  Worth,  Texas,  6-7  November  1958 
(Institute  of  the  Aeronautical  Sciences) 


14  E.  C.  Titchmarsh The  theory  of  functions. 

Oxford  University  Press 

15  N.  W.  McLachlan  Bessel  functions  for  engineers. 

Oxford  University  Press  (1934,  reprinted  1948) 

16  C.  W.  Clenshaw The  numerical  solution  of  linear  differential  equations  in 

Chebyshev  series. 

Proc.  Camb.  Phil.  Soc.  53,  pp.  134-149  (1957) 

17  D.  E.  Davies  On  the  use  of  Fortran  programs  for  evaluating  the  generalised 

airforces  and  aerodynamic  loading  on  a flat  plate  wing  oscillat- 
ing harmonically  in  subsonic  flow. 


RAE  Technical  Memorandum  Structures  881  (1976) 

t 


40 


T 


APPENDIX  A 


Integration  Formulae 


Integrand  with  Weight  Function  '/(T— 

Let  /r(f ),  r = 0, 1, 2, . . . , be  the  set  of  polynomials  of  degree  /•  in  f which  satisfy  the  orthogonality  relations 


C mm  Sr, 

h ^ c 


(A-l) 


where  is  Kronecker’s  delta 


Srs  = 


1 r = s 


0 r^s. 


(A-2) 


If  we  make  the  change  of  variables 


^ = 2(1 -cos  6) 


(A-3) 


in  the  integrand  of  (A-l)  we  get 


(A-4) 


from  which  it  follows  that 


/I -cos  ff\  [2  cos (r  + 2)0 
\ 2 IT  cos  2O 


(A-5) 


for  this  is  a polynomial  of  degree  r in  the  original  variable  ( and  both  (A-l)  and  (A-4)  are  satisfied.  Since 
is  positive  for  ^ in  (0, 1 ) the  polynomial  I,  (^ ) will  have  r zeros  in  (0, 1 ).  Let  these  zeros  be  denoted  by 
/ = 1, 2, . . . , r.The  locations  of  these  zeros  are  obtained  directly  from  (A-5)  and  (A-3)  and  are  given  by 


^l'>  = i(l- cos  »!'>),  / = l,2,...,r, 


(A-6) 


where 


, V / = 1,  2, . . . , r. 


2r-l-l 


(A-7) 


We  define  the  interpolation  polynomials  by  means  of  the  formulae 


ri  (4r%).  « = l,2,...,r. 

k-1  / 


(r) 


(A-8) 


These  polynomials  are  of  degree  (r  — 1)  in  f and  have  the  property 

They  are  given  also  by  the  alternative  formula 


(A-9) 


m 
/;(^h(^-fh 


j = 1, 2, . . . , r. 


(A-10) 


I 

J 
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we  get  a corresponding  approximation  to  I which  is  given  by 

/<o=  £ (A-16) 

1*1 

by  replacing  f{()  on  the  right  of  (A- 14)  by  Formula  (A-16)  is  the  Gaussian  numerical  integration 

formula  for  a weight  function  when  r integration  points  are  used.  Because  of  the  orthogonality 

relationship  (A- 1 ) we  have  that  the  approximation  is  exactly  equal  to  / whenever  /(f)  is  a polynomial  in  f of 
degree  «2r  — 1 . If  /(f)  is  not  such  a polynomial  then  we  can  give  an  estimate  for  the  error  / - 1^'^  by  using  the 
following  procedure. 


We  can  show  that  the  approximation  to  I given  by  formula  (A- 16)  is  also  given  by  the  formula 


I 


I 


/<'>  = ^ lim 

2p  p-*® 


f'  / 1 -cos  g\{l  + (-1)*’'^*  cos  (2r  + l)pg} 
1 ' 2 / {cos(r+l)tfp 


(cosifl)^  de. 


(A-17) 


where  p runs  through  the  positive  integers.  To  show  that  (A-16)  and  (A-17)  are  equivalent  formulae  for 
continuous  functions /(f),  the  range  (0,  ir)  of  integration  in  (A-17)  is  divided  into  r subranges  with  each  of  the 
subranges  containing  one  of  the  zeros  6 = eY^,  i = 1, 2, . . . , r of  cos  (r  + ^)0  internal  to  it.  For  each  of  the 
subranges  a new  variable  of  integration  is  introduced  by  translation  of  0 so  that  the  zero  of  the  new  variable  is 
at  = dY^  in  the  ith  subrange  and  Fejfr’s  integral  (see  Ref.  14)  is  thereby  obtained  for  that  subrange.  Then,  by 
applying  Fejer’s  theorem  (see  Ref.  14)  to  each  of  the  subranges  the  result  that  (A-16)  and  (A-17)  are 
equivalent  formulae  is  obtained. 

From  (A-17)  we  get 


/<'>=hjn|  /(^— f (-!)*(  1-^)  cos  (2r+l)yflj(cosifl)^ 

- i ' -f) 


de 

X 


x[cos  {(2r  + l)j  - 2}^  + cos  {(2r  + l)s  + ] cos  ko  d0 

ill  '-»■(■-?)  i'«)x 

X[/(2r+l)l-l(f)+/(2r+l)j(f)]^-^<ff 
= /+lim  Y,  (~l)*(l ^)[a(2r+l).-l+fl(2r+l)j] 

P-.00  V 2 1-1  V p/ 


where 


The  formula 


/<'’  = /+ lim 

p-*oo 


(— l)*[a(2,+  l),-l  +a(2r+l)i] 


(A-18) 


(A-19) 


(A-20) 


is  equivalent  to  the  formula  (A-18)  provided  that  the  limit  on  the  right  hand  side  of  (A-20)  exists.  The  formula 
(A-18)  is  valid  whether  the  limit  on  the  right  hand  side  of  (A-20)  exists  or  not.  Taking  the  limit  as  p -»ao  in 
(A-20)  corresponds  to  summing  an  infinite  series  directly  whereas  taking  the  limit  as  p-*<X)  in  (A-18) 
corresponds  to  summing  the  same  series  by  arithmetic  means  or  by  Cesikro’s  means  of  first  order. 

The  error  7-/^''  may  now  be  gauged  from  either  formula  (A- 1 8)  or  formula  (A-20).  The  coefficients  a„  tend 
to  zero  as  n tends  to  infinity,  as  is  shown  by  an  application  of  the  Riemann-Lebesgue  theorem  (see  Ref.  14)  to 
formula  (A-19)  after  changing  the  integration  variable  from  f to  0.  The  error  converges  to  zero  as  r 
increases  indefinitely  because  the  multipliers  of  formula  (A-16)  are,  according  to  formula  (A-13),  all 
positive. 

The  formula  (A- 1 8)  is  also  valid  for  more  general  /(f)  than  continuous  functions,  e.g.  for  /(f)  having  a finite 
number  of  jump  discontinuities  for  f in  (0, 1),  but  in  this  paper  we  are  interested  only  in  continuous  functions 

/(f). 
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Intesnuid  with  Weight  Fiuctioa 

Let  y,(ri),  r = 0, 1, 2, ...  be  the  set  of  polynomials  of  degree  r in  17  which  satisfy  the  orthogonality  relations 


J ^ yr(7j)y,(l7)N/l-7J^d7}  =«„. 

(A-21) 

If  we  make  the  change  of  variables 

17  =cos 

(A-22) 

in  the  integrand  of  (A-20)  we  get 

yr(cm  <t>f  d<l>  = S„ 

Jo 

(A-23) 

from  which  it  follows  that 

, ,,  f2%in(r  + l)4> 

yr(cos4)-\J^  . , . 

’ TT  Sin  0 

(A-24) 

for  this  is  a polynomial  of  degree  r in  the  original  variable  T7  and  both  (A-21)  and  (A-23)  are  satisfied.  Since 
is  positive  for  rj  in  (-1, 1)  the  polynomial  '>',(77 ) will  have  r zeros  in  (-1, 1).  Let  these  zeros  be  denoted 
by  rf]'\  J = 1,2, . . . , r.  The  locations  of  these  zeros  are  obtained  directly  from  (A-24)  and  (A-22)  and  are  given 
by 

77}'’  = cos  j = \,2,...,r. 

(A-25) 

where 


Then,  from  formula  (A-29)  we  get 


sin  {r  + \)<f> 
sin  <ii(cos  - cos 


(A-31) 


Therefore 


= f «r(cos  <^Xsin  <l>f  d<f> 

JO 


(-ly*.  ■ ^(r)^2  f"(cosr<^-cos(r  + 2)<^) 
2(r+l)^® 


771 

(cos  <t>  - COS  4>  ) 0 


2(r+l) 


Ml/  vf  ^ 


-^j|y(sin.>i'ycos(r  + l)^!" 

-(TfliW"*!")" 


We  are  interested  in  the  numerical  evaluation  of  the  integral 

•^=1  ^ dll 


(A-32) 


(A-33) 


where  /(17 ) is  an  arbitrary  continuous  function  of  tj  in  (- 1 , 1 ).  If  we  take  an  approximation  f^'Xrt ) to  /(t; ) which 
is  the  interpolation  polynomial  of  degree  (r  — 1)  in  tj 


f'\ri)=  I nvWXri), 


i-i 


(A-34) 


we  get  a corresponding  approximation  to  J which  is  given  by 


/'>=  I 

/-I 


(A-35) 


by  replacing  /(ij)  on  the  right  of  (A-33)  by  /'’(tj). 

Formula  (A-35)  is  the  Gaussian  numerical  integration  formula  for  a weight  function  when  r 

integration  points  are  used.  Because  of  the  orthogonality  relationship  (A-21)  we  have  that  the  approximation 
is  exactly  equal  to  / whenever  /(tj)  is  a polynomial  in  17  of  degree  «2r  - 1.  If  /(tj)  is  not  such  a polynomial 
then  we  can  give  an  estimate  for  the  error 
Instead  of  formula  (A-35)  we  can  write 


= lim  - j /(cos  (sin  4>?  d<t> 

p-®pJo  \sin(r  + l)<p/ 


(A-36) 


where  p runs  through  the  positive  integers.  That  (A-35)  and  (A-36)  are  equivalent  formulae  for  continuous 
functions  /(tj  ) is  demonstrated  if  we  divide  the  range  (0,  tt)  of  integration  in  (A-36)  into  r subranges  with  each 
of  the  subranges  containing  one  of  the  zeros  <l>  = <^S'\  / = 1, 2, . . . , r,  of  sin  (r+ 1)<^  internal  to  it,  and  then 
apply  Fejdr’s  integral  formula  (see  Ref.  14)  to  each  of  the  subranges. 
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From  (A-36)  we  get 


y^'*=  lim  f /(cos<^)|l  + 2 Y,  ( 1 ) cos  2(r+ l)i<^|(sin 
P-®  Jo  • »-i  ' P'  ^ 

= f /(cos  <^Xsin  0)^  + iim  £ (l~“)  f /(cos<^)x 

Jo  p/ Jo 

x(sin  {2(r  + 1 )j  + 1 }^  - sin  {2(r  + 1 )s  - 1 }^]  sin  0 

= j ^ /(7j)Vl-TJ^  dTt  + 

^i(^~p)  1 1 y2(r+i),-2]'/ 1 - dv 

= /+lim  \j—  £ (l — V</2(f+i)»  ~<^2(r+i)i-2]  (A-37) 

p-ao  » 2 1.1  ' p/ 

where 

</m=|  f(v)ym(v)^^-v^  dij.  (A-38) 

We  can  write,  instead  of  (A-37),  the  formula 

J^'^  = J+^~\im  Y [<^2(f+i)i ~d2(r+i)i-2)  (A-39) 

’ 2 p-»®  1.1 

provided  that  the  limit  on  the  right  hand  side  of  (A-39)  exists.  The  formula  (A-37)  is  valid  whether  the  right 
hand  side  of  (A-39)  exists  or  not. 

The  error  may  now  be  gauged  from  either  formula  (A-37)  or  (A-39).  The  coefficients  d„  tend  to  zero 

as  m tends  to  infinity,  as  is  shown  by  an  application  of  the  Riemann-Lebesgue  theorem  (see  Ref.  14)  to 
formula  (A-38)  after  changing  the  integration  variable  from  77  to  <(>.  The  error  converges  to  zero  as  r 

increases  indefinitely  because  the  multipliers  G)'^  of  formula  (A-35)  are,  according  to  formula  (A-32),  all 
positive. 

The  formula  (A-37)  is  also  valid  for  more  general  /(tj),  e.g.  /(tj)  having  a finite  number  of  jump 
discontinuities  for  t)  in  (-1, 1),  but  in  this  paper  we  are  interested  only  in  continuous  functions  /(rj). 

We  are  also  interested  in  the  analytical  evaluation  of  the  integral 


nnWER 


(A-40) 


Let  us  write 


M V f"'  yAvy^^-v^  . 
J-1  (v-C)  ’’  y'rint^  J-i  ^ 


(A-41) 


Now,  we  have 


' (cos  - cos  (r  + 2)<^) 


J-1  •^2ir  Jo  (cos  -cos  i/r) 

— ('■■*■2)^) 


d<t> 


2 sin  0 

= -V^cos  (r+ 1)0 


(A-42) 


where 


^ = cos  til. 


On  changing  f to  tj}  * in  (A-42),  and  consequently,  on  changing  tji  to  we  get  immediately 


^4-1  ^r\^) 

J j 1 - 77^  drj  = -yf^  cos  (r  + 

If  we  substitute  from  (A-42)  and  (A-44)  into  (A-41)  and  use  (A-30)  we  get 


1: 


*' gi'\vy^i-v^ . 7r(-iy 


(v-£)  (r  + l) 


r,  ...(,K2,(cos(r-(-l)<^-cos(r+l)<^;'>) 

(T^ • 


Hence,  if 


£=vV  k^i 


we  get  from  (A-45) 


and  if 


we  get  from  (A-45) 


1-, 


:<o (rKif-C-ty  ]. 


(r+i)(vr-vi'r 


r - Y,*'' 
i~Vi 


. ^(-lyr.  / (rVif  r 

rll -(»?/)  Ihj7  COS  (r+l ),;- 

) Ldt  Jf-V?' 


f * 

I s;  VI/”  * 'I  , '•  V */  I 

J-,  ‘'"'■(TMji 

_^^)/r^  , (,),2.rd/d(^cos(r  + l)<^1 

) ^ L (d  cos  ill)/ dl/l 


(^+1)' 

= 0. 


Also  from  (A-45),  on  differentiating  with  respect  to  £,  we  get 

i: 


iv-£f  d^J-i  (rt-tT) 
>r(-iy"' 


(»'+l) 


fY  /_(>K2i(cos(r+l)<^-cos(r+l)<^yV 
U (’ll)!  + 


Hence,  if 


we  get  from  (A-50) 


£=vV  k^i 


f"gr(r>y/T^ 

i.  ^ (r+ix»jr-»fyy^‘  ^ ^ ’’ 
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(A-43) 


(A-44) 


(A-45) 


(A-46) 


(A-47) 


(A-48) 


(A-49) 


i 

1 

1 


I 

I 


(A-50) 


(A-51) 


(A-52) 


i 

t 

i 


■f 


whereas  if 


-Vi 


(A-53) 


we  get  from  (A-50) 


r'gnv)^T^. 

i,  iv-v’fr"" 


7r(-iy,.  , „^,,r  \ d I d 


ii'-'-wiy 


(r+lJ  ,, 


2(r  + l)‘‘  ' 'Lsin  <fiA  sin  lA  rf<A 

iv*‘ri_r  '-Vtfr.  i , ^ (>•  + \)*^  sin(r  + l).A1 
2^  y [ (Vi  ) 1[(  (sin il>f  (sini^y 


= -f(r+l). 


(A-54) 


Collecting  together  the  results  (A-52)  and  (A-54)  into  one  formula  we  get 


I 


-f(r  + l) 

] Q / ly**! 


fc=y 

k^i 


(A-55) 


Integrand  with  Weight  Function  l/'/T— 7 

Let  T,(ti),  r = 0, 1, 2, ....  be  the  set  of  polynomials  of  degree  r in  17  which  satisfy  the  orthogonality  relations 


r Tr{v)Uv)-r==^=S„ 
J-i  vl-17 


(A-56) 


If  we  make  the  change  of  variables 


r)  = cosd> 


(A-57) 


we  get 


from  which  it  follows  that 


f T,(<x>s<f>)T,(cos4>)d<f>  = S„ 
Jo 


(A-58) 


r,(cos<^)  = 


f_L 


r = 0 


cos  r<i>  r^O 


(A-59) 


for  this  is  a polynomial  of  degree  r in  the  original  variable  17  and  both  (A-56)  and  (A-58)  are  satisfied.  Since 
l/Vl-ry*  is  positive  for  rj  in  (-1, 1)  the  polynomial  7,(77)  will  have  r zeros  in  (-1, 1).  Let  these  zeros  be 
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r 

f 

1 

\ 

denoted  by  / = 1. 2, . . 

are  given  by 

. , r.  The  locations  of  these  zeros  are  obtained  directly  from  (A-59)  and  (A-57)  and 

where 


2/-1 

<^,  - 2r  ^ 


We  define  interpolation  polynomials  gj'^v)  by  means  of  the  formulae 

n (JiirS)  /=1’2 r. 

t-1  'ATj  AT*  > 

k^i 

These  polynomials  are  of  degree  (r- 1)  in  tj  and  have  the  property 

They  are  also  given  by  the  alternative  formula 

.(r)/  X _ Tr{r\) 

where  the  dash  represents  differentiation  with  respect  to  the  argument  of  the  appropriate  function.  By 
differentiating  formula  (A'59)  with  respect  to  <t>  we  easily  establish  that 

’ rr  sm  q>j 

Then,  from  formula  (A-64)  we  get 

g/'’(^)  = 


(A-60) 


(A-61) 


(A-62) 


(A-63) 


(A-64) 


(0- 


COSr,^ 


(A-65) 


- sin  <A/ 


(cos  cos 


(A-66) 


Therefore 


= [ g]'\cos  <l>)  d<l> 

Jo 

(“ly""’  • 7(rlf’'  COS'’<^ 

= sm<^5'  7 

r Jo  (cos<^-cos<^;  0 


■ sin  r<^} 


7r(-iy^ 

r 

_ TT 

r * 

We  are  interested  in  the  numerical  evaluation  of  the  integral 

dr, 


(A-67) 


(A-68) 
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where  f{r\ ) is  an  arbitrary  continuous  function  of  t/  in  (—  1 , 1 ).  If  we  take  an  approximation  to  f(r\ ) which 

is  the  interpolation  polynomial  of  degree  (r- 1)  in  tj 


nv)=  i (A-69) 

/“I 

we  get  a corresponding  approximation  to  K which  is  given  by 

^ GfyCrD  (A-70) 

y-i 


by  replacing  /(tj)  on  the  right  of  fA-68)  by  f'X-n).  Formula  (A-70)  is  the  Gaussian  numerical  integration 
formula  for  a weight  function  when  r integration  points  are  used.  Because  of  the  orthogonality 

relationship  (A-56)  we  have  that  the  approximation  is  exactly  equal  to  K whenever /(17 ) is  a polynomial  in 
r\  of  degree  «2r-l.  If  /(rj)  is  not  such  a polynomial  then  we  can  give  an  estimate  for  the  error  K-K^'^. 

Instead  of  formula  (A-70)  we  can  write 


r(r) 


1 r’' 

= lim—  f(cos<f>) 

p-oo  2p  Jo 


{1  + (-!)'’  cos  2rp<f>} 
(cos  r<^)^ 


d4> 


(A-71) 


where  p runs  through  the  positive  integers.  That  (A-70)  and  (A-71)  are  equivalent  formulae  for  continuous 
functions  /(rj ) is  demonstrated  if  we  divide  the  range  (0,  n)  of  integration  in  (A-71)  into  r subranges  with  each 
of  the  subranges  containing  one  of  the  zeros  <t>  = / = 1, 2, . . . , r,  of  cos  r<t>  internal  to  it,  and  then  apply 

Fejer’s  integral  formula  (see  Ref.  14)  to  each  of  the  subranges. 

From  (A-71)  we  get 


f /(cos<^)jl-l-2  f (-!)*(  1--)  cos  2r5<^ly/<^ 

P-^Jo  I i-l  V pt  ) 

= r lim  I (-irfi-;;)  Kv)T2r,{v)-j=^ 

■'-1  vl— 77  P-"®  ' P'  J-1  vl  — T) 

= K+\xm^i  (-1)‘(  1 (A-72) 

P-OO  \ p! 

where 

(A-73) 

J-1  vl  — TJ 

We  can  write,  i 'stead  of  (A-72),  the  formula 

p 

A'‘'’  = Ar  + v^lim  S (-l)y2„  (A-74) 

p-*”°  j-i 

provided  that  the  limit  on  the  right  hand  side  of  (A-74)  exists.  The  formula  (A-72)  is  valid  whether  the  right 
hand  side  of  (A-74)  exists  or  not. 

The  error  K - may  now  be  gauged  from  either  formula  (A-72)  or  (A-74).  The  coefficients  f„  tend  to 
zero  as  m tends  to  infinity,  as  is  shown  by  an  application  of  the  Riemann-Lebesgue  theorem  (see  Ref.  14)  to 
formula  (A-73)  after  changing  the  integration  variable  from  tj  to  tj).  The  error  K - converges  to  zero  as  r 
increases  indefinitely  because  the  multipliers  of  formula  (A-70)  are,  according  to  formula  (A-67),  all 
positive. 

The  formula  (A-72)  is  also  valid  for  more  general  /(tj),  e.g.  /(tj)  having  a finite  number  of  jump 
discontinuities  for  tj  in  (-1, 1),  but  in  this  paper  we  are  interested  only  in  continuous  functions  /(tj). 
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Intetruds  with  Weight  Functions  V(1  + »})/(! -i})  and  n/(1-i|)/(1  + i|) 
We  obtain  numerical  formulae  of  integration  for  the  integrals 


■1, 


dri 


and 


B 


~L, 


by  application  of  the  numerical  formula  (A- 16).  If  we  make  the  transformation  of  variables 

i=ki-v) 

in  the  integrand  in  (A-75)  we  get 


(A-75) 


(A-76) 


(A-77) 


(A-78) 


Then,  on  applying  the  numerical  integration  formula  (A-16),  for  weight  function  V(l— ^)/fwith  r integration 
points,  to  (A-78)  we  get  the  approximation  to  A which  is  given  by 


A<'>  = - 


In 


(2r  + l)/tr, 


I (l  + ar)f{ar) 


(A-79) 


where 


af^=l-2^f 


If  we  make  the  transformation  of  variables 


^ = 2(l  + r7) 


in  the  integrand  in  (A-76)  we  get 


B=2 

Jo  ’ S 


(A-80) 


(A-81) 


(A-82) 


Then,  on  applying  the  numerical  integration  formula  (A-16),  for  weight  function  >/(l-f)/fwith  r integration 
points,  to  (A-82)  we  get  the  approximation  B*'*  to  B which  is  given  by 


5<')  = - 


27r 


(2/-+1); 


I (i-^ro/osn 


(A-83) 


where 


^r>=2f</-v,-i 

= cos(^7r),  /=1,2 r. 


(A-84) 


The  approximation  A is  exactly  equal  to  A and  the  approximation  B*''  is  exactly  equal  to  B whenever  f(r) ) 
is  a polynomial  in  r\  of  degree  «2r  - 1.  If  /(tj)  is  not  such  a polynomial  then  we  can  give  estimates  for  the  errors 
A - A ''’and  B-B*'’,  which  are  based  on  formula  (A- 18). 


51 


Spanwise  Integntion  Formulae 

Finally,  in  this  Appendix,  let  us  consider  several  different  evaluations  of  the  integral 

•^=1  dv  (A-85) 

which  are  based  on  the  numerical  integration  formulae  (A-35),  (A-70),  (A-79)  and  (A-83).  These  integration 
formulae  are  required  for  justifying  formula  (121)  of  the  main  text. 

Straightforwardly  we  have  the  numerical  estimate  from  formula  (A-35), 

= 7^  i /(T?r)[l  - (A-86) 

If  we  write 

^ fiv)0^  + v)\f^Z^d7i  (A-87) 

and  apply  the  formula  of  integration  (A-83)  we  get  the  estimate  for  J which  is  given  by 

•^'■’'  = rTrx  i /(T?g>\)[l-(rj^%\)^l  (A-88) 

since 

/=1,2 r.  (A-89) 

If  we  write 

/=f"'/(^)(l-’?)V^^’7  (A-90) 

J-l  ''1-7/ 

and  apply  the  formula  of  integration  (A-79)  we  get  the  estimate  for  / which  is  given  by 

= i (A-91) 

V ’^2)1-1 

since 

ar=T/Sl\  /=l,2,...,r.  (A-92) 

If  we  write 


/=fV(T7Xl-T/^)^j^  (A-93) 

and  apply  the  formula  of  integration  (A-70)  we  get  the  estimate  for  J which  is  given  by 

i /(T/i;:l)[i  -(7/‘2%VV1  (a-94) 

r /-I 

since 

; = 1.2 r (A-95) 

Formula  (121)  of  the  main  text  is  obtained  by  applying  one  of  the  estimates  (A-86),  (A-88),  (A-91)  and 
(A-94)  for  the  integral  appearing  on  its  left  hand  side. 


iJL 
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APPENDIX  B 


Derivation  of  Some  Identities  used  in  the  Main  Text 

In  this  Appendix  we  derive  the  formulae  (147)  and  (148)  of  the  main  text.  In  the  derivation  we  shall  use  the 
transformation  of  variables 

17=  cos  (B-1) 

and  the  numbers 

77S'^'  = cos(^^)  7=1,2,...,A/,  (B-2) 

where  M and  J are  positive  integers.  There  are  five  cases  to  be  considered.  The  method  of  derivation  is 
basically  that  of  Williams‘S. 

Cases  (i)  and  (U) 

The  function  T/Xr\)  is  defined  by  the  formula 

TA(y))=  J—fXis  A^l  (B-3) 

IT 

in  conformity  with  formula  (A-59).  The  function  Ta(t))  is  a polynomial  of  degree  A in  17  and  its  zeros 
p = 1, 2, . . . , A,  are  given  by 

;tp''’  = cos(^^^^  p=l,2,  ...,A,  (B-4) 


and 


^ iizk^L  V 

p-i  p-> 


(1-^7^)  I 


By  making  use  of  the  formula  (B-3)  for  T/i{rj)  we  establish  that 


and 


sinA(^ 

TaIv)  cosAd)sin</> 


IM 

niv) 


1 2 n(r?) 

Ta(v) 


cos 


<l> 


i- 


In  case  (i)  A is  given  by  (see  formula  (149)) 


A = a(M+l) 

where  a and  M are  any  positive  integers.  Then,  from  (B-11),  we  get 

, sin  [/A7r/(M+1)] 

Ta^tjS*^)  cos  [/Air/(M+ 1)]  sin  [/ir/(M+ 1)] 

sin  (aJir) 

^ cos  (oJtt)  sin  [Jn/{M+ 1)] 

= 0 

for  all  positive  integers  J from  1 to  M inclusive. 

In  case  (ii)  A is  given  by  (see  formula  (150)) 

A = i(2a-1)(M+1) 

where  a is  any  positive  integsr  and  M is  any  odd  positive  integer.  Then,  from  (B-11),  we  get 

l . sin(|(2fl-l)7Tr) 

Txiv?^  cos  (|(2a  - l)/7r)  sin  [J7r/(M+  1)] 

= 0 


(B-10) 


(B-11) 


(3-12) 


(B-13) 


(B-14) 


(B-15) 


(B-16) 


for  all  even  positive  integers  J from  2 to  (M  - 1)  inclusive.  From  (B-12)  we  then  get,  in  cases  (i)  and  (ii) 

.2 


n(T?n 

rA(T,n  {i-(t7S"y}' 


(B-17) 


Finally,  on  putting 


V -Vj 


fB-181 


into  formulae  (B-9)  and  (B-10),  using  results  (B-14),  (B-16)  and  (B-17),  and  rearranging  the  resulting 
expressions,  we  get 


(*o  1 f {i-Crnn 


= 0 


(B-19) 
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and 


, , 1 A , 


(B-20) 


in  both  cases  (i)  and  (ii).  It  must  be  remembered  however  that  in  case  (i)  the  quantity  M is  any  positive  integer 
and  J is  any  positive  integer  from  1 to  M inclusive,  whereas  in  case  (ii)  the  quantity  M is  any  odd  positive 
integer  and  J is  any  even  positive  integer  from  2 to  (Af  — 1)  inclusive. 

Formulae  (B-19)  and  (B-20)  agree  with  formulae  (147)  and  (148)  of  the  main  text  for  cases  (i)  and  (ii)  if  we 
take 


Pa  = A 


(B-21) 


because  the  p = 1, 2, . . . , A,  are  the  Gaussian  integration  points  corresponding  to  the  weight  function 

l/Vl-7,2. 


Case  (iii) 

We  define  the  function  Sa(i7)  by  means  of  the  formula 

5A(r7)  = cos(A  + 2)<^. 


(B-22) 


For  any  positive  integer  A it  is  the  product  of  vl  + 17  with  a polynomial  in  17  of  degree  A and  its  zeros  Op 
p = 1, 2, . . . , A,  are  given  by 


(B-23) 


which  is  in  conformity  with  (A-80).  We  note  that 


V -t 

L “p  -2- 

p-i 


(B-24) 


We  can  write  for  the  alternative  formula 


Sa{v)  = Sa  n (77-a{,''Vl  + n 


(B-25) 


where  5a  is  a constant  for  a given  integer  A.  From  (B-25)  and  the  equation  obtained  from  (B-25)  after 
differentiation  with  respect  to  77  we  get 


5^(77)  ^ 1 1 

s^M  P-.  (77 2(1  + 77)’ 


(B-26) 


and  on  differentiating  (B-26)  with  respect  to  77  we  get 


5^(77)  /5;(77)y  ^ 1 1 

5a(77)  V5a(77)/  pt,  (77-a<^V  2(1  + 77)" 


(B-27) 


Therefore 


A ..  7, 5^(77)  . 


(B-28) 
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For  any  positive  integer  A it  is  the  product  of  vl  — t;  with  a polynomial  in  r;  of  degree  A and  its  zeros 
p = 1, 2, ....  A,  are  given  by 

= cos  p = 1, 2, ....  A,  (B-40) 

which  is  in  conformity  with  (A-48).  We  note  that 

Z (B-41) 


We  can  write  for  R\{y\)  the  alternative  formula 


R^{r,)  = R^  n 


P“1 


(B-42) 


where  /?a  is  a constant  for  a given  integer  A. 

By  proceeding  in  exactly  the  same  manner  as  in  case  (iii)  from  formula  (B-25)  onwards  we  derive  the 
formulae 


and 


vr- 


1 ^ „ 


(B-43) 


1 + 


1 ^ 

(A+J)pri(,,S'^-/3^''y 


= A + ^ 


(B-44) 


provided  that  M is  any  even  positive  integer  and  J is  any  odd  positive  integer  from  1 to  (Af- 1)  inclusive. 
Formulae  (B-43)  and  (B-44)  agree  with  formulae  (147)  and  (148)  of  the  main  text  for  case  (iv)  if  we  take 

PA  = A-t-i,  (B-45) 

because  the  p = 1, 2, ....  A,  are  the  Gaussian  integration  points  corresponding  to  the  weight  function 
'/(l-i7o)/(l  + T>o). 

Case  (v) 

We  define  the  function  (?a(i7)  by  means  of  the  formula 


OA(r?)  = sin  (A-t-l)<^.  (B-46) 

For  any  positive  integer  A it  is  the  product  of  Vl  with  a polynomial  in  17  of  degree  A,  the  polynomial  being 
a multiple  of  y\(7f)  of  formula  (A-24).  The  zeros  p = 1, 2, ....  A,  of  C?a(»7)  are  given  by 


rjp''*  = cos 


p = 1, 2, ....  A, 


just  as  in  equation  (A-2S).  We  note  that 


Z n|,''’=0. 


We  can  write  for  C?A(rj)  the  alternative  formula 


0^{v)=0^  n (rj-rj'^yl^ 


P-1 


(B-47) 


(B-48) 


(B-49) 
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where  Of,  is  a constant  for  a given  integer  A.  From  (B-49)  and  the  equation  obtained  from  (B-49)  after 
differentiation  with  respect  to  17  we  get 


y —5— 

P-i  (Tj-rjh  (l-vy 


and  on  differentiating  (B-50)  with  respect  to  17  we  get 


Therefore 


and 


(0'a(v)y  f 1 (i+t7'‘) 

^OA(vy  p-i  (i-vY' 


P-1 


(1-r,^) 


QUt?) 

Oa(i7) 


+ ( A + 1 )t7 


(B-50) 


(B-51) 


(B-52) 


p-i  (i7-i7p''y 


q"a(v)  (O'Aivy 

iv)  ^Oa(v) 


ri 


+ 2»7 


QA(r7) 

Oa(tj) 


-(A  + 1). 


(B-53) 


By  making  use  of  the  formula  (B-46)  for  QaM  we  establish  that 


OaM^  (Ain 

OaM  ^ ^in(A+l)(Asin<f) 


and 


OX(v) 

Oa(v) 


(1 


In  case  (v)  A is  given  by  (see  formula  (153)) 

A = |(2a-l)(Af+l)-l 


(B-54) 


(B-55) 


(B-56) 


where  a is  any  positive  integer  and  A/  is  any  odd  positive  integer.  Then,  from  (B-54)  and  (B-55),  we  get 


O'Ainn 

04^) 


= 0 


(B-57) 


and 


QX(r?n  (A-H)^ 

OA(v^r)  {i-(7?ry} 


for  all  odd  positive  integers  / from  1 to  Af  inclusive.  Finally,  on  putting 


V = 


_(AO 

Vj 


(B-58) 


(B-59) 


into  formulae  (B-52)  and  (B-53),  using  results  (B-57)  and  (B-58),  and  rearranging  the  resulting  expressions, 
we  get 


„(Af) 

Vj 


1_  A {l-(virr}  n 

(A-H)p?,  (v'r-vi'^^ 


(B-60) 
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and 


provided  that  M is  any  odd  positive  integer  and  J is  any  odd  positive  integer  from  1 to  Af  inclusive. 
Formulae  (B-60)  and  (B-61)  agree  with  formulae  (147)  and  (148)  of  the  main  text  for  case  (v)  if  we  take 

Fa  = A + 1,  (B-62) 

because  the  p = 1, 2, ....  A,  are  the  Gaussian  integration  points  corresponding  to  the  weight  function 
v'l-tto. 


APPENDIX  C 


Noinerical  Evaluation  of  tbe  Function  i|,  i|o;  v,  Afw) 

Integration  with  respect  to  So 

In  this  Appendix  we  discuss  the  process  for  the  numerical  evaluation  of  the  function  tj,  tjo;  M*) 
from  the  formula 


/‘"’(f.  V,  Vo-,  i',  Woo)  (C- 1) 


where 


r,=^ 

s 


^o  = -7^[xo-x^.(yo)] 


c(yo) 


T. 

Vo 

s 


(C-2) 


-ittu/l 


du 


2Moc,(MooX  +R) 


, -/  f-jV/-x+MocJ?\) 


(C-3) 


and 


J?=V{x"  + (l-Mi)y"}.  (C-4) 

Formula  (C-1)  is  identical  with  formula  (114)  of  the  main  text. 

If  rjo  = 17  the  evaluation  of  the  integral  on  the  right  hand  side  of  (C-1)  is  straightforward  because,  from  (C-3) 
we  have 


and  therefore, 


X <0, 

x>0, 


(C-5) 


(C-6) 


The  integral  on  the  right  hand  side  of  (C-6)  can  be  evaluated  analytically  using  formulae  of  Appendix  A or  it 
can  be  evaluated  numerically.  To  evaluate  it  numerically,  divide  the  range  (0,  ^)  of  into  a number  of 
intervals  of  equal  length  and  apply  Gaussian  integration  formulae,  using  a small  number  of  integration  points, 
to  the  integral  over  each  of  these  intervals.  In  the  interval  abutting  on  fo  = 0 the  weight  function  1 /v^  is  used 
whereas  in  all  the  other  intervals  the  weight  function  1 is  used  for  the  Gaussian  integration.  The  number  of 
intervals  of  integration  which  it  is  necessary  to  take  to  obtain  a given  accuracy  will  depend  on  the  number  of 
integration  points  in  each  interval  and  on  the  value  of  n,  and  will  need  to  be  increased  if  n is  increased  because 
of  the  resulting  increase  in  the  number  of  undulations  of  the  function  hl"\(o).  The  actual  number  of  intervals 
required  can  be  assessed  by  experience. 

U Tto^V  fhe  range  (0, 1)  of  fo  in  the  integral  on  the  right  hand  side  of  (C-1)  is  divided  into  a number  of 
intervals,  not  necessarily  of  equal  length,  and  Gaussian  integration  formulae  using  a small  number  of 
integration  points  are  applied  to  the  integral  over  each  of  these  intervals.  In  the  interval  abutting  on  = 0 the 


weight  function  1/v^  is  used,  in  the  interval  abutting  on  = 1 the  weight  function  Vl— fo  is  used  whereas  in 
an  interval  completely  internal  to  the  range  (0, 1)  of  (o  the  weight  function  1 is  used  for  the  Gaussian 
integration.  At  least  two  intervals  of  integration  must  be  used  to  cover  the  range  (0, 1)  of  for  this  scheme.  If 
T/o  is  near  to  tj  the  intervals  of  integration  need  to  be  shorter  when  io  is  near  to  f than  when  is  far  from  i to 
attain  a given  accuracy  in  the  numerical  integration  because  of  the  rapid  change  in  values  of  [(y  - yo)//]  K[(x  - 
xq)/1,  (y  - yo)// ; v.  Moo]  and  its  derivatives  with  respect  to  when  fo  changes  in  the  neighbourhood  of  fo  = 
Also  the  nearer  tjo  is  to  tj  the  more  intervals  of  integration  for  in  (0, 1)  will  be  needed  to  attain  a given 
accuracy.  The  number  of  intervals  will  depend  on  n,  the  frequency  parameter  y and  the  Mach  number  M*. 

Expansions  for  the  Kernel  Function 

To  evaluate  numerically  the  integral  on  the  right  hand  side  of  (C- 1)  when  tjo  ^ 17,  we  must  be  able  to  obtain 
numerical  values  of 


for  a given  value  of  y ^ 0 and  a number  of  values  of  x.  We  shall  write 


where 


-r 


MUMccJc+R)y^ 

R(x^+y^) 


--‘yu/i 


(M^+y¥ 


(«"+!)* 


—X  +MaoR 
''  (1-Mi)  • 


(C-10) 


The  only  part  of  the  right  hand  side  of  (C-8)  that  is  at  all  complicated  to  evaluate  is  Jix/l,  y/\\v.  Moo).  The 
methods  of  evaluation  are  different  for  |X|/|y  | large  and  for  lX|/|y  | not  large.  We  shall  first  consider  |X|/|y  | to 
be  not  large. 

We  write  (C-9)  in  the  form 


4 f; 


-iv\y\u/l  r 

(u"+l)^  Jo 


(C-11) 


The  finite  integral  on  the  right  hand  side  of  (C- 1 1)  can  be  evaluated  numerically  in  a straightforward  manner. 
The  procedure  is  to  divide  the  interval  (0,  X/\y\)  of  the  integration  variable  u into  a number  of  subintervals  of 
equal  length,  the  length  of  each  subinterval  not  being  greater  than  any  of  the  quantities  1,  \X\/\y\  and  rr//r'|y  |. 
A Gaussian  numerical  integration  with  weight  function  1 is  then  carried  out  over  each  subinterval  using  a small 
number  of  integration  points.  This  procedure  is  satisfactory,  except  when  the  number  of  subintervals  becomes 
very  large,  as  happens  if  any  one  of  the  quantities  |X|/|y|  and  v\X\IitI  is  very  large.  We  do  not  anticipate 
dealing  with  such  a large  v that  v\X\/itI  is  so  large  that  the  above  procedure  is  unacceptable,  but  |X|/|y  | could 
be  large  enough  to  make  the  above  procedure  unacceptable.  We  shall  consider  the  demarcation  value  of 
|X|/|y|  later. 

Let  us  introduce  the  function  S(a),  for  a real,  by  means  of  the  formula 


, du 


(C-12) 
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Then  formula  (C-1 1)  can  be  written  as 


! 

i 

i 

) 


r 


( 

j 


(C-13) 

We  can  show  that  S(a)  satisfies  the  differential  equation 

fd"S(a)  J dSia)  . 

*1  da’  -*<“>!"  da  ■■ 

(C-14) 

• 

by  inserting  S(a)  from  (C-12)  into  (C-14).  If  we  put 

S(a)  = aP(o), 

(C-15) 

substitute  for  5(a)  from  (C-IS)  into  (C-14)  and  divide  the  resulting  equation  through  by  we 
differential  equation 

get  the 

P'(a)+-P'(«)-(l+-^W)  = T 

a \ a / a 

(C-16) 

for  P(a),  which  we  recognise  as  a modified  form  of  Bessel’s  differential  equation  with  a non-zero  right  hand 
side.  This  differential  equation  has  the  general  solution  (see  e.g.  Ref.  IS) 

P(a)  = ^{iLi(a)+CIt{a)+DK,(a)} 

p-  17) 

where  C and  D are  integration  constants,  Li(a)and  /iri(a)are  modified  Bessel  functions  of  order  1 and  of  the 
first  and  second  kinds  respectively  and  Li(a)  is  a modified  Struve  function  which  is  related  to  the  Struve 
function  //_i(ia)  by  means  of  the  formula  (see  Ref.  15) 

Liia)  = -H-x(ia). 

(C-18) 

For  small  values  of  |a|  we  may  write  (see  Ref.  15) 

Li(a)  = a]{a) 

(C-19) 

^i(« ) = ^ + «[  log  + r]  / («  ) - a/(a  ) 

(C-20) 

• 

Li(a)  = --g(a), 

17 

(C-21) 

where  /(a),  f(a)  and  g(a)  are  even  integral  functions  of  a which  have  the  power  series  expansions: 

(C-22) 

(C-23) 

and 

* 2^>'r' 

^^‘*^",?o(27)!(2r+l)!"  ’ 

(C-24) 
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For  real  positive  a the  branch  of  log  (a/2)  which  is  real  is  taken  in  formula  (C-20).  Therefore,  on  using  the 
above  expansions,  we  get,  for  small  values  of  a 


S(a)  = 2frD  — (ia)+4‘>>'Ca^  + 4vDa^^y  + log  — -j  + 0(a*).  (C-25) 


However,  we  can  obtain  an  expansion  of  S(a),  for  small  values  of  o,  directly  from  the  integral  representation 
(C-12).  In  fact,  for  small  real  positive  a,  we  can  show  that 


S(a)=l-(ia)+2a^j^>-+log^0-i+yj+O(cit’). 


(C-26) 


On  comparing  the  expansions  (C-25)  and  (C-26)  we  get  immediately 


(C-27) 


D = ^ 

TT 


(C-28) 


for  the  values  of  the  integration  constants  C and  D introduced  in  formula  (C-17). 

Having  obtained  the  values  of  the  integration  constants  C and  D we  can  substitute  for  P(a ) from  (C-17)  into 
(C-15)  to  get  S(a)  in  the  form 


S(a)  = F(a)+ iG(a) 


(C-29) 


where 


F(a  )=  1 -I-  y -(-  log  j a V(a ) - « V(« ) 


(C-30) 


G(a)  = -a^j(a)-ag{a) 


(C-31) 


are  real  functions  of  a for  real  positive  a. 

We  may  also  deduce  from  (C-12)  the  following  asymptotic  expansions  for  F(«)  and  G(a)of  formula  (C-29) 
for  large  real  positive  a.  With  the  integer  n^O  arbitrary  we  get 


F(a)=  \/^e-“jF„(a)-h^5„(a)} 


(C-32) 


G(a ) = G„  (a  ) + (a ) 


(C-33) 


where 


Fo(a)=l 


(C-34) 


(4-(2i-l)^)  n^l, 

,-tr!(8a)  »-i 


(C-35) 


1 " (2r)!(2r+f)! 


(C-36) 
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The  remainder  functions  S„(a),  e„(a),  for  any  n 3*0  have  the  behaviour 


5„(a)  = 0(l)  (C-37) 

e„(a)  = 0(l)  (C-38) 


for  a -»  +00. 

The  same  asymptotic  expansions  (C-32)  and  (C-33)  may  be  deduced  by  using  the  asymptotic  expansions  of 
Ii{a),  Ki(a)  and  Li(a),  as  given  for  example  in  Ref.  15,  in  formula  (C-17)  for  P(a)  with  C and  D given  by 
(C-27)  and  (C-28)  and  substituting  the  resulting  asymptotic  expansion  into  (C- 1 5).  It  is,  however,  quite  easy  to 
get  the  required  expansions  directly  from  (C-12). 

The  power  series  (C-22),  (C-23)  and  (C-24)  have  to  be  truncated  at  finite  values  of  r in  order  to  be  able  to 
evaluate  from  them  the  values  of  the  functions  j (a ),  f(a ) and  g(a ),  and  these  finite  values  of  r will  depend  upon 
the  accuracy  to  which  these  functions  are  required  and  on  the  value  of  a under  consideration.  If  we  work 
numerically  to  a given  number  of  significant  figures  the  accuracy  with  which  wc  can  evaluate  the  sums  of  the 
truncated  series  will  decrease  as  a increases.  Thus,  even  though  these  power  series  expansions  are  convergent 
for  any  finite  value  of  a,  they  cannot  be  used  to  give  accurate  values  of  the  functions  /(a ),  f(a ) and  g(a ) when  a 
becomes  indefinitely  large,  if  we  are  limited  in  the  number  of  significant  figures  used  in  the  arithmetical 
operations.  By  working  with  a given  number  of  significant  figures  there  is  a maximum  value  of  a for  which 
formula  (C-30),  with /(a)  and  /(a)  obtained  from  the  power  series  expansions  (C-22)  and  (C-23)  respectively, 
can  be  used  to  obtain  F{a)  to  within  some  prescribed  accuracy  e > 0.  Similarly  there  is  a maximum  value  of  a 
for  which  formula  (C-31),  with  y(a)  and  g(a)  obtained  from  the  power  series  expansions  (C-22)  and  (C-24) 
respectively,  can  be  used  to  obtain  C7(a)  with  the  accuracy  e.  The  smaller  a is,  the  fewer  terms,  in  general,  will 
be  needed  in  the  truncation  of  the  power  series  expansions  (C-22),  (C-23),  (C-24)  to  obtain  F(a)  and  G(a)  to 
within  the  accuracy  e. 

For  very  large  values  of  a we  can  use  the  asymptotic  formulae  (C-23)  and  (C-33)  to  evaluate  the  values  of 
F(a)  and  G(o).  Because  of  formula  (C-32)  we  can,  for  given  e >0,  find  Oi(n,  e)>0  such  that 


F(a)- 


e'‘F„(a) 


< e 


(C-39) 


whenever 


a>at(n,  e).  (C-40) 

If  we  take  a i (n,  e ) to  be  the  minimum  quantity  for  which  (C-39)  is  true  under  the  condition  (C-40)  we  find  that, 
for  fixed  e,ai(n,  e)  decreases  in  general  as  n is  increased  from  zero  up  to  a certain  value  of  n and  then  increases 
as  n is  increased  beyond  this  certain  value.  The  minimum  value  of  ai(n,  e),  for  all  the  values  of  n,  is  then  the 
minimum  value  of  a for  which  F(a)  may  be  obtained  to  within  accuracy  e from  formula  (C-32).  This  accuracy 
may  be  somewhat  reduced  if  we  work  numerically  to  a given  number  of  significant  figures.  The  higher  a is, 
beyond  the  minimum  value,  the  smaller  will  the  value  of  n need  to  be,  in  general,  for  (C-39)  to  be  true. 
Likewise,  because  of  formula  (C-33)  we  can,  for  given  c > 0,  find  a2(n,  c)>  0 such  that 

|G(a)-G„(a)|<e  (C-41) 


whenever 


a>02(n,  e).  (C-42) 

Again  there  is  a minimum  value  of  a for  which  G(a)  may  be  obtained  to  within  accuracy  e from  formula 
(C-33). 

The  maximum  values  of  a for  which  F(a)  and  G(a)  can  be  evaluated  to  the  given  accuracy  from  formulae 
(C-30)  and  (C-31)  with  /(a),  /(a)  and  g(a)  obtained  from  the  power  series  expansions  (C-22),  (C-23)  and 
(C-24)  respectively,  depends  strongly  on  the  number  of  significant  figures  used  in  the  arithmetic,  whereas  the 
minimum  values  of  a for  which  F(o)  and  G(a)  can  be  evaluated  to  the  given  accuracy  from  formulae  (C-32) 
and  (C-31)  is  hardly  dependent  on  the  number  of  significant  figures  used  in  the  arithmetic,  provided  that  this 
number  is  greater  than  the  number  of  significant  figures  required  in  the  values  of  the  functions.  The  functions 
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F(a)  and  G(a)  cannot  be  evaluated  to  the  given  accuracy  e,  for  all  a in  (0,  oo),  using  formulae  (C-30),  (C-31), 
(C-32)  and  (C-33)  as  described  above,  if  the  number  of  significant  figures  used  in  the  arithmetic  is  not 
sufficiently  high.  In  other  words,  if  we  work  to  a given  number  of  significant  figures,  then  e must  be  greater  than 
a certain  lower  bound  in  order  that  F{a ) and  G (o ) may  be  evaluated  in  the  above  manner  to  the  given  accuracy 
e for  all  a in  (0,  c»).  If  e is  less  than  this  bound  then  some  other  means  of  evaluating  the  functions  F(a)  and 
G(a)  must  be  used,  at  least  over  the  ranges  of  a for  which  the  above  method  does  not  yield  the  required 
accuracy  e.  We  shall  expand  the  functions  in  series  of  Chebyshev  polynomials  rather  than  power  series. 
Although  e will  still  have  to  be  greater  than  a certain  lower  bound  in  order  that  F{a)  and  G(a)  may  be 
evaluated  to  the  given  accuracy  e,  this  lower  bound  should  be  less  than  the  former  lower  bound,  thus  rendering 
the  procedure  involving  expansion  of  factions  in  series  of  Chebyshev  polynomials  of  wider  application  than 
that  involving  expansion  of  functions  in  power  series.  It  may  be  true  that  with  the  number  of  significant  figures 
available  on  a particular  computing  machine,  functions  may  be  evaluated  to  a sufficient  accuracy  for  some 
applications  using  expansions  in  power  series,  but  it  would  seem  to  be  good  practice  to  use  another  procedure 
which  is  capable  of  giving  superior  accuracy  and  is  no  more  difficult  to  apply. 

Use  of  Chebyshev  Polynomials 

The  formulae  (C-30),  (C-31),  (C-32)  and  (C-33)  are  valid  for  complex  values  of  a provided  the  branch  line 
of  log  (a/2)  lies  in  the  half-plane  Im  (a)>0,  but  we  are  here  concerned  only  with  real  positive  values  of  a. 
There  are  convergent  expansions  for  /(a),  /(a)  and  g(a),  which  are  valid  only  for  real  a but  which  are  more 
suitable  for  numerical  computation  than  are  (C-22),  (C-23)  and  (C-24)  when  a is  in  some  restricted  interval 
(0,  A)  where  A is  some  positive  finite  number.  We  can  also  use  expansions  for  F(a)  and  G(a)  which  are  valid 
only  for  real  a but  which  are  convergent  when  a is  in  the  restricted  range  (A,  oo)  as  opposed  to  the  expansions 
(C-32)  and  (C-33)  which  are  asymptotic  expansions.  These  convergent  expansions  in  series  of  orthogonal 
polynomials  are  again  more  suitable  for  numerical  computation  than  are  the  asymptotic  expansions.  Because 
of  their  simple  properties  we  shall  use  the  Chebyshev  polynomials  rather  than  other  orthogonal  polynomials, 
but  it  is  well  to  remember  that  other  orthogonal  polynomials  may  be  more  appropriate  to  use  in  some 
circumstances.  The  Chebyshev  polynomial  T„{x)  defined  by 

r„(Ar)  = cos  («  cos“’ x)  «=  0,1,2,...,  (C-43) 

is  a polynomial  of  degree  n in  x.  Because  j[a),  f(a)  and  g(a)  in  formulae  (C-30)  and  (C-31)  are  integral 
functions  of  a ^ we  can  write,  for  0«a  ^ A, 

F<«)- 1; DAA )r,(^- 1) log ( J)  I; l)  (C.44) 

and 

G(a)  = -a  I'  E,{A)Tri^- l)  + Ja^  £ a{A)z(^- l)  (C-45) 

where  the  dash ' on  the  summation  sign  X indicates  that  the  quantity  under  the  summation  sign  for  r = 0 is  to  be 
multiplied  by  i The  coefficients  C,(A),  D,{A)  and  E,{A)  in  the  formulae  (C-44)  and  (C-45)  may  be 
determined  numerically  for  r = 0, 1, 2, . . . , for  a given  value  of  A by  applying  the  method  of  Clenshaw'*. 

It  is  apparent  from  equations  (C-32)  and  (C-33)  that  the  function  c°F(a)/'Ja  is  a function  of  bounded 
variation  in  1/a  and  the  function  aG(a)  is  a function  of  bounded  variation  1/a^  for  large  a.  Hence  we  can 
write,  for  A « a « 00 


/=i:a)  = V^e-“  I'  A,(A)ri— -l)  (C-46) 

r-o  V a / 

and 

G(a)  = -!-  r B,(A)ri^-l),  (C-47) 

a r-o  'a  / 
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where,  again,  the  dash  ’ on  the  summation  sign  £ indicates  that  the  quantity  under  the  summation  sign  for  r = 0 
is  to  be  multiplied  by  2.  The  coefficients  Ar(A)  and  B,(A)  in  the  formulae  (C-46)  and  (C-47)  may  be 
determined  numerically  for  r = 0, 1, 2, . . . , for  a given  value  of  A by  applying  the  method  of  Clenshaw'*. 

If  we  determine  the  coefficients  A,(A ),  B,(A ),  Cr(A ),  Dr(A ) and  Er(A ) numerically  for  /4  = 7 we  get  the 
following  results; 


/io(7)=  2-57172 
Ai(7)=  0-032  21 
/42(7)= -0-000  32 
^3(7)=  0-000  01 


(C-48) 


Bo(7)  = 
B,(7)  = 
B2(7)  = 
B3(7)  = 

B,(7)  = 
Bs(l)  = 
Be{7)  = 

B7(7)  = 
B8(7)  = 
B^(7)  = 
Bio(7)  = 
Bu(7)  = 
Bu(7)  = 


2-080  503  0 
0-043  102  0 
0-002  271  0 
-0-000  790  2 
-0.000  144  5 
0-000  075  0 
-0-000  002  4 
-0-000  008  1 
0-000  003  5 
-0-000  000  4 
-0-000  000  4 
0-000  000  3 
-0-000  000  1 


Co(7)  = 

31-850 

997 

610 

^(7)  = 

20-478 

045 

274 

C2(7)  = 

6-700 

050 

944 

C,(7)  = 

1-298 

530 

059 

C4(7)  = 

0-164 

782 

643 

C5(7)  = 

0-014 

674 

957 

C,(7)= 

0-000 

964 

651 

Cy{7)  = 

0-000 

048 

629 

C*(7)  = 

0-000 

001 

937 

C,(7)  = 

0-000 

000 

062 

C,o(7)  = 

0-000 

000 

002 
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(C-49) 


(C-50) 


J 


Z?o(7)=  48-347  809  077 
£),(7)=  14-168  616  622 
Z?2{7)  = -19-622  120  487 
£>3(7)  = -14-015  111  195 
Z?4(7)=  -3-981  038  753 
05(7)=  -0-646  707  438 
06(7)=  -0-068  862  606 
0,(7)=  -0-005  195  562 
08(7)=  -0-000  292  608 
0,(7)=  -0-000  012  776 
O,o(7)=  -0-000  000  445 
0,i(7)=  -0-000  000  013 

£o(7)=  156-661  749  51 
E,(7)=  11 1-982  807  86 
E2(7)=  43-246  690  85 
£3(7)=  9-928  919  83 

£4(7)=  1-474  963  03 

£5(7)=  0-151515  57 

£6(7)  = 0-011  328  83 

£7(7)=  0-000  64165 

£8(7)=  0-000  028  41 

£9(7)=  0-000  00101 

£,o(7)=  0-000  000  03 


(C-51) 


(C-52) 


If  we  now  use  the  values  of  the  coefficients  A,(7),  £,(7),  Cr(7),  D,(l)  and  £,(7)  from  (C-48),  (C-49),  (C-50), 
, (C-51)  and  (C-52)  in  formulae  (C-44),  (C-45),  (C-46),  (C-47),  and  neglect  the  remaining  higher  order 

coefficients,  we  can  evaluate  F(a)  and  G(a)  for  any  given  value  of  a.  By  doing  this  for  all  the  integers  a from  0 
to  25,  using  1 1 significant  figures  in  the  arithmetic,  we  find  the  values  £(a)  and  G(a)  of  Table  C-1  (overleaf). 

The  function  values  tabulated  in  Table  C-1  should  all  be  correct  to  seven  decimal  places.  To  achieve  this 
accuracy  the  coefficients  >4,(7)  needed  to  be  given  only  to  five  decimal  places  because  of  the  factor  -Ja  e”“  in 
formula  (C-46),  and  the  coefficients  £,(7)  needed  to  be  given  only  to  seven  decimal  places  because  of  the  factor 
1/a  in  formula  (C-47).  The  coefficients  C,(7),  D,(7)  and  £,(7)  could  have  been  given  to  one  fewer  decimal 
place  each,  but  there  is  not  much  gain  in  this  because  these  coefficients  reduce  so  rapidly  as  r is  increased. 

It  is  also  possible  to  evaluate  F(a ) to  seven  places  of  decimals  for  a = 7 from  the  asymptotic  approximation 
•J trail  e“"£„(o)  with  « = 5.  For  higher  values  of  a the  value  of  n required  to  get  this  accuracy  may  be  less  than 
5.  However  it  is  possible  to  evaluate  G(a)  accurate  to  only  three  places  of  decimals  for  o = 7 from  the 
asymptotic  approximation  G„(a)  and  to  get  the  highest  accuracy  for  a = 7 we  must  take  n = 2.  Thus  with  the 
demarcation  value  >4  = 7 we  must  turn  to  some  other  formula,  such  as  (C-47),  in  order  to  be  able  to  evaluate 
G(a)  to  seven  places  of  decimals  for  a ^ 7.  The  asymptotic  approximation  G„(a)  cannot  be  used  to  evaluate 
G(a)  to  seven  decimal  places  unless  a > 17. 


If  F(a)  and  G(o)  were  required  to  higher  accuracy  than  seven  decimal  places  then  the  demarcation  value  >4 
in  formulae  (C-44),  (C-45),  (C-46)  and  (C-47)  would  have  to  be  taken  to  be  some  value  less  than  7,  unless  the 


TABLE  C-1 


a 

F(a) 

G(a) 

0 

1- 000  000  0 

0-000  000  0 

1 

0-601  907  2 

-0-468  450  8 

2 

0-279  731  8 

-0-467  289  0 

3 

0-120  469  3 

-0-376  343  2 

4 

0-049  934  0 

-0-291  743  6 

5 

0-020  223  1 

-0-229  284  5 

6 

0-008  063  5 

-0-185  609  1 

7 

0-003  179  3 

-0-154  965  3 

8 

0-001  243  0 

-0-132  881  5 

9 

0-000  482  7 

-0-116  408  0 

10 

0-000  186  5 

-0-103  692  7 

11 

0-000  071  7 

-0-093  577  4 

12 

0-000  027  5 

-0-085  324  8 

13 

0-000  010  5 

-0-078  451  3 

14 

0-000  004  0 

-0-072  629  2 

15 

0-000  001  5 

-0-067  628  7 

16 

0-000  000  6 

-0-063  283  7 

17 

0-000  000  2 

-0-059  471  1 

18 

0-000  000  0 

-0-056  097  1 

19 

0-000  000  0 

-0-053  089  4 

20 

0-000  000  0 

-0-050  390  6 

21 

0-000  000  0 

-0-047  955  1 

22 

0-000  000  0 

-0-045  745  8 

23 

0-000  000  0 

-0-043  732  4 

24 

0-000  000  0 

-0-041  889  7 

25 

0-000  000  0 

-0-040  196  9 

number  of  significant  figures  used  in  the  arithmetic  were  increased  beyond  11.  The  coefficients  A,(A)  and 
Br(A ) would  decrease  more  slowly  as  r increased  and  the  coefficients  CAA ),  DAA ) and  £,(A ) would  decrease 
more  rapidly  as  r is  increased  with  this  lower  value  of  A than  was  the  case  with  A=l.  For  high  enough 
required  accuracy  it  will  no  longer  be  possible  to  obtain  F(a)  from  the  asymptotic  approximation 
A-najl  e^“F„(a)  at  the  demarcation  value  a = A. 

The  numbers  of  terms  which  need  to  be  retained  in  the  infinite  series  in  (C-44),  (C-45),  (C-46)  and  (C-47) 
depend  only  on  A and  not  on  the  value  of  a.  On  the  other  hand,  the  numbers  of  terms  which  need  to  be 
retained  in  the  infinite  series  (C-22),  (C-23)  and  (C-24)  do  depend  on  a and  are  very  small  when  a is  very 
small.  Thus  for  very  small  a it  is  less  work  numerically  to  evaluate  F(a ) and  G(a ) from  (C-30)  and  (C-3 1 ) using 
the  series  (C-22),  (C-23)  and  (C-24)  fory(a),/(o)  and  g(a),  than  it  is  to  evaluate  F(a)  and  G(a)  from  (C-44) 
and  (C-45),  but  over  the  whole  range  (0,  A)  the  formulae  (C-44)  and  (C-45)  are  more  economical. 

The  numerical  evaluation  of  the  summations  in  formulae  (C-44),  (C-45),  (C-46)  and  (C-47)  is  easily  carried 
out  by  using  the  scheme  described  by  Clenshaw''’. 

From  (C-43)  we  get,  on  using  elementary  properties  of  the  cosine  function, 


To(x)=l  (C-53) 

r,(jt)=jc  (C-54) 


and  the  reduction  formula 


T„^.i(x)-2xT„(x)  + T„-,(x)  = 0. 


(C-55) 


68 


Suppose  now  that  we  wish  to  evaluate  k(x)  where  k(x)  is  given  by  the  formula 


k{x)=  X'  ajrix). 

r-0 


(C-56) 


Clenshaw’s  scheme  is  to  put 


bn+2  — 0 

bii+i  ~ 0 

bn—r  f +1  bn—r+2  "1”  ^n~ri  r 0,  1,  2,  . . . , /I. 


(C-57) 

(C-58) 

(C-59) 


k{x)  = kbo-b2) 


(C-60) 


as  may  be  shown  by  application  of  (C-53),  (C-54)  and  (C-55). 

The  function  S(a)  is  now  obtained  from  (C-29)  and  the  result  used  in  (C-13)  to  obtain  /(at//,  y/l ; v,  Moc) 
when  |X|/|y|  is  not  large. 

Procedures  when  |X|/|y|  is  Large 

To  evaluate  Jix/l,  y/l ; v,  Mco)  when  X/jy  | is  large  and  positive  we  replace  1/(m  ^ + 1)^  in  the  integrand  on  the 
right  hand  side  of  (C-9)  by  its  expansion  as  a power  series  in  1/m  and  integrate  term  by  term  to  get 


A/’/’  ’A^“j  2'>!r!  Ixj  / j 


(C-61) 


Jk  OO  1 

e”"’— 7,  r = 3,  5, 7, .... 

a V 


(C-62) 


The  expansion  on  the  right  hand  side  of  (C-61)isconvergentforX/|y|s  1 and,  the  larger  X/|y|  is,  the  faster 
does  the  expansion  (C-61)  converge,  i.e.  the  fewer  terms  in  a finite  truncation  of  it  are  necessary  to  get 
J(x/l,  y/l ; p,  Mco)  to  a given  accuracy.  We  may  arbitrarily  choose  to  use  the  expansion  on  the  right  hand  side  of 
(C-61)  for  obtaining  J(x/l,  y/l ; p,  Mco)  when  X/\y  1 3=  2. 

To  evaluate  £2r+3(«)  we  first  use  the  recurrence  relationships 


:2p>3(a)  [np  + 2) 


(2p  + 2)  (2p  + 2X2p  + l)J  (2p-l-2)(2p  + l) 


(2p-l-2)(2p  + l)^^'’''‘^“^  ^ 


to  express  E2r+3(o)  in  terms  of  £3(0). 

We  may  express  £3(0)  by  means  of  the  formula 


£3(«  ) = 2«^^log  a + -I-  H(a  ) - laK  {a  ), 


(C-64) 


where  the  branch  of  the  function  log  a which  is  real  for  a real  positive  is  taken  and  H(a)  and  K(a)  are  even 
integral  functions  of  a.  We  may  deduce  the  power  series  expansions 


fi{a)  = 2-2i2-y)a^-2a*  I 


; (-l)g^* 

o(s  + l)(2s  + 4)! 


(C-65) 


and 


K(a)  \+a  I^(2j  + ix2s  + 3)! 


(C-66) 


from  the  integral  representation  (C-62)  with  r = 3 for  Ei(a). 

For  large  values  of  o the  power  series  expansions  (C-65)  and  (C-66)  are  of  little  value.  It  is  then  better  to 
represent  E3(a)  by  means  of  the  formula 


E3(a)  = e ‘“{h(a)+ik(a)}. 


(C-67) 


We  may  then  deduce  from  (C-62)  with  r = 3 the  following  asymptotic  expansions  for  /i(a)  and  k(a)  for  real 
positive  a.  With  the  integer  n ^ 0 arbitrary  we  get 


h(a)=  /in(a)+-j;r^/i„(a) 
a 


(C-68) 


and 


fc  (a ) = (a ) + I'n  (a  ) 
a 


(C-69) 


where 


1 " (-l)-(2r  + 3)! 

hn{a)—  T 2 L 2r 

2a  r-o  01 


(C-70) 


and 


. _ 1 ^ (-l)'(2r  + 2)! 

kn(a)=-T-  h Tr • 

2a  r-o  a 


(C-71) 


The  remainder  functions  v„(a),  for  any  « >0  have  the  behaviour 

Ain(a)  = 0(l) 

*'»(«)  = 0(1) 


(C-72) 

(C-73) 


for  a -»  +00. 

Because  //(a)  and  ^r(a)  in  (C-64)  are  integral  functions  of  a^  we  can  write,  for  0=sa 

where  the  dash ' on  the  summation  sign  X indicates  that  the  quantity  under  the  summation  sign  for  r = 0 is  to  be 
multiplied  by  i The  coefficients  A/,(i4)  and  ^,(^4)  in  the  formula  (C-74)  may  be  determined  numerically  for 
r = 0, 1, 2, . . . , for  a given  value  of  A by  applying  the  method  of  Clenshaw*®. 

The  functions  a^h(a)  and  ak(a),  where  h(a)  and  it  (a)  are  defined  in  formula  (C-67),  are  functions  of 
bounded  variation  in  1/a^  for  large  a.  Hence  we  can  write,  for  >4  «a  «oo 

e3(«)  = ^ r Er(A )tJ^- l)  +^  r GAA )t{^- l)  (C-75) 

where,  again,  the  dash ' on  the  summation  sign  X indicates  that  the  quantity  under  the  summation  sign  for  r = 0 
is  to  be  multiplied  by  i The  coefficients  F,(A)  and  G,(,A)  in  the  formulae  (C-75)  may  be  determined 
numerically  for  r = 0, 1,  2, . . . , for  a given  value  of  A by  applying  the  method  of  Clenshaw'*. 
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If  we  determine  the  coefficients  F,{A ),  G,{A ),  M,{,A ) and  Nr(,A ) numerically  tor  A =1  we  get  the  following 
results: 


Fo(7) 

F,(7) 

F2CI) 

F3(7) 

F4(7) 

F,{1) 

F^d) 

F-rd) 

Fg(7) 

F9d) 

F,od) 


= +5-180  814 
= -0-360  529 
= +0-041  133 
= -0-006  359 
= +0-001210 
= -0-000  268 
= +0-000  067 
= -0-000  018 
= +0-000  005 
= -0-000  002 
= +0-000  001 


(C-76) 


) = -l-8198580 
)= +0-081  768  8 
)= -0-007  181  8 
>=+0-000  926  5 
>=-0-000153  8 
>=+0-000  030  6 
> = -0-000  007  0 
>=+0-000  0017 
>=-0-000  000  5 
>=+0-000  0001 


(C-77) 


Mo(7) 

M.(7) 

M2(7) 

M3(7) 

M4(7) 

Msd) 

M^d) 

Mid) 

M^d) 

M,(7) 


= +4-768  526  27 
= -0-677  567  89 
= -2-089  539  69 
= +0-398  136  24 
= -0-065  821  07 
= +0-007  658  89 
= -0-000  634  37 
= +0-000  038  88 
= -0-000  001  83 
= +0-000  000  07 


(C-78) 


No(7)= +7-151  718  43 
N,(7)= +2-201  177  93 
N2(7)= -0-310  914  27 
^^3(7)= +0-055  598  29 
N4(7)= -0-007  412  33  . 
N5(7)= +0-000  704  63 
^6(7)= -0-000  049  07 
^7(7)= +0-000  002  59 
^8(7)= -0-000  000  11 , 


(C-79) 


If  now  we  use  the  values  of  the  coefficients  iv(7),  G,(7),  Afr(7)and  A/r(7)from  (C-76),  (C-77),  (C-78),  (C-79) 
in  formulae  (C-74)  and  (C-75),  and  neglect  the  remaining  higher  order  coefficients,  we  can  evaluate  £3(0)  for 
any  given  value  of  a.  By  doing  this  for  all  the  integers  a from  0 to  25,  using  nine  significant  figures  in  the 
arithmetic,  we  find  the  values  Ei(a)  of  Table  C-2. 

TABLE  C-2 


a 

Re  £3(0) 

Im  £3(0) 

0 

+0-500  000  0 

+0-000  000  0 

1 

+0-018117  6 

-0-378  530  2 

2 

-0-271  409  2 

-0-107  735  2 

3 

-0-168  342  2 

+0-164  075  8 

4 

+0-058  9296 

+0-186  434  0 

5 

+0-163  769  9 

+0-031  120  2 

6 

+0-093  301  2 

-0-110  845  2 

7 

-0-043  467  6 

-0-120  258  2 

8 

-0-112  298  8 

-0-021  174  8 

9 

-0-068  523  3 

+0-077  155  1 

10 

+0-027  748  1 

+0-089  804  8 

11 

+0-083  589  3 

+0-021  269  7 

12 

+0-057  204  0 

-0-055  431  2 

13 

-0-015  793  7 

-0-072  263  6 

14 

-0-065  040  0 

-0-023  107  1 

15 

-0-0506515 

+0-040  227  3 

16 

+0-006  772  4 

+0-060  467  0 

17 

+0  051729  4 

+0-024  944  2 

18 

+0-046  076  9 

-0-028  863  5 

19 

+0-000  160  8 

-0-051616  2 

20 

-0-041  447  3 

-0-026  368  1 

21 

-0-042  393  6 

+0-019  961  0 

22 

-0-005  685  7 

+0-044  4416 

23 

+0-033  085  1 

+0-027  298  5 

24 

+0-0391214 

-0-012  760  0 

25 

+0-009  811  7 

-0-038  304  3 

The  function  values  tabulated  in  Table  C-2  should  all  be  correct  to  seven  decimal  places.  To  achieve  this 
accuracy  the  coefficients  F,(7)  needed  to  be  given  only  to  six  decimal  places  because  of  the  factor  e~'“/a^ 
multiplying  the  first  series  on  the  right  of  (C-75)  and  the  coefficients  G,(7)  needed  to  be  given  only  to  seven 
decimal  places  because  of  the  factor  e'‘”/o  multiplying  the  second  series  or.  the  right  of  (C-75). 

It  is  possible  to  evaluate  £3(0 ) only  to  three  places  of  decimals  for  a = 7 using  the  asymptotic  approximation 
hn (a)  to  h(a)  and  k„  (a)tok(a) and  to  get  the  highest  accuracy  at  or  = 7 we  must  take  n = 2 for  evaluating  both 
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h(a)  and  k (a ).  Thus  with  the  demarcation  value  -4  = 7 we  must  turn  to  a formula,  such  as  (C-75 ),  in  order  to  be 
able  to  evaluate  £3(0)  to  seven  places  of  decimals  for  a ^ 7.  The  asymptotic  approximations  /»,(o)  and  kn(a) 
cannot  be  used  to  get  £3(0)  to  seven  places  of  decimals  unless  a > 17. 

We  may  now  evaluate  /(jt/ /,  y//;  i>.  A/®)  from  the  infinite  series  (C-61 ) when  A/|y  | > 2 if  we  use  the  relations 
(C-63)  to  express  the  occurring  in  (C-61 ),  in  terms  of  Ei(yXll)  and  then  evaluate  Ei(vXH)  from 

either  formula  (C-74)  or  (C-75)  depending  on  whether  vXjl  < i4  or  vXjl  > A. 

To  evaluate  /(x//,  yll\v.  A/®)  when  X/\y  \ is  large  and  negative,  we  write,  from  (C-9), 


Xf-r 


du 


-[ 


X/lvl 


du 


(u^+\f 


g'-lvlu/l 


du 


(C-80) 


The  evaluation  of  aKi(a)  for  0«a  «cx)  has  already  been  considered.  The  evaluation  of 


du 

(M^  + 1)* 


for  -Xl\y\ > 2 is  carried  out  in  exactly  the  same  manner  as  the  evaluation  of  J{x/l,  y/l-,  v,  Af®)  for  Xl\y\ > 2 
was  carried  out,  apart  from  a change  in  the  sign  of  the  imaginary  part. 
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TABLE  1 


Nunerkal  Values  of  ApproximatioBS  <)iy,  I = 1, 2, 3;  y = 1, 2, 3,  to  die  Geaeralised  Airforces  on  a Tapered 
Swept  Winf  of  Aspect  Ratio  2,  Oscillating  in  Heave,  Pitch  and  Control  Snrface  Rotation  with  v - 0*32560, 

A#oc  = 0*7806 


m 

n 

M 

N 

d;. 

OT. 

On 

d'i3 

15 

2 

15 

2 

1 

0*056  532 

-2*4867 

-2*4586 

-4*5856 

-0*59672 

0*047  442 

15 

2 

15 

2 

2 

0059  524 

-2*5112 

-2*4806 

-4*6469 

-0*59274 

0*058  493 

15 

2 

15 

2 

4 

0059  303 

-2*5094 

-2*4791 

-4*6423 

-0*59301 

0*058  082 

15 

2 

15 

2 

8 

0*059  245 

-2*5092 

-2*4790 

-4*6412 

-0*59296 

0*058  138 

15 

3 

15 

3 

1 

0*063  709 

-2*5395 

-2*5066 

-4*7566 

-0*59599 

0*092  102 

15 

3 

15 

3 

2 

0*060  337 

-2*5124 

-2*4836 

-4*6788 

-0*58650 

0*084  073 

15 

3 

15 

3 

4 

0*061  863 

-2*5223 

-2*4921 

-4*7058 

-0*58572 

0*083  562 

15 

3 

15 

3 

8 

0*061  750 

-2*5214 

-2*4913 

-4*7032 

-0*58585 

0*083  567 

15 

4 

15 

4 

1 

0*072  016 

-2*5627 

-2*5224 

-4*8710 

-0*59746 

0*081  139 

15 

4 

15 

4 

2 

0*061  005 

-2*5217 

-2*4921 

-4*7037 

-0*58622 

0*087  618 

15 

4 

15 

4 

4 

0*061  854 

-2*5220 

-2*4921 

-4*7079 

-0*58420 

0*084  590 

15 

4 

15 

4 

8 

0*061  849 

-2*5227 

-2*4927 

-4*7084 

-0*58426 

0*084  585 

15 

5 

15 

5 

1 

0*069  369 

-2*5493 

-2*5125 

-4*8268 

-0*58719 

0*085  382 

15 

5 

15 

5 

2 

0*065  804 

-2*5421 

-2*5080 

-4*7790 

-0*58680 

0*082  546 

15 

5 

15 

5 

4 

0*060  595 

-2*5179 

-2*4891 

-4*6924 

-0*58435 

0*082  547 

15 

5 

15 

5 

8 

0*062  152 

-2*5243 

-2*4940 

-4*7140 

-0*58371 

0*085  162 

15 

6 

15 

6 

1 

0*068  997 

-2*5496 

-2*5130 

-4*8208 

-0*58459 

0*083  682 

15 

6 

15 

6 

2 

0*067  449 

-2*5440 

-2*5086 

-4*7980 

-0*58530 

0*082405 

15 

6 

15 

6 

4 

0*061  730 

-2*5239 

-2*4938 

-4*7125 

-0*58410 

0*084  300 

15 

6 

15 

6 

8 

0*061  857 

-2*5225 

-2*4926 

-4*7093 

-0*58364 

0*084  845 

15 

7 

15 

7 

1 

0*068  974 

-2*5507 

-2*5142 

-4*8176 

-0*58518 

0*085  373 

15 

7 

15 

7 

2 

0*065  972 

-2*5380 

-2*5042 

-4*7755 

-0*58626 

0*082  964 

15 

7 

15 

7 

4 

0*063  669 

-2*5315 

-2*4996 

-4*7419 

-0*58488 

0*084  483 

15 

7 

15 

7 

8 

0*061  305 

-2*5207 

-2*4912 

-4*7021 

-0*58346 

0*085  272 

15 

8 

15 

8 

1 

0*067  845 

-2*5448 

-2*5095 

-4*7959 

-0*58257 

0*087  852 

15 

8 

15 

8 

2 

0*065  255 

-2*5356 

-2*5024 

-4*7649 

-0*58684 

0*085  725 

15 

8 

15 

8 

4 

0*064  781 

-2*5348 

-2*5019 

-4*7571 

-0*58619 

0*085  461 

15 

8 

15 

8 

8 

0*061416 

-2*5218 

-2*4921 

-4*7052 

-0*58416 

0*084  262 

15 

9 

15 

9 

1 

0*067  288 

-2*5435 

-2*5087 

-4*7863 

-0*57992 

0*086  495 

15 

9 

15 

9 

2 

0*065  524 

-2*5388 

-2*5053 

-4*7699 

-0*58603 

0*086  130 

15 

9 

15 

9 

4 

0*064  664 

-2*5336 

-2*5010 

-4*7544 

-0*58664 

0*085  478 

15 

9 

15 

9 

8 

0*062  050 

-2*5245 

-2*4943 

-4*7152 

-0*58462 

0*084  989 

15 

10 

15 

10 

1 

0*066  940 

-2*5415 

-2*5070 

-4*7783 

-0*57914 

0*085  793 

15 

10 

15 

10 

2 

0*065  865 

-2*5382 

-2*5045 

-4*7711 

-0*58455 

0*085  210 

15 

10 

15 

10 

4 

0*064  066 

-2*5307 

-2*4987 

-4*7450 

-0*58584 

0*084  590 

15 

10 

15 

10 

8 

0*062  761 

-2*5272 

-2*4963 

-4*7258 

-0*58466 

0*084  965 

15 

2 

15 

3 

1 

0*058  199 

-2*5305 

-2*5003 

-4*6620 

-0*59146 

0*049  011 

15 

2 

15 

3 

2 

0*056910 

-2*4996 

-2*4713 

-4*6048 

-0*58893 

0*043  323 

15 

2 

15 

3 

4 

0*058  384 

-2*5112 

-2*4814 

-4*6337 

-0*58732 

0*049  199 

15 

2 

15 

3 

8 

0*058  298 

-2*5103 

-2*4806 

-4*6313 

-0*58745 

0*048  993 
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TABLE  1 — (continued) 


m 

n 

M 

N 

<1 

dll 

Oh 

d'i2 

Oh 

dh 

15 

2 

15 

4 

1 

0 066  883 

-2-5526 

-2-5159 

-4-7788 

-0-59040 

0-072  441 

15 

2 

15 

4 

2 

0056  460 

-2-5081 

-2-4799 

-4-6138 

-0-58954 

0-044  034 

15 

2 

15 

4 

4 

0058  345 

-2-5092 

-2-4796 

-4-6308 

-0-58752 

0-048  557 

15 

2 

15 

4 

8 

0-058  262 

-2-5101 

-2-4806 

-4-6310 

-0-58735 

0-048  878 

15 

3 

15 

4 

1 

0-071518 

-2-5608 

-2-5219 

-4-8634 

-0-59576 

0-083  508 

15 

3 

15 

4 

2 

0-060  772 

-2-5210 

-2-4917 

-4-6987 

-0-58872 

0-083  986 

15 

3 

15 

4 

4 

0-061712 

-2-5211 

-2-4914 

-4-7036 

-0-58510 

0-081  262 

15 

3 

15 

4 

8 

0-061  692 

-2-5219 

-2-4922 

-4-7042 

-0-58510 

0-081  363 

15 

3 

15 

5 

1 

0-067  870 

-2-5428 

-2-5078 

-4-8026 

-0-59178 

0-083  598 

15 

3 

15 

5 

2 

0-065  649 

-2-5406 

-2-5070 

-4-7738 

-0-59044 

0-082  491 

15 

3 

15 

5 

4 

0-060  261 

-2-5166 

-2-4882 

-4-6850 

-0-58587 

0-082  290 

15 

3 

15 

5 

8 

0-061  963 

-2-5231 

-2-4932 

-4-7086 

-0-58504 

0-081  277 

15 

3 

15 

6 

1 

0-067  323 

-2-5448 

-2-5100 

-4-7963 

-0-58833 

0-083  298 

15 

3 

15 

6 

2 

0-066  517 

-2-5394 

-2-5053 

-4-7812 

-0-59007 

0-081  535 

15 

3 

15 

6 

4 

0-061  439 

-2-5229 

-2-4932 

-4-7058 

-0-58691 

0-082  380 

15 

3 

15 

6 

8 

0-061  634 

-2-5213 

-2-4917 

-4-7031 

-0-58509 

0-081406 

15 

4 

15 

5 

1 

0-069  064 

-2-5449 

-2-5081 

-4-8175 

-0-59207 

0-082  452 

15 

4 

15 

5 

2 

0-065  842 

-2-5414 

-2-5072 

-4-7783 

-0-59036 

0-082  924 

15 

4 

15 

5 

4 

0-060  581 

-2-5176 

-2-4888 

-4-6917 

-0-58422 

0-086  010 

15 

4 

15 

5 

8 

0-062  137 

-2-5240 

-2-4938 

-4-7134 

-0-58434 

0-084  148 

m 

n 

M 

N 

« 

Oh 

^22 

0^3 

023 

15 

2 

15 

2 

1 

-0-083  903 

2-7077 

2-6375 

5-7238 

0-95694 

0-050  595 

15 

2 

15 

2 

2 

-0-088  217 

2-7098 

2-6344 

5-7720 

0-95471 

0-044  836 

15 

2 

15 

2 

4 

-0-087  931 

2-7090 

2-6340 

5-7680 

0-95500 

0-045  119 

15 

2 

15 

2 

8 

-0-087  884 

2-7086 

2-6337 

5-7667 

0-95494 

0-045  098 

15* 

3 

15 

3 

1 

-0-100  36 

2-7431 

2-6537 

5-9967 

0-98120 

-0-028  454 

15 

3 

15 

3 

2 

-0-094  451 

2-7286 

2-6470 

5-9099 

0-96384 

-0-027  201 

15 

3 

15 

3 

4 

-0-096  970 

2-7265 

2-6419 

5-9322 

0-96210 

-0-026  717 

15 

3 

15 

3 

8 

-0-096  806 

2-7259 

2-6415 

5-9295 

0-96234 

-0-026  787 

15 

4 

15 

4 

1 

-0-111  53 

2-7054 

2-6033 

6-0884 

0-97587 

-0-019  808 

15 

4 

15 

4 

2 

-0-09o315 

2-7395 

2-6555 

5-9499 

0-96552 

-0-033  136 

15 

4 

15 

4 

4 

-0-097  049 

2-7281 

2-6436 

5-9408 

0-95943 

-0-030  097 

15 

4 

15 

4 

8 

-0-097  155 

2-7285 

2-6439 

5-9415 

0-95952 

-0-029  990 

15 

5 

15 

5 

1 

-0-106  49 

2-7071 

2-6116 

6-0428 

0-96472 

-0-024  649 

15 

5 

15 

5 

2 

-0  103  00 

2-7259 

2-6341 

6-0109 

0-96199 

-0-025  359 

15 

5 

15 

5 

4 

-0-095  435 

2-7354 

2-6528 

5-9339 

0-95856 

-0-029  318 

15 

5 

15 

5 

8 

-0-097  602 

2-7283 

2-6431 

5-9469 

0-95864 

-0-030  893 

15 

6 

15 

6 

1 

-0-106  64 

2-7106 

2-6148 

6-0408 

0-95871 

-0-025  191 

15 

6 

15 

6 

2 

-0-104  82 

2-7152 

2-6215 

6-0219 

0-95983 

-0-024  189 

15 

6 

15 

6 

4 

-0-097  225 

2-7344 

2-6495 

5-9530 

0-95930 

-0-029  144 

15 

6 

15 

6 

8 

-0-097  122 

2-7291 

2-6446 

5-9437 

0-95848 

-0-030  833 

.75 
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TABLE  1 — (continued) 


m 

n 

M 

N 

<? 

Oil 

OS. 

diz 

^22 

Oi3 

0-23 

15 

7 

15 

7 

1 

-0- 106  55 

2-7060 

2-6104 

6-0281 

0-95920 

-0-028  807 

15 

7 

15 

7 

2 

-0  102  44 

2-7191 

2-6283 

6-0024 

0-96111 

-0-026  344 

15 

7 

15 

7 

4 

-0  099  863 

2-7280 

2-6400 

5-9760 

0-96023 

-0-028  960 

15 

7 

15 

7 

8 

-0096  415 

2-7320 

2-6483 

5-9398 

0-95872 

-0-031  418 

15 

8 

15 

8 

1 

-0- 104  87 

2-7059 

2-6125 

6-0057 

0-95623 

-0-032  553 

15 

8 

15 

8 

2 

-0- 101 72 

2-7223 

2-6323 

5-9956 

0-96328 

-0-030  357 

15 

8 

15 

8 

4 

-0101  25 

2-7229 

2-6334 

5-9867 

0-96242 

-0-029  838 

15 

8 

15 

8 

8 

-0  096685 

2-7328 

2-6488 

5-9438 

0-95905 

-0-030459 

15 

9 

15 

9 

1 

-0-104  17 

2-7066 

2-6140 

5-9944 

0-95113 

-0-031  572 

15 

9 

15 

9 

2 

-0-102  14 

2-7215 

2-6310 

5-9965 

0-96176 

-0-030  934 

15 

9 

15 

9 

4 

-0-100  92 

2-7222 

2-6332 

5-9836 

0-96308 

-0-030  275 

15 

9 

15 

9 

8 

-0-097  600 

2-7312 

2-6460 

5-9522 

0-96010 

-0-030  638 

15 

10 

15 

10 

1 

-0-103  63 

2-7053 

2-6134 

5-9838 

0-94942 

-0-031  183 

15 

10 

15 

10 

2 

-0-102  47 

2-7172 

2-6265 

5-9936 

0-95880 

-0-029  849 

15 

10 

15 

10 

4 

-0-100  07 

2-7240 

2-6360 

5-9764 

0-96142 

-0-029  413 

15 

10 

15 

10 

8 

-0-098  555 

2-7286 

2-6423 

5-9602 

0-96007 

-0-030  257 

15 

2 

15 

3 

1 

-0-086112 

2-7368 

2-6638 

5-7803 

0-94719 

0-049  598 

15 

2 

15 

3 

2 

-0-083  509 

2-7138 

2-6440 

5-7246 

0-94308 

0-056  959 

15 

2 

15 

3 

4 

-0-086  132 

2-7129 

2-6397 

5-7488 

0-94305 

0-053  779 

15 

2 

15 

3 

8 

-0-086  013 

2-7122 

2-6392 

5-7466 

0-94319 

0-053  965 

15 

2 

15 

4 

1 

-0-10125 

2-6985 

2-6083 

5-9073 

0-95868 

0-041  854 

15 

2 

15 

4 

2 

-0-082  895 

2-7257 

2-6564 

5-7324 

0-94343 

0-054  703 

15 

2 

15 

4 

4 

-0-085  998 

2-7116 

2-6387 

5-7467 

0-94308 

0-053  934 

15 

2 

15 

4 

8 

-0-085  955 

2-7125 

2-6396 

5-7466 

0-94290 

0-053  571 

15 

3 

15 

4 

1 

-0-11063 

2-7075 

2-6079 

6-0811 

0-97631 

-0-026  072 

15 

3 

15 

4 

2 

-0-095  684 

2-7389 

2-6561 

5-9385 

0-96706 

-0-023  964 

15 

3 

15 

4 

4 

-0-096643 

2-7268 

2-6432 

5-9319 

0-96006 

-0-023  619 

15 

3 

15 

4 

8 

-0-096  725 

2-7274 

2-6437 

5-9325 

0-96012 

-0-023  528 

15 

3 

15 

5 

1 

-0-104  36 

2-7114 

2-6190 

6-0227 

0-97059 

-0-028  846 

15 

3 

15 

5 

2 

-0-102  56 

2-7245 

2-6338 

6-0014 

0-96833 

-0-024  178 

15 

3 

15 

5 

4 

-0-094  744 

2-7343 

2-6529 

5-9208 

0-96194 

-0-023  842 

15 

3 

15 

5 

8 

-0-097  121 

2-7267 

2-6426 

5-9362 

0-95990 

-0-023  570 

15 

3 

15 

6 

1 

-0-104  06 

2-7148 

2-6224 

6-0145 

0-96499 

-0-027  106 

15 

3 

15 

6 

2 

-0-103  23 

2-7163 

2-6251 

6-0039 

0-96726 

-0-025  185 

15 

3 

15 

6 

4 

-0-096  636 

2-7335 

2-6497 

5-9409 

0-96343 

-0-023  819 

15 

3 

15 

6 

8 

-0-096  583 

2-7274 

2-6439 

5-9318 

0-96011 

-0-023  699 

15 

4 

15 

5 

1 

-0-106  38 

2-7056 

2-6100 

6-0371 

0-96894 

-0-024  067 

15 

4 

15 

5 

2 

-0-103  05 

2-7249 

2-6329 

6-0097 

0-96737 

-0-025  570 

15 

4 

15 

5 

4 

-0-095  427 

2-7352 

2-6526 

5-9332 

0-96108 

-0-031  900 

15 

4 

15 

5 

8 

-0-097  576 

2-7280 

2-6429 

5-9459 

0-95930 

-0-029  343 

76 


TABLE  1 — (continued) 


m 

n 

M 

N 

9 

63. 

0^2 

Oh 

O33 

O33 

15 

2 

15 

2 

1 

-0000  532  58 

0-004  866  3 

0-004  244  7 

0-018  135 

0-009  0414 

0-001  303  0 

15 

2 

15 

2 

2 

-0000  592  28 

0-004  193  7 

0-003  478  1 

0-017  756 

0-009  115  8 

0-001  474  1 

15 

2 

15 

2 

4 

-0000  587  63 

0-004  253  2 

0-003  545  4 

0-017  800 

0-009  1126 

0-001462  0 

15 

2 

15 

2 

8 

-0  000  587  55 

0-004  253  6 

0-003  546  0 

0-017  800 

0-009  112  6 

0-001  462  0 

15 

3 

15 

3 

1 

-0-000747  41 

0-003  061  3 

0-001  754  3 

0-017  798 

0-014  066 

0-004  315  3 

15 

3 

15 

3 

2 

-0-000  749  42 

0-003  109  8 

0-001810  9 

0-017  860 

0-013  810 

0-004  252  3 

15 

3 

15 

3 

4 

-0-000  740  01 

0-003  105  1 

0-001817  0 

0-017  738 

0-013  788 

0-004  251  2 

15 

3 

15 

3 

8 

-0-000  741  d4 

0-003  107  0 

0-001817  5 

0-017  756 

0-013  795 

0-004  252  5 

15 

4 

15 

4 

1 

-0-000  769  11 

0-002  632  7 

0-001  274  2 

0-017  394 

0-012  822 

0-006  553  8 

15 

4 

15 

4 

2 

-0-000  708  77 

0-003  065  0 

0-001  794  3 

0-017  260 

0-012  665 

0-006  316  3 

15 

4 

15 

4 

4 

-0-000  707  38 

0-003  016  3 

0-001  750  8 

0-017127 

0-012  543 

0-006  240  2 

15 

4 

15 

4 

8 

-0-000  708  01 

0-003  010  7 

0-001  744  0 

0-017  130 

0-012  548 

0-006  242  9 

15 

5 

15 

5 

1 

-0-000  71663 

0-002  941 1 

0-001  653  3 

0-017  273 

0-010  740 

0-007  433  1 

15 

5 

15 

5 

2 

-0-000  718  13 

0-002  934  6 

0-001650  4 

0-017  209 

0-010  931 

0-007  201  0 

15 

5 

15 

5 

4 

-0-000  702  06 

0-003  0804 

0-001  826  8 

0-017  184 

0-010  933 

0-007  019  9 

15 

5 

15 

5 

8 

-0-000  702  45 

0-003  014  8 

0-001  760  3 

0-017  076 

0-010  929 

0-006  996  0 

15 

6 

15 

6 

1 

-0-000  715  48 

0-002  983  1 

0-001  705  9 

0-017  247 

0-010436 

0-007  430  0 

15 

6 

15 

6 

2 

-0-000717  90 

0-002  984  1 

0-001700  7 

0-017  290 

0-010  667 

0-007  376  3 

15 

6 

15 

6 

4 

-0-00070613 

0-003  048  7 

0-001  787  4 

0-017  204 

0-010  761 

0-007  220  3 

15 

6 

15 

6 

8 

-0-000  702  07 

0-003  028  9 

0-001  775  5 

0-017  097 

0-010  784 

0-007  146  8 

15 

7 

15 

7 

1 

-0-000  710  60 

0-002  839  6 

0-001  573  9 

0-016903 

0-010  363 

0-007  203  2 

15 

7 

15 

7 

2 

-0-000  719  90 

0-002  968  9 

0-001  684  5 

0-017  291 

0-010565 

0-007  331  6 

15 

7 

15 

7 

4 

-0-000  711  16 

0-003  001  2 

0-001  7314 

0-017  206 

0-010  594 

0-007  272  3 

15 

7 

15 

7 

8 

-0-000  701  44 

0-003  055  0 

0-001  802  3 

0-017136 

0-010  590 

0-007  197  7 

15 

8 

15 

8 

1 

-0-000  697  83 

0-002  835  8 

0-001  591  7 

0-016  684 

0-009  960  5 

0-007  092  5 

15 

8 

15 

8 

2 

-0-000  717  94 

0-002  938  8 

0-001  658  3 

0-017  191 

0-010  416 

0-007  258  9 

15 

8 

15 

8 

4 

-0-000  713  29 

0-002  978  8 

0-001  705  2 

0-017  200 

0-010  406 

0-007  282  0 

15 

8 

15 

8 

8 

-0-000  703  39 

0-003  0510 

0-001  795  2 

0-017  162 

0-010  327 

0-007  232  7 

15 

9 

15 

9 

1 

-0-00069417 

0-002  788  0 

0-001  551  2 

0-016  526 

0-009  701  3 

0-007  078  1 

15 

9 

15 

9 

2 

-0-000  704  23 

0-002  989  3 

0-001730  4 

0-017  086 

0-010  332 

0-007  218  9 

15 

9 

15 

9 

4 

-0-000  709  97 

0-003  013  0 

0-001  744  5 

0-017  211 

0-010412 

0-007  279  8 

15 

9 

15 

9 

8 

-0-000705  05 

0-003  037  1 

0-001  778  2 

0-017  166 

0-010320 

0-007  2514 

15 

10 

15 

10 

1 

-0-000685  53 

0-002  811  6 

0-001  588  9 

0-016  437 

0-009  577  7 

0-007  043  7 

15 

10 

15 

10 

2 

-0-000  707  57 

0-002  915  6 

0-001653  8 

0-016  982 

0-010  125 

0-007  221  9 

15 

10 

15 

10 

4 

-0-000  713  28 

0-002  9864 

0-001  714  2 

0-017  204 

0-010  265 

0-007  2994 

15 

10 

15 

10 

8 

-0-000  706  86 

0-003  018  8 

0-001  756  8 

0-017  164 

0-010  241 

0-007  272  8 

15 

2 

15 

3 

1 

-0-000  638  71 

0-003  708  3 

0-002  922  5 

0-017  632 

0-009  239  5 

0-001  600  1 

15 

2 

15 

3 

2 

-0-000  530  33 

0-004  772  9 

0-004  151  8 

0-017  972 

0-009  043  1 

0-001  396  2 

15 

2 

15 

3 

4 

-0-000  582  01 

0-004  2713 

0-003  5709 

0-017  786 

0-009  117  8 

0-001  505  0 

15 

2 

15 

3 

8 

-0-000  578  53 

0-004  309  8 

0-003  615  0 

0-017  807 

0-009  114  4 

0-001  497  1 

77 


TABLE  1 — (concluded) 


m 

n 

M 

N 

<1 

d'ix 

<^32 

<5^3 

<^33 

15 

2 

15 

4 

1 

-0-000  787  12 

0-002  694  9 

0-001  694  3 

0-017  780 

0-009  4936 

0-001  801  0 

15 

2 

15 

4 

2 

-0-000  557  74 

0-004  461  1 

0-003  796  5 

0-017  825 

0-009  1015 

0-001  468  8 

15 

2 

15 

4 

4 

-0-000  570  95 

0-004  390  8 

0-003  707  7 

0-017  842 

0-009  1025 

0-001  479  5 

15 

2 

15 

4 

8 

-0-000  577  16 

0-004  322  8 

0-003  630  0 

0-017  811 

0-0091134 

0-001  495  1 

15 

3 

15 

4 

1 

-0-000  578  90 

0-004  402  9 

0-003  317  0 

0-017  892 

0-014  239 

0-004  327  4 

15 

3 

15 

4 

2 

-0-000  759  44 

0-003  095  1 

0-001783  2 

0-017  999 

0-013  917 

0-004  280  5 

15 

3 

15 

4 

4 

-0-000  744  61 

0-003  102  0 

0-001  809  2 

0-017  805 

0-013  817 

0-004  266  8 

15 

3 

15 

4 

8 

-0-000  741  43 

0-003  128  9 

0-001  840  5 

0-017  810 

0-013  821 

0-004  266  5 

15 

3 

15 

5 

1 

-0-000  517  96 

0-004  790  6 

0-003  785  4 

0-017  784 

0-014  207 

0-004  330  8 

15 

3 

15 

5 

2 

-0-000  667  89 

0-003  656  4 

0-002  457  8 

0-017  782 

0-013  999 

0-004  299  3 

15 

3 

15 

5 

4 

-0-000  785  41 

0-002  879  6 

0-001  534  3 

0-017  955 

0-013  837 

0-004  274  4 

15 

3 

15 

5 

8 

-0-000  734  17 

0-003  1694 

0-001  8904 

0-017  784 

0-013  821 

0-004  267  4 

15 

3 

15 

6 

1 

-0-000  55915 

0-004  372  9 

0-003  3116 

0-017  567 

0-014  050 

0-004  305  9 

15 

3 

15 

6 

2 

-0-000  61417 

0-004  030  5 

0-002  900  6 

0-017  729 

0-014  025 

0-004  301  4 

15 

3 

15 

6 

4 

-0-000  749  38 

0-003  127  3 

0-001  829  2 

0-017  922 

0-013  878 

0-004  279  6 

15 

3 

15 

6 

8 

-0-000  747  98 

0-003  083  9 

+0-001  786  4 

0-017  817 

0-013  822 

0-004  269  3 

15 

4 

15 

5 

1 

-0-000  925  61 

0-001  522  9 

-0-000  047  544 

0-017  482 

0-012  581 

0-006  693  6 

15 

4 

15 

5 

2 

-0-000  76152 

0-002  674  8 

+0-001  330  0 

0-017  315 

0-012  733 

0-006  394  5 

15 

4 

15 

5 

4 

-0-000  689  51 

0-003  183  7 

0-001939  6 

0-017  185 

0-012  642 

0-006  235  5 

15 

4 

15 

5 

8 

-0-000  71248 

0-002  980  5 

0-001  707  2 

0-017  134 

0-012  578 

0-006  235  6 

[ 


78 


T 


TABLE  2 

Nmnerical  Values  of  Approximatioiis  taix,  y)  to  the  Loading  on  a Tapered  Swept  Wing  of  Aspect  Ratio  2 
when  the  Control  Surface  is  Oscillating  about  its  Hinge, 
m = 15,  n = 10,  M = 15,  N = 10,  =8,  / = 1, 2, . . . , 10,  r =0-32560,  M*  = 0-78060 


h(x,y) 


V 


4 

0 00000 

0-19509 

0-38268 

0-55557 

0-70711 

0-83147 

0-92388 

0-98079 

001 

-0-00074 

-0-00189 

-0-04349 

+0-03650 

-0-07581 

-0-43460 

-0-95712 

-1-39046 

0-03 

-0-00755 

-0-00257 

-0-03097 

+0-11061 

+0-06418 

-0-14523 

-0-51628 

-0-78938 

006 

-0-00203 

-0-00715 

-0-01899 

-0-10096 

-0-21939 

-0-39160 

-0-68185 

-0-63627 

010 

+0-00095 

-0-01028 

-0-02427 

-0-19141 

-0-30860 

-0-47131 

-0-67715 

-0-45242 

015 

-0-00389 

-0-01220 

-0-04030 

-0-07286 

-0-12525 

-0-29648 

-0-44353 

-0-24766 

0-20 

-0-01010 

-0-01618 

-0-04694 

+0-01337 

-0-04112 

-0-23565 

-0-36681 

-0-20765 

0-25 

-0-01365 

-0-02281 

-0-04667 

-0-03830 

-0-16006 

-0-37652 

-0-48970 

-0-28164 

0-30 

-0-01568 

-0-03045 

-0-05097 

-0-15836 

-0-34213 

-0-56731 

-0-62918 

-0-34285 

0-35 

-0-01914 

-0-03783 

-J-06600 

-0-23743 

-0-43422 

-0-65711 

-0-64111 

-0-32812 

0-40 

-0-02594 

-0-04521 

-0-08841 

-0-22306 

-0-40129 

-0-62260 

-0-53781 

-0-27671 

0-45 

-0-03603 

-0-05410 

-0-10946 

-0-14912 

-0-33515 

-0-56491 

-0-45755 

-0-28526 

0-50 

-0-04787 

-0-06625 

-0-12268 

-0-10915 

-0-38132 

-0-62569 

-0-55382 

-0-42883 

0-55 

-0-05963 

-0-08251 

-0-12996 

-0- 19997 

-0-64113 

-0-88676 

-0-88252 

-0-69828 

0-60 

-0-07051 

-0-10232 

-0-14260 

-0-46456 

-1-10029 

-1-30704 

-1-34765 

-0-98741 

0-65 

-0-08130 

-0-12396 

-0-17616 

-0-85816 

-1-61686 

-1-72711 

-1-73382 

-1-13935 

0-68 

-0-08857 

-0-13699 

-0-21109 

-1-10536 

-1-86525 

-1-89247 

-1-82991 

-1-11053 

0-70 

-0-09407 

-0-14550 

-0-24070 

-1-25281 

-1-97644 

-1-94227 

-1-81539 

-1-03474 

0-72 

-0-10024 

-0-15377 

-0-27458 

-1-37358 

-2-03157 

-1-93528 

-1-73395 

-0.91715 

0-75 

-0-11093 

-0-16571 

-0-33037 

-1-48342 

-1-99595 

-1-81480 

-1-49554 

-0-68112 

0-80 

-0-13218 

-0-18417 

-0-41834 

-1-42519 

-1-63370 

-1-37169 

-0-90088 

-0-24695 

0-85 

-0-15446 

-0-19936 

-0-46281 

-1-07807 

-1-04704 

-0-82157 

-0-34836 

+0-01419 

0-90 

-0-16898 

-0-20407 

-0-42714 

-0-61656 

-0-53606 

-0-43508 

-0-12911 

-0-02855 

0-94 

-0-16326 

-0-18677 

-0-33586 

-0-36251 

-0-34633 

-0-30926 

-0-15724 

-0-13087 

0-97 

-0-13517 

-0-14597 

-0-23535 

-0-29872 

-0-30318 

-0-23479 

-0-12230 

-0-06467 

0-99 

-0-08648 

-0-08813 

-0-13772 

-0-24236 

-0-22438 

-0-12760 

-0-01873 

+0-05462 

h(x,y) 


V 


0-00000 

0-19509 

0-38268 

0-55557 

0-70711 

0-83147 

0-92388 

0-98079 

0-01 

0-01419 

0-03523 

0-07638 

0-17930 

0-27935 

0-37386 

0-43515 

0-40115 

0-03 

0-00714 

0-02440 

0-04585 

0-12139 

0-17180 

0-21293 

0-24798 

0-21555 

0-06 

0-00999 

0-01971 

0-04192 

0-05792 

0-08133 

0-11276 

0-13914 

0-10881 

0-10 

0-01318 

0-01935 

0-03879 

0-03238 

0-05149 

0-07843 

0-09330 

0-04887 

0-15 

0-01402 

0-02197 

0-03377 

0-04896 

0-07269 

0-08848 

0-08465 

0-02136 

0-20 

0-01488 

0-02428 

0-03326 

0-06320 

0-08265 

0-08761 

0-06718 

+0-00401 

0-25 

0-01732 

0-02573 

0-03654 

0-05840 

0-06797 

0-06633 

0-03517 

-0-01132 

0-30 

0-02075 

0-02720 

0-03966 

0-04474 

0-04758 

0-04179 

+0-00530 

-0-01997 

0-35 

0-02408 

0-02933 

0-04023 

0-03569 

0-03774 

0-02668 

-0-01161 

-0-02024 

0-40 

0-02671 

0-03194 

0-03867 

0-03661 

0-03998 

0-02041 

-0-01733 

-0-01658 

0-45 

0-02868 

0-03434 

0-03698 

0-04414 

0-04510 

0-01428 

-0-02057 

-0-01569 

0-50 

0-03043 

0-03582 

0-03681 

0-05087 

0-04240 

+0-00019 

-0-02878 

-0-02172 

0-55 

0-03230 

0-03608 

0-03825 

0-05021 

+0-02626 

-0-02396 

-0-04388 

-0-03413 

79 


TABLE  2 — (concluded) 


h(x,  y) 


V 


( 

000000 

0-19509 

0-38268 

0-55557 

0-70711 

0-83147 

0-92388 

0-98079 

0-60 

003428 

0-03532 

0-03982 

0-03884 

-0-00278 

-0-05457 

-0-06323 

-0-04890 

0-65 

003598 

0-03391 

0-03922 

+0-01652 

-0-04098 

-0-08622 

-0-08347 

-0-06159 

0-68 

003660 

0-03285 

0-03697 

-0-00162 

-0-06677 

-0-10417 

-0-09543 

-0-06712 

0-70 

0 03674 

0-03203 

0-03451 

-0-01546 

-0-08478 

-0-11552 

-0-10336 

-0-06985 

0-72 

003662 

0-03110 

0-03123 

-0-03056 

-0-10325 

-0-12636 

-0-11126 

-0-07186 

0-75 

0-03587 

0-02935 

0-02478 

-0-05515 

-0-13119 

-0-14154 

-0-12283 

-0-07357 

0-80 

0-03282 

0-02503 

+0-01028 

-0-09878 

-0-17424 

-0-16223 

-0-13852 

-0-07270 

0-85 

0-02731 

0-01820 

-0-00764 

-0-13810 

-0-20113 

-0-17022 

-0-14048 

-0-06537 

0-90 

0-01949 

0-00883 

-0-02550 

-0-15820 

-0-19454 

-0-15493 

-0-11781 

-0-05046 

0-94 

0-01210 

+0-00077 

-0-03415 

-0-14672 

-0-15865 

-0-12400 

-0-08496 

-0-03775 

0-97 

0-00646 

-0-00378 

-0-03140 

-0-11278 

-0-11470 

-0-09209 

-0-06022 

-0-03212 

0-99 

0-00281 

-0-00418 

-0-02012 

-0-06734 

-0-06932 

-0-05866 

-0-04003 

-0-02534 

80 


TABLE  3 


Numerical  Values  of  Approximatious  d„/  = l,2,3;/  = l,2,3,  to  the  Generalised  Airforces  on  a Tapered 
Wing  of  Aspect  Ratio  6,  Oscillating  in  Heave,  Pitch  and  Controi  Surface  Rotation  with  i'  = 3'15(r9, 

Af*  = 0'4 


m 

n 

M 

N 

O'u 

On 

d'i2 

Oh 

Oh 

Oh 

0^2 

Oh 

14 

6 

14 

6 

12 

37-780 

-13-417 

32-235 

-25-150 

55-368 

-16-452 

60-928 

-35-208 

14 

6 

27 

6 

12 

36-990 

-13-526 

31-155 

-25-022 

54-831 

-16-602 

59-937 

-35-190 

14 

6 

28 

6 

12 

37-148 

-13-514 

31-330 

-25-046 

54-959 

-16-582 

60-104 

-35-201 

15 

6 

15 

6 

12 

36-604 

-13-628 

30-576 

-24-988 

54-702 

-16-741 

59-529 

-35-275 

15 

6 

29 

6 

12 

37-108 

-13-606 

31-040 

-25-121 

55-015 

-16-653 

59-968 

-35-287 

15 

6 

30 

6 

12 

37-156 

-13-589 

31-107 

-25-118 

55-035 

-16-633 

60-012 

-35-273 

22 

6 

22 

6 

8 

37-456 

-13-633 

31-312 

-25-232 

55-391 

-16-652 

60-377 

-35-380 

22 

6 

43 

6 

8 

37-134 

-13-645 

30-919 

-25-166 

55-074 

-16-676 

59-926 

-35-327 

22 

6 

44 

6 

8 

37-173 

-13-642 

30-962 

-25-171 

55-105 

-16-671 

59-966 

-35-329 

23 

6 

23 

6 

8 

36-983 

-13-684 

30-708 

-25-177 

54-984 

-16-737 

59-762 

-35-381 

23 

6 

45 

6 

8 

37-168 

-13-653 

30-936 

-25-180 

55-111 

-16-682 

59-953 

-35-342 

23 

6 

46 

6 

8 

37-167 

-13-645 

30-945 

-25-173 

55-101 

-16-673 

59-952 

-35-331 

30 

6 

30 

6 

8 

37-186 

-13-646 

30-969 

-25-183 

55-128 

-16-677 

59-994 

-35-347 

31 

6 

31 

6 

8 

37-189 

-13-663 

30-952 

-25-198 

55-155 

-16-694 

60-002 

-35-369 

m 

n 

M 

N 

Oi3 

di3 

Om 

O31 

Oh 

Oh 

Oh 

Oh 

14 

6 

14 

6 

12 

-1-9962 

-0-36897 

0-39310 

-0-041884 

-3-4483 

-0-84922 

0-67051 

-0-14343 

14 

6 

27 

6 

12 

-1-9856 

-0-37948 

0-38768 

-0-040999 

-3-4323 

-0-85459 

0-66562 

-0-14118 

14 

6 

28 

6 

12 

-1-9867 

-0-37862 

0-38801 

-0-041066 

-3-4334 

-0-85371 

0-66577 

-0-14132 

15 

6 

15 

6 

12 

-1-9823 

-0-39329 

0-38081 

-0-042441 

-3-4224 

-0-86763 

0-65651 

-0-14093 

15 

6 

29 

6 

12 

-1-9693 

-0-38541 

0-38485 

-0-042582 

-3-4107 

-0-85736 

0-65844 

-0-14201 

15 

6 

30 

6 

12 

-1-9634 

-0-38517 

0-38512 

-0-042490 

-3-4100 

-0-85704 

0-65884 

-0-14194 

22 

6 

22 

6 

8 

-1-9833 

-0-39212 

0-39046 

-0-043107 

-3-4242 

-0-86718 

0-66427 

-0-14396 

22 

6 

43 

6 

8 

-1-9788 

-0-38920 

0-38894 

-0-042988 

-3-4211 

-0-86169 

0-66146 

-0-14347 

22 

6 

44 

6 

8 

-1-9792 

-0-38908 

0-38894 

-0-042987 

-3-4215 

-0-86154 

0-66140 

-0-14346 

23 

6 

23 

6 

8 

-1-9763 

-0-39093 

0-39053 

-0-043204 

-3-4188 

-0-86521 

0-66399 

-0-14407 

23 

6 

45 

6 

8 

-1-9791 

-0-38918 

0-38901 

-0-042994 

-3-4215 

-0-86153 

0-66144 

-0-14348 

23 

6 

46 

6 

8 

-1-9790 

-0-38916 

0-38894 

-0-042996 

-3-4213 

-0-86148 

0-66129 

-0-14347 

30 

6 

30 

6 

8 

-1-9790 

-0-38933 

0-38959 

-0-042961 

-3-4217 

-0-86316 

0-66308 

-0-14360 

31 

6 

31 

6 

8 

-1-9785 

-0-38957 

0-38934 

-0-C:2982 

-3-4212 

-0-86319 

0-66275 

-0-14356 

m 

N 

M 

N 

<?  OhOh 

14 

6 

14 

6 

12  -0-023448-0-035733 

14 

6 

27 

6 

12  -0-023851  -0-035654 

14 

6 

28 

6 

12  -0-023865-0-035641 

15 

6 

15 

6 

12  -0-024078-0-036449 

15 

6 

29 

6 

12  -0-023986-0-036189 

15 

6 

30 

6 

12  -0-024001  -0-036180 

81 


TABLE  3— (concluded) 


m N M N q O33O33 

22  6 22  6 8 -0-023692 -0  037150 

22  6 43  6 8 -0-023994  -0-036940 

22  6 44  6 8 -0-024002  -0-036934 

23  6 23  6 8 -0-023757-0-037180 

23  6 45  6 8 -0-024011-0-036952 

23  6 46  6 8 -0-024017-0-036946 

30  6 30  6 8 -0-023652-0-037585 

31  6 31  6 8 -0-023689-0-037647 


TABLE  4 

Numerical  Values,  Evaluated  by  Lehrian  and  Garner*,  of  Approximations  Qij,  i = l,2;  j = 1, 2,  to  the 
Generalised  Airforces  on  a Tapered  Swept  Wing  of  Aspect  Ratio  6,  Oscillating  in  Heave  and  Pitch  with 

i'=s  3-1569,  Moo  = 0-4 


m 

n 

a 

d'n 

Ola 

Ol'a 

Oi, 

O21 

O22 

Oh 

14 

6 

37-745 

-13-391 

32-129 

-25-110 

55-284 

-16-399 

60-733 

-35-123 

15 

6 

36-413 

-13-586 

30-391 

-24-871 

54-398 

-16-664 

59-203 

-35-087 

22 

6 

37-370 

-13-610  1 

31-180 

-25-169 

55-214 

-16-610 

60-094 

-35-270 

23 

6 

36-841 

-13-646 

30-545 

-25-084 

54-737 

-16-675 

59-438 

-35-230 

82 


TABLE  5 


Nninerical  Values  of  Approximations  Qij,  / = 1, 2;  / = 1, 2,  to  the  Generalised  Airforces  on  a Rectangular 
Wing  of  Aspect  Ratio  2 Oscillating  in  Heave  and  Pitch  with  i/  = 1*0,  Mco  — 0*8,  qi=q=  32, 1 = 1,2,...,N 


m 

n 

M 

N 

Oi. 

On 

On 

On 

02. 

Oh 

Oh 

^22 

4 

2 

4 

2 

0-90594 

-3-0921 

-2-9891 

-3-2292 

0-79937 

-0-77193 

-0-43632 

-1-8191 

4 

2 

9 

2 

0-90593 

-3-0920 

-2-9890 

-3-2291 

0-79937 

-0-77192 

-0-43631 

-1-8190 

4 

2 

14 

2 

0-90593 

-3-0920 

-2-9890 

-3-2292 

0-79937 

-0-77192 

-0-43631 

-1-8190 

4 

2 

19 

2 

0-90593 

-3-0920 

-2-9890 

-3-2292 

0-79937 

-0-77192 

-0-43630 

-1-8190 

4 

2 

4 

4 

0-91725 

-3-1079 

-3-0054 

-3-2545 

0-81469 

-0-77069 

-0-42610 

-1-8405 

4 

2 

9 

4 

0-91730 

-3-1079 

-3-0054 

-3-2546 

0-81475 

-0-77069 

-0-42608 

-1-8406 

4 

2 

14 

4 

0-91730 

-3-1079 

-3-0054 

-3-2546 

0-81475 

-0-77069 

-0-42608 

-1-8406 

4 

2 

19 

4 

0-91730 

-3-1079 

-3-0054 

-3-2546 

0-81474 

-0-77069 

-0-42608 

-1-8406 

4 

2 

4 

6 

0-91731 

-3-1081 

-3-0056 

-3-2547 

0-81480 

-0-77077 

-0-42615 

-1-8407 

4 

2 

9 

6 

0-91738 

-3-1081 

-3-0055 

-3-2547 

0-81486 

-0-77068 

-0-42603 

-1-8407 

4 

2 

14 

6 

0-91735 

-3-1080 

-3-0055 

-3-2547 

0-81483 

-0-77066 

-0-42602 

-1-8407 

4 

2 

19 

6 

0-91734 

-3-1080 

-3-0055 

-3-2546 

0-81481 

-0-77065 

-0-42601 

-1-8406 

4 

2 

4 

8 

0-91742 

-3-1083 

-3-0059 

-3-2550 

0-81494 

-0-77085 

-0-42619 

-1-8409 

4 

2 

9 

8 

0-91735 

-3-1080 

-3-0055 

-3-2547 

0-81483 

-0-77065 

-0-42600 

-1-8406 

4 

2 

14 

8 

0-91728 

-3-1079 

-3-0053 

-3-2545 

0-81475 

-0-77060 

-0-42597 

-1-8405 

4 

2 

19 

8 

0-91723 

-3-1078 

-3-0052 

-3-2543 

0-81470 

-0-77056 

-0-42595 

-1-8404 

4 

4 

4 

4 

0-90950 

-3-2618 

-3-3188 

-3-3228 

0-96652 

-0-84864 

-0-49919 

-2-1919 

4 

4 

9 

4 

0-90957 

-3-2618 

-3-3188 

-3-3229 

0-96660 

-0-84865 

-0-49915 

-2-1920 

4 

4 

14 

4 

0-90956 

-3-2618 

-3-3188 

-3-3229 

0-96660 

-0-84865 

-0-49915 

-2-1920 

4 

4 

19 

4 

0-90956 

-3-2618 

-3-3188 

-3-3229 

0-96660 

-0-84865 

-0-49915 

-2-1920 

4 

4 

4 

6 

0-90972 

-3-2620 

-3-3190 

-3-3232 

0-96675 

-0-84866 

-0-49911 

-2-1923 

4 

4 

9 

6 

0-90983 

-3-2620 

-3-3190 

-3-3233 

0-96681 

-0-84858 

-0-49900 

-2-1924 

4 

4 

14 

6 

0-90981 

-3-2620 

-3-3190 

-3-3233 

0-96680 

-0-84858 

-0-49901 

-2-1923 

4 

4 

19 

6 

0-90981 

-3-2620 

-3-3190 

-3-3232 

0-96679 

-0-84858 

-0-49901 

-2-1923 

4 

4 

4 

8 

0-90991 

-3-2622 

-3-3191 

-3-3236 

0-96694 

-0-84873 

-0-49909 

-2-1927 

4 

4 

9 

8 

0-90987 

-3-2620 

-3-3190 

-3-3234 

0-96685 

-0-84858 

-0-49899 

-2-1925 

4 

4 

14 

8 

0-90983 

-3-2620 

-3-3190 

-3-3233 

0-96680 

-0-84857 

-0-49899 

-2-1924 

4 

4 

19 

8 

0-90981 

-3-2620 

-3-3190 

-3-3232 

0-96679 

-0-84856 

-0-49900 

-2-1923 

4 

6 

4 

6 

0-90947 

-3-2621 

-3-3194 

-3-3228 

0-96684 

-0-84890 

-0-49963 

-2-1928 

4 

6 

9 

6 

0-90958 

-3-2621 

-3-3195 

-3-3229 

0-96689 

-0-84881 

-0-49950 

-2-1929 

4 

6 

14 

6 

0-90957 

-3-2621 

-3-3194 

-3-3228 

0-96688 

-0-84881 

-0-49950 

-2-1929 

4 

6 

19 

6 

0-90956 

-3-2621 

-3-3194 

-3-3228 

0-96687 

-0-84881 

-0-49950 

-2-1929 

4 

6 

4 

8 

0-90966 

-3-2622 

-3-3195 

-3-3232 

0-96703 

-0-84898 

-0-49961 

-2-1932 

4 

6 

9 

8 

0-90963 

-3-2621 

-3-3194 

-3-3230 

0-96694 

-0-84882 

-0-49949 

-2-1930 

4 

6 

14 

8 

0-90959 

-3-2621 

-3-3194 

-3-3229 

0-96690 

-0-84881 

-0-49949 

-2-1929 

4 

6 

19 

8 

0-90959 

-3-2621 

-3-3194 

-3-3229 

0-96689 

-0-84881 

-0-49949 

-2-1929 

4 

8 

4 

8 

0-90966 

-3-2622 

-3-3195 

-3-3232 

0-96703 

-0-84898 

-0-49961 

-2-1932 

4 

8 

9 

8 

0-90963 

-3-2621 

-3-3194 

-3-3230 

0-96694 

-0-84882 

-0-49949 

-2-1930 

4 

8 

14 

8 

0-90960 

-3-2621 

-3-3194 

-3-3229 

0-96690 

-0-84881 

-0-49949 

-2-1929 

4 

8 

19 

8 

0-90959 

-3-2621 

-3-3194 

-3-3229 

0-96689 

-0-84880 

-0-49949 

-2-1929 

83 


TABLE  5 — (continued) 


r 


m 

n 

M 

N 

On 

0t2 

Oh 

Oil 

dh 

Oh 

022 

9 

2 

9 

2 

0-90653 

-3-0924 

-2-9891 

-3-2302 

0-79973 

-0-77198 

-0-43620 

-1-8197 

9 

2 

14 

2 

0-90593 

-3-0920 

-2-9890 

-3-2292 

0-79938 

-0-77192 

-0-43631 

-1-8191 

9 

2 

19 

2 

0-90608 

-3-0921 

-2-9890 

-3-2295 

0-79946 

-0-77194 

-0-43628 

-1-8192 

9 

2 

9 

4 

0-91817 

-3-1088 

-3-0060 

-3-2562 

0-81539 

-0-77081 

-0-42596 

-1-8417 

9 

2 

14 

4 

0-91732 

-3-1079 

-3-0054 

-3-2546 

0-81477 

-0-77068 

-0-42606 

-1-8406 

9 

2 

19 

4 

0-91753 

-3-1081 

-3-0056 

-3-2550 

0-81492 

-0-77071 

-0-42604 

-1-8409 

9 

2 

9 

6 

0-91838 

-3-1091 

-3-0064 

-3-2567 

0-81567 

-0-77085 

-0-42590 

-1-8421 

9 

2 

14 

6 

0-91738 

-3-1080 

-3-0055 

-3-2547 

0-81486 

-0-77065 

-0-42600 

-1-8407 

9 

2 

19 

6 

0-91760 

-3-1083 

-3-0057 

-3-2552 

0-81504 

-0-77068 

-0-42597 

-1-8410 

9 

2 

9 

8 

0-91845 

-3-1093 

-3-0066 

-3-2569 

0-81577 

-0-77086 

-0-42588 

-1-8423 

9 

2 

14 

8 

0-91730 

-3-1079 

-3-0053 

-3-2545 

0-81478 

-0-77058 

-0-42595 

-1-8405 

9 

2 

19 

8 

0-91752 

-3-1081 

-3-0055 

-3-2549 

0-81495 

-0-77060 

-0-42590 

-1-8408 

9 

4 

9 

4 

0-91062 

-3-2623 

-3-3187 

-3-3246 

0-96719 

-0-84853 

-0-49870 

-2-1931 

9 

4 

14 

4 

0-90966 

-3-2619 

-3-3188 

-3-3230 

0-96667 

-0-84864 

-0-49911 

-2-1922 

9 

4 

19 

4 

0-90990 

-3-2620 

-3-3188 

-3-3234 

0-96680 

-0-84861 

-0-49900 

-2-1924 

9 

4 

9 

6 

0-91112 

-3-2626 

-3-3190 

-3-3254 

0-96760 

-0-84844 

-0-49846 

-2-1938 

9 

4 

14 

6 

0-90997 

-3-2620 

-3-3190 

-3-3235 

0-96692 

-0-84854 

-0-49892 

-2-1926 

9 

4 

19 

6 

0-91025 

-3-2622 

-3-3190 

-3-3239 

0-96709 

-0-84851 

-0-49881 

-2-1928 

9 

4 

9 

8 

0-91129 

-3-2628 

-3-3192 

-3-3258 

0-96776 

-0-84844 

-0-49840 

-2-1941 

9 

4 

14 

8 

0-90999 

-3-2620 

-3-3190 

-3-3235 

0-96693 

-0-84853 

-0-49891 

-2-1926 

9 

4 

19 

8 

0-91029 

-3-2622 

-3-3190 

-3-3240 

0-96711 

-0-84849 

-0-49878 

-2-1929 

9 

6 

9 

6 

0-91086 

-3-2627 

-3-3195 

-3-3250 

0-96766 

-0-84869 

-0-49896 

-2-1943 

9 

6 

14 

6 

0-90972 

-3-2621 

-3-3194 

-3-3231 

0-96699 

-0-84878 

-0-49941 

-2-1931 

9 

6 

19 

6 

0-91000 

-3-2623 

-3-3194 

-3-3236 

0-96715 

-0-84876 

-0-49930 

-2-1934 

9 

6 

9 

8 

0-91106 

-3-2629 

-3-3196 

-3-3254 

0-96784 

-0-84870 

-0-49891 

-2-1946 

9 

6 

14 

8 

0-90977 

-3-2622 

-3-3194 

-3-3232 

0-96703 

-0-84878 

-0-49940 

-2-1932 

9 

6 

19 

8 

0-91007 

-3-2623 

-3-3195 

-3-3237 

0-96722 

-0-84875 

-0-49928 

-2-1935 

9 

8 

9 

8 

0-91106 

-3-2629 

-3-3196 

-3-3254 

0-96784 

-0-84869 

-0-49890 

-2-1946 

9 

8 

14 

8 

0-90977 

-3-2622 

-3-3194 

-3-3232 

0-96703 

-0-84877 

-0-49939 

-2-1932 

9 

8 

19 

8 

0-91007 

-3-2623 

-3-3194 

-3-3237 

0-96721 

-0-84875 

-0-49926 

-2-1935 

14 

2 

14 

2 

0-90593 

-3-0920 

-2-9890 

-3-2292 

0-79938 

-0-77192 

-0-43631 

-1-8191 

14 

2 

19 

2 

0-90593 

-3-0920 

-2-9890 

-3-2292 

0-79938 

-0-77192 

-0-43631 

-1-8191 

14 

2 

14 

4 

0-91732 

-3-1079 

-3-0054 

-3-2546 

0-81477 

-0-77068 

-0-42606 

-1-8406 

14 

2 

19 

4 

0-91732 

-3-1079 

-3-0054 

-3-2546 

0-81477 

-0-77068 

-0-42606 

-1-8406 

14 

2 

14 

6 

0-91738 

-3-1080 

-3-0055 

-3-2547 

0-81486 

-0-77065 

-0-42600 

-1-8407 

14 

2 

19 

6 

0-91736 

-3-1080 

-3-0055 

-3-2547 

0-81484 

-0-77064 

-0-42599 

-1-8407 

14 

2 

14 

8 

0-91730 

-3-1079 

-3-0053 

-3-2545 

0-81478 

-0-77058 

-0-42595 

-1-8405 

14 

2 

19 

8 

0-91725 

-3-1078 

-3-0053 

-3-2544 

0-81473 

-0-77054 

-0-42593 

-1-8404 

84 


TABLE  5 — (concluded) 


m 

n 

M 

N 

dll 

dll 

dn 

di2 

dii 

dSi 

dzz 

dh 

14 

4 

14 

4 

0-90966 

-3-2619 

-3-3188 

-3-3230 

0-96667 

-0-84864 

-0-49911 

-2-1922 

14 

4 

19 

4 

0-90966 

-3-2619 

-3-3188 

-3-3230 

0-96667 

-0-84864 

-0-49911 

-2-1922 

14 

4 

14 

6 

0-90997 

-3-2620 

-3-3190 

-3-3235 

0-96692 

-0-84854 

-0-49892 

-2-1926 

14 

4 

19 

6 

0-90996 

-3-2620 

-3-3190 

-3-3235 

0-96692 

-0-84854 

-0-49892 

-2-1925 

14 

4 

14 

8 

0-90999 

-3-2620 

-3-3190 

-3-3235 

0-96693 

-0-84853 

-0-49891 

-2-1926 

14 

4 

19 

8 

0-90997 

-3-2620 

-3-3190 

-3-3235 

0-96692 

-0-84852 

-0-49891 

-2-1925 

14 

6 

14 

6 

0-90972 

-3-2621 

-3-3194 

-3-3231 

0-96699 

-0-84878 

-0-49941 

-2-1931 

14 

6 

19 

6 

0-90972 

-3-2621 

-3-3194 

-3-3231 

0-96699 

-0-84878 

-0-49941 

-2-1931 

14 

6 

14 

8 

0-90977 

-3-2622 

-3-3194 

-3-3232 

0-96703 

-0-84878 

-0-49940 

-2-1932 

14 

6 

19 

8 

0-90976 

-3-2622 

-3-3194 

-3-3231 

0-96703 

-0-84878 

-0-49940 

-2-1931 

14 

8 

14 

8 

0-90977 

-3-2622 

-3-3194 

-3-3232 

0-96703 

-0-84877 

-0-49939 

-2-1932 

14 

8 

19 

8 

0-90976 

-3-2621 

-3-3194 

-3-3231 

0-96702 

-0-84877 

-0-49939 

-2-1931 

19 

2 

19 

2 

0-90608 

-3-0921 

-2-9890 

-3-2295 

0-79946 

-0-77194 

-0-43628 

-1-8192 

19 

2 

19 

4 

0-91753 

-3-1081 

-3-0056 

-3-2550 

0-81492 

-0-77071 

-0-42604 

-1-8409 

19 

2 

19 

6 

0-91760 

-3-1083 

-3-0057 

-3-2552 

0-81504 

-0-77068 

-0-42597 

-1-8410 

19 

2 

19 

8 

0-91752 

-3-1081 

-3-0055 

-3-2549 

0-81495 

-0-77060 

-0-42590 

-1-8408 

19 

4 

19 

4 

0-90990 

-3-2620 

-3-3188 

-3-3234 

0-96680 

-0-84861 

-0-49900 

-2-1924 

19 

4 

19 

6 

0-91025 

-3-2622 

-3-3190 

-3-3239 

0-96709 

-0-84851 

-0-49881 

-2-1928 

19 

4 

19 

8 

0-91025 

-3-2622 

-3-3190 

-3-3240 

0-96710 

-0-84851 

-0-49882 

-2-1929 

19 

6 

19 

6 

0-91000 

-3-2623 

-3-3194 

-3-3236 

0-96715 

-0-84876 

-0-49930 

-2-1934 

19 

6 

19 

8 

0-91007 

-3-2623 

-3-3195 

-3-3237 

0-96722 

-0-84875 

-0-49928 

-2-1935 

19 

8 

19 

8 

0-91007 

-3-2623 

-3-3194 

-3-3237 

0-96721 

-0-84875 

-0-49926 

-2-1935 
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TABLE  6 


Numerical  Values  of  Approximations  Qif,  i = 1, 2;  / = 1, 2,  to  the  Generalised  Airforces  on  a Rectangular 
Wing  of  Aspect  Ratio  8 Oscillating  in  Heave  and  Pitch  with  v = 1-0,  Afoo  = 0*8,  qi  =q  = 32, 1 = 1,2,.  ..,N 


m 

n 

M 

N 

O'n 

6?. 

d'i2 

Oh 

Oil 

dh 

Oh 

dh 

4 

2 

4 

2 

-1-2372 

-15-939 

-19-123 

-9-5648 

1-9318 

-5-1566 

-5-0567 

-7-3567 

4 

2 

9 

2 

-1-2323 

-15-956 

-19-143 

-9-5824 

1-9402 

-5-1633 

-5-0615 

-7-3746 

4 

2 

14 

2 

-1-2323 

-15-956 

-19-143 

-9-5824 

1-9402 

-5-1632 

-5-0614 

-7-3746 

4 

2 

19 

2 

-1-2322 

-15-956 

-19-142 

-9-5824 

1-9402 

-5-1632 

-5-0614 

-7-3746 

4 

2 

4 

4 

-1-1413 

-15-986 

-19-113 

-9-7481 

1-9823 

-5-0979 

-4-9262 

-7-3833 

4 

2 

9 

4 

-1-1294 

-16-001 

-19-125 

-9-7725 

1-9939 

-5-0996 

-4-9226 

-7-4018 

4 

2 

14 

4 

-1-1289 

-16-000 

-19-123 

-9-7724 

1-9938 

-5-0989 

-4-9215 

-7-4011 

4 

2 

19 

4 

-1-1288 

-16-000 

-19-123 

-9-7724 

1-9938 

-5-0987 

-4-9213 

-7-4009 

4 

2 

4 

6 

-1-1416 

-15-991 

-19-120 

-9-7490 

1-9859 

-5-0979 

-4-9246 

-7-3876 

4 

2 

9 

6 

-1-1308 

-16-005 

-19-132 

-9-7733 

1-9946 

-5-1018 

-4-9255 

-7-4048 

4 

2 

14 

6 

-1-1295 

-16-002 

-19-126 

-9-7729 

1-9941 

-5-0999 

-4-9228 

-7-4024 

4 

2 

19 

6 

-1-1290 

-16-001 

-19-124 

-9-7727 

1-9940 

-5-0991 

-4-9218 

-7-4016 

4 

2 

4 

8 

-1-1363 

-15-982 

-19-100 

-9-7464 

1-9870 

-5-0885 

-4-9105 

-7-3798 

4 

2 

9 

8 

-1-1312 

-16-007 

-19-135 

-9-7735 

1-9958 

-5-1022 

-4-9255 

-7-4065 

4 

2 

14 

8 

-1-1301 

-16-003 

-19-129 

-9-7732 

1-9946 

-5-1009 

-4-9240 

-7-4039 

4 

2 

19 

8 

-1-1294 

-16-002 

-19-126 

-9-7728 

1-9942 

-5-0998 

-4-9226 

-7-4024 

4 

4 

4 

4 

-1-9903 

-16-192 

-20-312 

-8-3273 

2-1285 

-5-8804 

-6-2842 

-8-3902 

4 

4 

9 

4 

-1-9848 

-16-204 

-20-325 

-8-3398 

2-1356 

-5-8812 

-6-2822 

-8-4013 

4 

4 

14 

4 

-1-9849 

-16-203 

-20-324 

-8-3387 

2-1352 

-5-8804 

-6-2812 

-8-3998 

4 

4 

19 

4 

-1-9849 

-16-203 

-20-323 

-8-3384 

2-1351 

-5-8802 

-6-2809 

-8-3395 

4 

4 

4 

6 

-1-9886 

-16-196 

-20-316 

-8-3284 

2-1317 

-5-8776 

-6-2799 

-8-3924 

4 

4 

9 

6 

-1-9838 

-16-208 

-20-331 

-8-3452 

2-1377 

-5-8828 

-6-2828 

-8-4067 

4 

4 

14 

6 

-1-9840 

-16-206 

-20-327 

-8-3426 

2-1367 

-5-8808 

-6-2804 

-8-4029 

4 

4 

19 

6 

-1-9841 

-16-205 

-20-326 

-8-3415 

2-1363 

-5-8800 

-6-2795 

-8-4015 

4 

4 

4 

8 

-1-9886 

-16-187 

-20-303 

-8-3155 

2-1295 

-5-8654 

-6-2657 

-8-3761 

4 

4 

9 

8 

-1-9834 

-16-210 

-20-332 

-8-3453 

2-1387 

-5-8824 

-6-2823 

-8-4079 

4 

4 

14 

8 

-1-9838 

-16-207 

-20-329 

-8-3440 

2-1374 

-5-8818 

-6-2815 

-8-4051 

4 

4 

19 

8 

-1-9840 

-16-206 

-20-327 

-8-3426 

2-1368 

-5-8808 

-6-2804 

-8-4031 

4 

6 

4 

6 

-1-9935 

-16-189 

-20-310 

-8-3134 

2-1268 

-5-8806 

-6-2905 

-8-3868 

4 

6 

9 

6 

-1-9884 

-16-202 

-20-325 

-8-3306 

2-1330 

-5-8863 

-6-2941 

-8-4019 

4 

6 

14 

6 

-1-9885 

-16-200 

-20-321 

-8-3284 

2-1322 

-5-8845 

-6-2919 

-8-3986 

4 

6 

19 

6 

-1-9855 

-16-199 

-20-320 

-8-3275 

2-1319 

-5-8838 

-6-2910 

-8-3974 

4 

6 

4 

8 

-1-9929 

-16-182 

-20-298 

-8-3024 

2-1254 

-5-8696 

-6-2772 

-8-3728 

4 

6 

9 

8 

-1-9880 

-16-204 

-20-326 

-8-3309 

2-1340 

-5-8859 

-6-2936 

-8-4031 

4 

6 

14 

8 

-1-9883 

-16-202 

-20-323 

-8-3297 

2-1328 

-5-8854 

-6-2930 

-8-4006 

4 

6 

19 

8 

-1-9884 

-16-200 

-20-321 

-8-3285 

2-1323 

-5-8846 

-6-2920 

-8-3988 

4 

8 

4 

8 

-1-9929 

-16-182 

-20-298 

-8-3024 

2-1254 

-5-8696 

-6-2772 

-8-3728 

4 

8 

9 

8 

-1-9880 

-16-204 

-20-326 

-8-3309 

2-1340 

-5-8860 

-6-2937 

-8-4032 

4 

8 

14 

8 

-1-9883 

-16-202 

-20-323 

-8-3297 

2-1328 

-5-8854 

-6-2930 

-8-4007 

4 

8 

19 

8 

-1-9884 

-16-200 

-20-321 

-8-3285 

2-1323 

-5-8846 

-6-2920 

-8-3989 
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TABLE  6 — (continued) 


1 


m 

n 

M 

N 

On 

Q'i2 

Ol2 

021 

d'22 

022 

9 

2 

9 

2 

-1-2419 

-15-950 

-19-138 

-9-5696 

1-9332 

-5-1624 

-5-0659 

-7-3735 

9 

2 

14 

2 

-1-2429 

-15-938 

-19-126 

-9-5601 

1-9311 

-5-1573 

-5-0601 

-7-3650 

9 

2 

19 

2 

-1-2426 

-15-941 

-19-129 

-9-5625 

1-9317 

-5-1585 

-5-0615 

-7-3671 

9 

2 

9 

4 

-1-1379 

-15-998 

-19-128 

-9-7635 

1-9904 

-5-0986 

-4-9264 

-7-4075 

9 

2 

14 

4 

-1-1355 

-15-980 

-19-103 

-9-7545 

1-9869 

-5-0901 

-4-9160 

-7-3939 

9 

2 

19 

4 

-1-1359 

-15-984 

-19-108 

-9-7567 

1-9878 

-5-0918 

-4-9181 

-7-3969 

9 

2 

9 

6 

-1-1413 

-16-006 

-19-141 

-9-7637 

1-9920 

-5-1028 

-4-9320 

-7-4133 

9 

2 

14 

6 

-1-1360 

-15-982 

-19-106 

-9-7550 

1-9873 

-5-0911 

-4-9173 

-7-3953 

9 

2 

19 

6 

-1-1365 

-15-986 

-19-111 

-9-7568 

1-9882 

-5-0927 

-4-9192 

-7-3983 

9 

2 

9 

8 

-1-1433 

-16-010 

-19-149 

-9-7634 

1-9938 

-5-1046 

-4-9339 

-7-4170 

9 

2 

14 

8 

-1-1367 

-15-983 

-19-109 

-9-7552 

1-9877 

-5-0921 

-4-9185 

-7-3967 

9 

2 

19 

8 

-1-1373 

-15-987 

-19-114 

-9-7568 

1-9885 

-5-0937 

-4-9205 

-7-3996 

9 

4 

9 

4 

-2-0062 

-16-200 

-20-327 

-8-3118 

2-1204 

-5-8855 

-6-2962 

-8-3944 

9 

4 

14 

4 

-2-0091 

-16-184 

-20-310 

-8-2965 

2-1165 

-5-8790 

-6-2889 

-8-3802 

9 

4 

19 

4 

-2-0084 

-16-188 

-20-313 

-8-2998 

2-1173 

-5-8802 

-6-2902 

-8-3831 

9 

4 

9 

6 

-2-0050 

-16-207 

-20-337 

-8-3187 

2-1236 

-5-8886 

-6-2991 

-8-4031 

9 

4 

14 

6 

-2-0080 

-16-187 

-20-313 

-8-3006 

2-1180 

-5-8793 

-6-2880 

-8-3835 

9 

4 

19 

6 

-2-0074 

-16-190 

-20-316 

-8-3035 

2-1188 

-5-8803 

-6-2892 

-8-3860 

9 

4 

9 

8 

-2-0047 

-16-211 

-20-341 

-8-3198 

2-1253 

-5-8894 

-6-3003 

-8-4068 

9 

4 

14 

8 

-2-0078 

-16-188 

-20-315 

-8-3020 

2-1187 

-5-8803 

-6-2.'’92 

-8-3856 

9 

4 

19 

8 

-2-0073 

-16-192 

-20-319 

-8-3048 

2-1194 

-5-8814 

-6-2905 

-8-3882 

9 

6 

9 

6 

-2-0099 

-16-201 

-20-331 

-8-3038 

2-1187 

-5-8921 

-6-3103 

-8-3980 

9 

6 

14 

6 

-2-0126 

-16-181 

-20-307 

-8-2863 

2-1135 

-5-8831 

-6-2997 

-8-3793 

9 

6 

19 

6 

-2-0120 

-16-184 

-20-311 

-8-2893 

2-1143 

-5-8842 

-6-3008 

-8-3819 

9 

6 

9 

8 

-2-0097 

-16-205 

-20-336 

-8-3050 

2-1205 

-5-8928 

-6-3113 

-8-4016 

9 

6 

14 

8 

-2-0125 

-16-182 

-20-309 

-8-2876 

2-1142 

-5-8840 

-6-3007 

-8-3813 

9 

6 

19 

8 

-2-0119 

-16-186 

-20-313 

-8-2905 

2-1149 

-5-8852 

-6-3021 

-8-3839 

9 

8 

9 

8 

-2-0097 

-16-205 

-20-336 

-8-3051 

2-1205 

-5-8929 

-6-3114 

-8-4017 

9 

8 

14 

8 

-2-0125 

-16-182 

-20-309 

-8-2876 

2-1142 

-5-8841 

-6-3008 

-8-3813 

9 

8 

19 

8 

-2-0119 

-16-186 

-20-313 

-8-2905 

2-1149 

-5-8852 

-6-3021 

-8-3840 

14 

2 

14 

2 

-1-2429 

-15-938 

-19-126 

-9-5601 

1-9311 

-5-1573 

-5-0601 

-7-3650 

14 

2 

19 

2 

-1-2429 

-15-938 

-19-126 

-9-5601 

1-9312 

-5-1573 

-5-0601 

-7-3650 

14 

2 

14 

4 

-1-1355 

-15-980 

-19-103 

-9-7545 

1-9869 

-5-0901 

-4-9160 

-7-3939 

14 

2 

19 

4 

-1-1354 

-15-980 

-19-102 

-9-7544 

1-9869 

-5-0899 

-4-9157 

-7-3937 

14 

2 

14 

6 

-1-1360 

-15-982 

-19-106 

-9-7550 

1-9873 

-5-0911 

-4-9173 

-7-3953 

14 

2 

19 

6 

-1-1356 

-15-980 

-19-103 

-9-7547 

1-9871 

-5-0903 

-4-9162 

-7-3944 

14 

2 

14 

8 

-1-1367 

-15-983 

-19-109 

-9-7552 

1-9877 

-5-0921 

-4-9185 

-7-3967 

14 

2 

19 

8 

-1-1360 

-15-981 

-19-105 

-9-7549 

1-9873 

-5-0910 

-4-9171 

-7-3952 

TABLE  6 — (concluded) 


-2  0091 
-20091 

-20080 
-2  0081 

-2-0078 
-2  0079 

-2-0126 

-2-0126 

-2-0124 

-2-0125 

-2-0124 

-2-0125 

-1-2426 

-1-1358 

-1-1365 

-1-1373 

-2-0084 

-2-0074 

-2-0072 

-2-0119 

-2-0118 


n 

M 

N 

4 

14 

4 

4 

19 

4 

4 

14 

6 

4 

19 

6 
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14 

8 
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19 

8 

6 

14 

6 

6 

19 

6 

6 

14 

8 

6 

19 

8 

8 

14 

8 

8 

19 

8 
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19 

2 

2 

19 

4 

2 

19 

6 

2 

19 

8 

4 

19 

4 

4 

19 

6 

19  4 19  8 - 

19  6 19  6 - 

19  6 19  8 - 


19  8 19  8 -2-0118  - 


dii 

&i2 

&I2 

O21 

16-184 

-20-310 

-8-2965 

2-1165 

16-184 

-20-309 

-8-2963 

2-1164 

16-187 

-20-313 

-8-3007 

2-1181 

16-186 

-20-312 

-8-2997 

2-1177 

16-188 

-20-315 

-8-3021 

2-1188 

16-187 

-20-313 

-8-3008 

2-1181 

16-181 

-20-307 

-8-2864 

2-1135 

16-180 

-20-306 

-8-2855 

2-1133 

16-182 

-20-309 

-8-2877 

2-1142 

16-181 

-20-307 

-8-2865 

2-1137 

16-182 

-20-309 

-8-2877 

2-1142 

16-181 

-20-307 

-8-2865 

2-1137 

15-941 

-19-129 

-9-5625 

1-9317 

15-984 

-19-108 

-9-7567 

1-9878 

15-986 

-19-111 

-9-7568 

1-9882 

15-987 

-19-114 

-9-7569 

1-9885 

16-188 

-20-313 

-8-2998 

2-1173 

16-190 

-20-316 

-8-3035 

2-1189 

16-192 

-20-319 

-8-3049 

2-1195 

16-184 

-20-311 

-8-2893 

2-1144 

16-186 

-20-313 

-8-2906 

2-1149 

16-186 

-20-313 

-8-2906 

2-1149 

-5-8789 

-5-8788 

-5-8793 

-5-8785 

-5-8803 

-5-8793 

-5-8831 

-5-8824 

-5-8840 

-5-8832 

-5-8840 

-5-8832 

-5-1585 

-5-0918 

-5-0927 

-5-0937 

-5-8802 

-5-8803 

-5-8814 

-5-8814 

-5-8852 


-6-2889 

-6-2887 

-6-2880 

-6-2871 

-6-2892 

-6-2880 

-6-2997 

-6-2988 

-6-3007 

-6-2997 

-6-3008 

-6-2997 

-5-0615 

-4-9181 

-4-9192 

-4-9205 

-6-2902 

-6-2892 

-6-2905 

-6-3008 

-6-3021 


-8-3803 

-8-3799 

-8-3835 

-8-3821 

-8-3857 

-8-3837 

-8-3794 

-8-3781 

-8-3814 

-8-3796 

-8-3814 

-8-3796 

-7-3671 

-7-3969 

-7-3938 

-7-3996 

-8-3831 

-8-3860 

-8-3882 

-8-3820 

-8-3840 


-5-8852  -6-3021  -8-3840 
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1 
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0 000 

0000 
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014 

017 

0-21 

0-25 
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012 
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004 
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on 

017 
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0000 

0-000 

012 

0-09 
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0-18 

012 

0-12 

0001 

0-001 
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0-023 

0031 

0-033 
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Aq  = ir/3C 
s = c 

Co  = { (3  + 2 /T)c  = 1-6160254  c 

C,  = -1(5  -2  /l)  c = 0-3839746  c 
e,  s v/Tc  = 1-7320508  c 

X,  = I (3-^2  = 1-61602  54  c 


X2  = ^^9-^4/^)c  = 1-9910254  c 
Y2  = c 

Yp  = sin  (ir/IO)  c = 0-1950903c 
Fio.  2.  Tapered  swept  wing  of  aspect  ratio  2. 
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L 


X 


s = 3 c 
Co  = 3c/2 
Cl  = c/2 

e,  = (/3  +1/2)c  = 2-232051  c 

Yi  = 6c/5 

Yj  = 21  c/10 

j Vr  = s sin  = 0-1950903s 

1 

[ Fig.  3.  Tapered  swept  wing  of  aspect  ratio  6. 
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100 
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0-32560, 


0 00080 


Fig.  5m.  Generalised  airforce  coefficients  on  a tapered  sv/ept  wing  of  aspect  ratio  2.  */  = 0-32560, 

Af  = 0-78060. 


Fio.  5n.  Generalised  airforce  coefficients  on  a tapered  swept  wing  of  aspect  ratio  2.  »<  = 0-32560,  M 

0-78060. 


Generalised  airforce  coefficients  on  a tapered  swept  wing  of  aspect  ratio  2.  j'  = 0"32560,  M 

0-78060. 
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taking  more  upwash  points  than  the  number  of  basic  loading  functions  taking  more  upwash  points  than  the  number  of  basic  loading  functions 

more  accurate  values  should  be  obtained  for  the  generalised  airforce  more  accurate  values  should  be  obtained  for  the  generalised  airforce 

coefficients  for  a given  number  of  basic  functions.  coefficients  for  a given  number  of  basic  functions. 

The  procedure  has  been  programmed  in  ICL 1900  FORTRAN.  The  The  procedure  has  been  programmed  in  ICL  1900  FORTRAN.  The 

program  has  been  used  to  obtain  results  for  tapered  swept  and  program  has  been  used  to  obtain  results  for  tapered  swept  and 


